Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 343 records · Page 19

A spectroscopic, structural, and computational study of Ag–oxo interactions in Ag + /UO 2 2+ complexes

For this work, twelve novel Ag + /UO 2 2+ heterometallic complexes have been prepared and characterized via structural, spectroscopic, and computational methods to probe the effects of Ag–oxo interactions on bonding and photophysical properties of the uranyl cation. Structural characterization reveals Ag–oxo interaction distances ranging from 2.475(3) Å to 4.287(4) Å. These interactions were probed using luminescence and Raman spectroscopy which displayed little effect on the luminescence intensity and the energy of the Raman active U═O symmetric stretch peak as compared to previously reported Pb–oxo interactions. Computational efforts via density functional theory-based natural bond orbital analysis revealed that the highest stabilization energy associated with the Ag–oxo interaction had a value of only 11.03 kcal mol –1 and that all other energy values fell at 7.05 kcal mol –1 or below indicating weaker interactions relative to those previously reported for Pb 2+ /UO 2 2+ heterometallic compounds. In contrast, quantum theory of atoms in molecules analysis of bond critical point electron density values indicated higher electron density in Ag–oxo interactions as compared to Pb–oxo interactions which suggests more covalent character with the Ag + . Overall, this data indicates that Ag + has a less significant effect on UO 2 2+ bonding and photophysical properties as compared to other Pb 2+ , likely due to the high polarizability of the cation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning reveals factors that control ion mobility in anti-perovskite solid electrolytes

Solid-state batteries are projected to exhibit improved energy densities and safety compared to liquid-electrolyte-based systems. Consequently, the development of solid electrolytes that can conduct cations at rates comparable to that of liquid electrolytes is an active area of research. To facilitate this search it is desirable to understand the chemical and structural features that control ion transport through a crystalline lattice. Here, machine learning (ML) is used to identify features that influence ion mobility in solids, quantify the relative importance of those features, and demonstrate how the variation of individual features alters ion mobility. Adopting the anti-perovskite lattice as a model system, ML algorithms were trained to predict migration barriers from a data set containing hundreds of barriers calculated using density functional theory. In this study, more than 10 6 feature sets of non-redundant descriptors were evaluated for their ability to predict migration barriers. The present analysis correctly reiterates the importance of several features that are known to influence mobility – for example, anion polarizability – but also identifies others whose importance is not widely recognized. Lattice properties such as hopping distance and channel width exhibit the greatest influence on cation mobility. Individual conditional expectation analyses performed on these features shows that barriers decrease as hopping distance decreases and as channel width increases. The defect formation energy – which is commonly associated with the concentration of defects, but not with their mobility – is identified as an important feature for interstitial migration. In sum, this analysis aids in the design of optimal solid electrolytes by simplifying the multi-dimensional design space to a sub-set of properties that are the most important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural mapping and tuning of mixed halide ions in amorphous sulfides for fast Li-ion conduction and high deformability

Amorphous sulfides are among the most promising candidates for solid electrolytes (SEs) owing to their excellent deformability and acceptable Li-ion conductivity (σion) at room temperature. However, the complex atomic structure of these amorphous materials without long-range ordering results in a lack of structural understanding and difficulties in tuning material properties. Here, in this study, we performed structural mapping of glassy sulfide SEs composed of various anion clusters using a combinatorial atomic level analysis of synchrotron X-ray-based pair distribution function (PDF) and reverse Monte-Carlo (RMC) methods, demonstrating the potential to tune σ ion in glass SEs. First, we prepared a new glassy sulfide with a mixed anion framework of two halogens (Br and I) as atomic anions, with PS 4 3- molecular anions in a Li 2 S–P 2 S 5 -based glass SE. At a specific Br content, [(Li 2 S) 0.658 (LiI 0.9 LiBr 0.1 ) 0.342 ] 0.825 [P 2 S 5 ] 0.175 recorded a σ ion of 2.27 mS cm -1 , the highest value for any glass sulfides reported to date; however, its elastic modulus was still suppressed to 14.48 GPa. PDF and RMC calculations successfully provided structural mapping of anion clusters, including two halogens. Molecular dynamics simulations of each composition confirmed that flexible coordination caused by the rattling of small polarizable Br ions in the mixed halogens of glassy SEs contributed to the superior σ ion . Our results may provide new insights into the design of superior glassy SEs that play key roles in all-solid-state batteries requiring fast Li-ion conduction and high deformability.

36 MATERIALS SCIENCE↗

Improved quantum yield in geometrically constrained tetraphenylethylene-based metal–organic frameworks

In this paper we report the synthesis of a novel, tetraphenylethylene-based ligand for metal–organic frameworks (MOFs). Incorporation of this ligand into a Zn- or Eu-based MOF increased the quantum yield (QY) by almost 2.5× compared to the linker alone. Furthermore, the choice of guest solvent impacted the QY and solvatochromatic response. These shifts are consistent with solvent dielectric constant as well as molecular polarizability.

36 MATERIALS SCIENCE↗

Structural chemistry of penta- and hexanitrato thorium($\tiny{IV}$) complexes isolated using N–H donors

In this study, a series of fifteen tetravalent thorium phases were prepared. The compounds were isolated from acidic aqueous nitrate solutions using protonated nitrogen heterocycles of varying hydrogen-bond donation strength. Structural analysis via single crystal X-ray diffraction showed that the structures are built from pentanitrato, [Th(NO 3 ) 5 (H 2 O) 2 ] 1- , and hexanitrato, [Th(NO 3 ) 6 ] 2- , molecular units, with the latter being far more prevalent in the solid state. The vibrational properties of the compounds were examined using Raman and IR spectroscopy; the spectra are dominated by stretches characteristic of nitrate and the organic ions. The relative energetics of nitrate complexation was examined using electronic structure theory. These results confirmed that there are clear thermodynamic sinks for the penta- and hexanitrato structural units that were observed experimentally. Additionally, electrostatic surface potentials (ESPs) were calculated in an effort to better understand the counterion stabilization of the complexes. The ESP surfaces showed that the position of the water and nitrate molecules and the coordination geometry of the metal complex had a clear effect on the polarizability of the two structural motifs. Despite limited speciation of the Th–nitrate structural units, the compounds exhibit rich supramolecular chemistry resulting from hydrogen bonding of the Th complexes with the organic N–H donors and π–π stacking interactions from the protonated N-heterocycles.

36 MATERIALS SCIENCE↗

Magnetic coupling between Fe(NO) spin probe ligands through diamagnetic Ni II , Pd II and Pt II tetrathiolate bridges

Reaction of the nitrosylated-iron metallodithiolate ligand, paramagnetic (NO)Fe(N 2 S 2 ), with [M(CH 3 CN) n ][BF 4 ] 2 salts (M = Ni II , Pd II , and Pt II ; n = 4 or 6) affords di-radical tri-metallic complexes in a stairstep type arrangement ([FeMFe] 2+ , M = Ni, Pd, and Pt), with the central group 10 metal held in a MS 4 square plane. These isostructural compounds have nearly identical ν(NO) stretching values, isomer shifts, and electrochemical properties, but vary in their magnetic properties. Despite the intramolecular Fe∙∙∙Fe distances of ca. 6 Å, antiferromagnetic coupling is observed between {Fe(NO)} 7 units as established by magnetic susceptibility, EPR, and DFT studies. The superexchange interaction through the thiolate sulfur and central metal atoms is on the order of Ni II < Pd II $\ll$ Pt II with exchange coupling constants (J) of -3, -23, and -124 cm -1 , consistent with increased covalency of the M–S bonds (3d < 4d < 5d). This trend is reproduced by DFT calculations with molecular orbital analysis providing insight into the origin of the enhancement in the exchange interaction. Specifically, the magnitude of the exchange interaction correlates surprisingly well with the energy difference between the HOMO and HOMO-1 orbitals of the triplet states, which is reflected in the central metal's contribution to these orbitals. These results demonstrate the ability of sulfur-dense metallodithiolate ligands to engender strong magnetic communication by virtue of their enhanced covalency and polarizability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

“Freedom of design” in chemical compound space: towards rational in silico design of molecules with targeted quantum-mechanical properties

The rational design of molecules with targeted quantum-mechanical (QM) properties requires an advanced understanding of the structure–property/property–property relationships (SPR/PPR) that exist across chemical compound space (CCS). In this work, we analyze these fundamental relationships in the sector of CCS spanned by small (primarily organic) molecules using the recently developed QM7-X dataset, a systematic, extensive, and tightly converged collection of 42 QM properties corresponding to ≈4.2M equilibrium and non-equilibrium molecular structures containing up to seven heavy/non-hydrogen atoms (including C, N, O, S, and Cl). By characterizing and enumerating progressively more complex manifolds of molecular property space—the corresponding high-dimensional space defined by the properties of each molecule in this sector of CCS—our analysis reveals that one has a substantial degree of flexibility or “freedom of design” when searching for a single molecule with a desired pair of properties or a set of distinct molecules sharing an array of properties. To explore how this intrinsic flexibility manifests in the molecular design process, we used multi-objective optimization to search for molecules with simultaneously large polarizabilities and HOMO–LUMO gaps; analysis of the resulting Pareto fronts identified non-trivial paths through CCS consisting of sequential structural and/or compositional changes that yield molecules with optimal combinations of these properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

K–Co–Mo–S x chalcogel: high-capacity removal of Pb 2+ and Ag + and the underlying mechanisms

Chalcogenide-based aerogels, known as chalcogels, represent a novel class of nanoparticle-based porous amorphous materials characterized by high surface polarizability and Lewis base properties, exhibiting promising applications in clean energy and separation science. This work presents a K–Co–Mo–S x (KCMS) chalcogel as a highly efficient sorbent for heavy metal ions and details its sorption mechanisms. Its incoherent structure comprises Mo 2 V (S 2 ) 6 and Mo 3 IV S(S 6 ) 2 anion-like clusters with four- and six-coordinated Co–S polyhedra, forming a Co–Mo–S covalent network that hosts K + ions through electrostatic attraction. The interactions of KCMS with heavy metal ions, particularly Pb 2+ and Ag + , reveal that KCMS is exceptionally effective in removing these ions from ppm concentrations down to trace levels (≤5 ppb). KCMS rapidly removes Ag + (≈81.7%) and Pb 2+ (≈99.5%) within five minutes, achieving >99.9% removal within an hour, with a distribution constant K d ≥10 8 mL g -1 . KCMS exhibits an impressive removal capacity of 1378 mg g -1 for Ag + and 1146 mg g -1 for Pb 2+ , establishing it as one of the most effective materials known to date for heavy metal removal. This material is also effective for the removal of Ag + and Pb 2+ along with Hg 2+ , Ni 2+ , Cu 2+ , and Cd 2+ from various water sources even in the presence of highly concentrated and chemically diverse cations, anions, and organic species. Analysis of the post-interacted KCMS by synchrotron X-ray pair distribution function (PDF), X-ray photoelectron spectroscopy (XPS) and energy dispersive X-ray spectroscopy (EDS) revealed that the sorption of Pb 2+ , Ag + , and Hg 2+ mainly occurs by the exchange of K + and Co 2+ . Despite being amorphous, this material exhibits unprecedented ion-exchange mechanisms both for the ionically and covalently bound K + and Co 2+ , respectively. In conclusion, this discovery advances our knowledge of amorphous gels and guides material synthesis principles for the highly selective and efficient removal of heavy metal ions from water.

54 ENVIRONMENTAL SCIENCES↗

Atlas of putative minima and low-lying energy networks of water clusters n = 3 - 25

We report a database consisting of the putative minima and ~3.2 million local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3 – 25 using an improved version of the Monte Carlo Temperature Basin Paving (MCTBP) global optimization procedure in conjunction with the ab-initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as lowlying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation MP2/aug-cc-pVTZ level of theory. In total, we have identified 3,138,303 networks corresponding to local minima of the clusters n = 3 – 25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1 – 3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but no limited to neural networks, machine and deep learning) to understand the properties of water, nature’s most important substance.

Rakshit, Avijit↗

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗

Simulations of activities, solubilities, transport properties, and nucleation rates for aqueous electrolyte solutions

This article reviews recent molecular simulation studies of "collective" properties of aqueous electrolyte solutions, specifically free energies and activity coefficients, solubilities, nucleation rates of crystals, and transport coefficients. These are important fundamental properties for biology and geoscience, but also relevant for many technological applications. Their determination from molecular-scale calculations requires large systems and long sampling times, as well as specialized sampling algorithms. As a result, such properties have not typically been taken into account during optimization of force field parameters; thus, they provide stringent tests for the transferability and range of applicability of proposed molecular models. There has been significant progress on simulation algorithms to enable the determination of these properties with good statistical uncertainties. Comparisons of simulation results to experimental data reveal deficiencies shared by many commonly used models. Moreover, there appear to exist specific tradeoffs within existing modeling frameworks, so that good prediction of some properties is linked to poor prediction for specific other properties. For example, non-polarizable models that utilize full charges on the ions generally fail to predict accurately both activity coefficients and solubilities; the concentration dependence of viscosity and diffusivity for these models is also incorrect. Scaled-charge models improve the dynamic properties and could also perform well for solubilities, but fail in the prediction of nucleation rates. Even models that do well at room temperature for some properties generally fail to capture their experimentally observed temperature dependence. Finally, the main conclusion from the present review is that qualitatively new physics will need to be incorporated in future models of electrolyte solutions to allow description of collective properties for broad ranges of concentrations, temperatures, and solvent conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanoscale structure of a hybrid aqueous–nonaqueous electrolyte

A new class of electrolytes have been reported, hybridizing aqueous with non-aqueous solvents, which combines non-flammability and non-toxicity characteristics of aqueous electrolytes with the superior electrochemical stability of non-aqueous systems. Here, we report measurements of the structure of an electrolyte composed of an equal-mass mixture of 21 m LiTFSI-water and 9 m LiTFSI-dimethyl carbonate using high-energy x-ray diffraction and polarized neutron diffraction with isotope substitution. Neutron structure factors from partially and fully deuterated samples exhibit peaks at low scattering vector Q that we ascribe to long-range correlations involving both solvent molecules and TFSI − anions. Here we compare both sets of measurements with results of molecular dynamics simulations based on a polarizable force field. The structures derived from simulations are generally in agreement with those measured, except that neutron structure factors predicted for two partially deuterated samples show very intense scattering increasing up to the low- Q limit of simulation, indicating a partial segregation between the two solvents not observed in experimental measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Gaussian field approach to the planar electric double layer structures in electrolyte solutions

Here, in this work, the planar, electric, double-layer structures of non-polarizable electrodes in electrolyte solutions are studied with Gaussian field theory. A response function with two Yukawa functions is used to capture the electrostatic response of the electrolyte solution, from which the modified response function in the planar symmetry is derived analytically. The modified response function is further used to evaluate the induced charge density and the electrostatic potential near an electrode. The Gaussian field theory, combined with a two-Yukawa response function, can reproduce the oscillatory decay behavior of the electric potentials in concentrated electrolyte solutions. When the exact sum rules for the bulk electrolyte solutions and the electric double layers are used as constraints to determine the parameters of the response function, the Gaussian field theory could at least partly capture the nonlinear response effect of the surface charge density. Comparison with results for a planar electrode with fixed surface charge densities from molecular simulations demonstrates the validity of Gaussian field theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗