Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical calculations, Polarization”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Full-Scale Ab Initio Simulation of Magic-Angle-Spinning Dynamic Nuclear Polarization

Theoretical models aimed at describing magic-angle-spinning (MAS) dynamic nuclear polarization (DNP) NMR have great potential in facilitating the in silico design of DNP polarizing agents and formulations. These models must typically face a trade-off between the accuracy of a strict quantum mechanical description and the need for using realistically large spin systems, for instance, using phenomenological models. In this paper, we show that the use of aggressive state-space restrictions and an optimization strategy allows full-scale ab initio MAS-DNP simulations of spin systems containing thousands of nuclei. Our simulations are shown to reproduce experimental DNP enhancements quantitatively, including their MAS rate dependence, for both frozen solutions and solid materials. They also reveal the importance of a previously unrecognized structural feature found in some polarizing agents that helps minimize the sensitivity losses imposed by the spin diffusion barrier.

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

High-throughput computation of electric polarization in solids via Berry flux diagonalization

Electric polarization in the absence of an externally applied electric field is a key property of polar materials, but the standard interpolation-based ab initio approach to compute polarization differences within the modern theory of polarization presents challenges for automated high-throughput calculations. Berry flux diagonalization [J. Bonini et al., Phys. Rev. B 102, 045141 (2020)] has been proposed as an efficient and reliable alternative, though it has yet to be widely deployed. Here, we assess Berry flux diagonalization using ab initio calculations of a large set of materials, introducing and validating heuristics that ensure branch alignment with a minimal number of intermediate interpolated structures. Our automated implementation of Berry flux diagonalization succeeds in cases where prior interpolation-based workflows fail due to band-gap closures or branch ambiguities. Benchmarking with ab initio calculations of 176 candidate ferroelectrics, we demonstrate the efficacy of the approach on a broad range of insulating materials and obtain accurate effective polarization values with fewer interpolated structures than prior automated interpolation-based workflows. Our real-space heuristics that can predict gauge stability a priori from ionic displacements enable a general automated framework for reliable polarization calculations and efficient high-throughput screening of chemically and structurally diverse polar insulators. These results establish Berry flux diagonalization as a robust and efficient method to compute the effective polarization of solids and to accelerate the data-driven discovery of functional polar materials.

Poteshman, Abigail N. [University of Chicago, IL (↗
DOE OSTI2023

Signatures of Non-universal Quantum Dynamics of Ultracold Chemical Reactions of Polar Alkali Dimer Molecules with Alkali Metal Atoms: Li( 2 S) + NaLi( a 3 Σ + ) → Na( 2 S) + Li 2 ( a 3 Σ u + )

Ultracold chemical reactions of weakly bound triplet-state alkali metal dimer molecules have recently attracted much experimental interest. Here we perform rigorous quantum scattering calculations with a new ab initio potential energy surface to explore the chemical reaction of spin-polarized NaLi(a 3 Σ + ) and Li( 2 S) to form Li 2 (a 3 Σ u + ) and Na( 2 S). The reaction is exothermic and proceeds readily at ultralow temperatures. Significantly, we observe strong sensitivity of the total reaction rate to small variations of the three-body part of the Li 2 Na interaction at short range, which we attribute to a relatively small number of open Li 2 (a 3 Σ u + ) product channels populated in the reaction. This provides the first signature of highly non-universal dynamics seen in rigorous quantum reactive scattering calculations of an ultracold exothermic insertion reaction involving a polar alkali dimer molecule, opening up the possibility of probing microscopic interactions in atom+molecule collision complexes via ultracold reactive scattering experiments.

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

Solvent Induced Proton Polarization within the Nuclear–Electronic Orbital Framework

To explicitly account for nuclear quantum effects and solvent environments in simulations of chemical processes, the nuclear-electronic orbital approach is coupled with a polarizable continuum model (PCM). This NEO–PCM approach is used to explore the influence of solvation on nuclear polarization through applications to a water dimer and a set of protonated water tetramers. Nuclear polarization in these species is analyzed in terms of changes in proton density and oxygen–hydrogen bond length. Solvation is shown to enhance nuclear polarization with increasing dielectric constant. For the water dimer, the internal, hydrogen-bonded proton is shown to polarize more than the external, free proton. Moreover, proton quantization leads to greater solvent polarization through their mutual polarization. Furthermore, these calculations highlight the complex interplay among electronic, nuclear, and solvent polarization in chemical systems.

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

Computational and experimental study of different brines in temperature swing solvent extraction desalination with amine solvents

Rapid global urbanization and high-salinity wastewater disposal from industrial activities have exerted significant pressure on water resources. Over the past few years, temperature swing solvent extraction (TSSE) has been identified as a promising technique to desalinate hypersaline brines. Despite its potential, the TSSE desalination literature has been mainly based on empirical insights, and the limited molecular simulation studies have primarily focused on NaCl brines. Herein, we use molecular dynamics (MD) simulations to study the TSSE desalination of four different brines, namely, KCl, KBr, NaCl, and NaBr using diisopropylamine as the solvent. Based on both bulk and interfacial brine-diisopropylamine MD results, here we investigate the qualitative and quantitative performance of the simulations by benchmarking these results against our experimental evaluations of these same systems. MD results provide satisfactory qualitative agreement with the experimental data of water solubilization in the organic phase and amine solubilization in the aqueous phase for the KBr, KCl, and NaBr brines. Also, the molecular mechanism of solvation of ionic species by water molecules over diisopropylamine suggested by the MD simulations is in agreement with our experimental data. However, larger qualitative and quantitative deviations were observed for the NaCl brines, and this is likely due to polarization and charge transfer effects, as quantified by our quantum chemical calculations.

42 ENGINEERING↗
DOE OSTI2022

Numerical Simulation of Vibrational Sum Frequency Generation Intensity for Non-Centrosymmetric Domains Interspersed in an Amorphous Matrix: A Case Study for Cellulose in Plant Cell Wall

Vibrational sum frequency generation (SFG) spectroscopy can specifically probe molecular species non-centrosymmetrically arranged in a centrosymmetric or isotropic medium. This capability has been extensively utilized to detect and study molecular species present at the two-dimensional (2D) interface at which the centrosymmetry or isotropy of bulk phases is naturally broken. The same principle has been demonstrated to be very effective for the selective detection of non-centrosymmetric crystalline nanodomains interspersed in three-dimensional (3D) amorphous phases. However, the full spectral interpretation of SFG features has been difficult due to the complexity associated with the theoretical calculation of SFG responses of such 3D systems. This paper describes a numerical method to predict the relative SFG intensities of non-centrosymmetric nanodomains in 3D systems as functions of their size and concentration as well as their assembly patterns, i.e., the distributions of tilt, azimuth, and rotation angles with respect to the lab coordinate. We applied the developed method to predict changes in the CH and OH stretch modes characteristic to crystalline cellulose microfibrils distributed with various orders, which are relevant to plant cell wall structures. As a result, the same algorithm can also be applied to any SFG-active nanodomains interspersed in 3D amorphous matrices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Proximity ferroelectricity in wurtzite heterostructures

Proximity ferroelectricity is an interface-associated phenomenon in electric-field-driven polarization reversal in a non-ferroelectric polar material induced by one or more adjacent ferroelectric materials. Here, in this paper, we report proximity ferroelectricity in wurtzite ferroelectric heterostructures. In the present case, the non-ferroelectric layers are AlN and ZnO, whereas the ferroelectric layers are Al 1−x B x N, Al 1−x Sc x N and Zn 1−x Mg x O. The layered structures include nitride–nitride, oxide–oxide and nitride–oxide stacks that feature two-layer (asymmetric) and three-layer (symmetric) configurations. Ferroelectric switching in both layers is validated by multimodal characterization methods, including polarization hysteresis, anisotropic chemical etching, second harmonic generation, piezo response force microscopy, electromechanical testing and atomic resolution polarization orientation imaging in real space by scanning transmission electron microscopy. We present a physical switching model in which antipolar nuclei originate in the ferroelectric layer and propagate towards the internal non-ferroelectric interface. The domain wall leading edge produces elastic and electric fields that extend beyond the interface at close proximity, reducing the switching barrier in the non-ferroelectric layer, and allowing complete domain propagation without breakdown. Density functional theory calculations of polymorph energies, reversal barriers and domain wall energies support this model. Proximity ferroelectricity enables polarization reversal in wurtzites without the chemical or structural disorder that accompanies elemental substitution, opening new questions and opportunities regarding interface-based ferroelectricity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics of Aqueous Electrolyte Solutions: Challenges for Simulations

This Perspective article focuses on recent simulation work on the dynamics of aqueous electrolytes. It is well-established that full-charge, nonpolarizable models for water and ions generally predict solution dynamics that are too slow in comparison to experiments. Models with reduced (scaled) charges do better for solution diffusivities and viscosities but encounter issues describing other dynamic phenomena such as nucleation rates of crystals from solution. Polarizable models show promise, especially when appropriately parametrized, but may still miss important physical effects such as charge transfer. First-principles calculations are starting to emerge for these properties that are in principle able to capture polarization, charge transfer, and chemical transformations in solution. Finally, while direct ab initio simulations are still too slow for simulations of large systems over long time scales, machine-learning models trained on appropriate first-principles data show significant promise for accurate and transferable modeling of electrolyte solution dynamics.

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

Electron–Nucleus Hyperfine Coupling Calculated from Restricted Active Space Wavefunctions and an Exact Two-Component Hamiltonian

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab-initio wavefunction calculations at the multi-reference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron – nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on non-relativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had only limited applicability. Furthermore, results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex, demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

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

Lone Pair Rotational Dynamics in Solids

Traditional classifications of crystalline phases focus on nuclear degrees of freedom. Through examination of both electronic and nuclear structure, we introduce the concept of an electronic plastic crystal. Such a material is classified by crystalline nuclear structure, while localized electronic degrees of freedom — here lone pairs — exhibit orientational motion at finite temperatures. This orientational motion is an emergent phenomenon arising from the coupling between electronic structure and polarization fluctuations generated by collective motions, such as phonons. Using ab initio molecular dynamics simulations, we predict the existence of electronic plastic crystal motion in halogen crystals and halide perovskites, and suggest that such motion may be found in a broad range of solids with lone pair electrons. Finally, such fluctuations in the charge density should be observable, in principle via synchrotron scattering.

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

Surface reservoirs dominate dynamic gas-surface partitioning of many indoor air constituents

Human health is affected by indoor air quality. One distinctive aspect of the indoor environment is its very large surface area that acts as a poorly characterized sink and source of gas-phase chemicals. In this work, air-surface interactions of 19 common indoor air contaminants with diverse properties and sources were monitored in a house using fast-response, on-line mass spectrometric and spectroscopic methods. Enhanced-ventilation experiments demonstrate that most of the contaminants reside in the surface reservoirs and not, as expected, in the gas phase. They participate in rapid air-surface partitioning that is much faster than air exchange. Phase distribution calculations are consistent with the observations when assuming simultaneous equilibria between air and large weakly polar and polar absorptive surface reservoirs, with acid-base dissociation in the polar reservoir. Chemical exposure assessments must account for the finding that contaminants that are fully volatile under outdoor air conditions instead behave as semivolatile compounds indoors.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

Data-driven computational prediction and experimental realization of exotic perovskite-related polar magnets

Rational design of technologically important exotic perovskites is hampered by the insufficient geometrical descriptors and costly and extremely high-pressure synthesis, while the big-data driven compositional identification and precise prediction entangles full understanding of the possible polymorphs and complicated multidimensional calculations of the chemical and thermodynamic parameter space. Here we present a rapid systematic data-mining-driven approach to design exotic perovskites in a high-throughput and discovery speed of the A 2 BB ’O 6 family as exemplified in A 3 TeO 6 . The magnetoelectric polar magnet Co 3 TeO 6 , which is theoretically recognized and experimentally realized at 5 GPa from the six possible polymorphs, undergoes two magnetic transitions at 24 and 58 K and exhibits helical spin structure accompanied by magnetoelastic and magnetoelectric coupling. We expect the applied approach will accelerate the systematic and rapid discovery of new exotic perovskites in a high-throughput manner and can be extended to arbitrary applications in other families.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic Scale Crystal Field Mapping of Polar Vortices in Oxide Superlattices

Polar vortices in oxide superlattices can be utilized as potential candidates for data storage applications due to their unique polarization topologies. The structure and dipole arrangement in polar vortices has been studied via X-ray scattering techniques, (scanning) transmission electron microscopy ((S)/TEM) and computational calculations. However, the fundamental correlation between the atomic structure and the electronic structure (which is manifested in the chemical bonding) has heretofore not been explored. Further, The hybridization between nominally empty $\textit{d}$ orbitals on the B-site with the occupied O 2$\textit{p}$ orbitals favors the condensation of a polar (ferroelectric) state in $\textit{AB}$O 3 perovskite oxides. The complex, continuously rotating local polarization texture of the vortices, in turn, can result in especially intricate $\textit{d}$- orbital interactions. Soft X-ray spectroscopy can probe these interactions at the transition metal $\textit{L}$-edge, but these techniques do not have the spatial resolution to resolve variations within one vortex (~5 nm region). Electron energy loss spectroscopy (EELS) in the STEM mode uses inelastically scattered electrons to probe the core-shell excitations (empty density of states) of transition metals at atomic resolution.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Structural, electronic, and polarization properties of YN and LaN

ScN has attracted great attention for its electronic properties and its ability to enhance polarization of AlN; however, its sister compounds, YN and LaN, remain much less studied. Here, we use first-principles calculations to evaluate YN and LaN in their cubic and hexagonal phases. Rocksalt YN and LaN are semiconductors, although we show that LaN differs from ScN and YN in having a direct band gap, which we attribute to its weaker p-p coupling. Both have low electron effective masses. In addition to their rocksalt structures, we evaluate the layered hexagonal and wurtzite phases of YN and LaN. For YN, the wurtzite phase cannot be stabilized, and hexagonal YN is higher in energy than rocksalt YN. In contrast, for LaN, the wurtzite phase is favored, and it is comparable in energy to rocksalt LaN. Wurtzite LaN has a polarization of 0.608 C/m 2 (referenced to the centrosymmetric layered hexagonal phase), and a high piezoelectric coefficient e 33 = 1.78 C/m 2 . Interestingly, we find that the polarization of wurtzite LaN may be reversible; we find a relatively small switching barrier of 0.06 eV per formula unit, offering the potential for its use as a ferroelectric. Since wurtzite LaN is closely lattice matched to InP, we investigate a heterostructure between (0001) wurtzite LaN and (111) zinc-blende InP, and find the polarization discontinuity would yield a bound charge of 1.3 × 10 14 e/cm 2 , offering the potential for novel electronic applications such as tunnel junctions. Furthermore, our results compare and contrast ScN, YN, and LaN, and highlight the potential of these materials for adoption in electronic and ferroelectric devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Giant Magnetic and Optical Anisotropy in Cerium-Substituted M-Type Strontium Hexaferrite Driven by $4f$ Electrons

By performing density-functional calculations, we find a giant magnetocrystalline anisotropy (MCA) constant in abundant element cerium $\mathrm{Ce}$ substituted M-type hexaferrite, in the energetically favorable strontium site, assisted by a quantum confined electron transfer from $\mathrm{Ce}$ to a specific iron $\mathrm{2a}$ site. Remarkably, the calculated electronic structure shows that the electron transfer leads to the formation of $\mathrm{Ce}^{3+}$ and $\mathrm{Fe}^{2+}$ at the $\mathrm{2a}$ site producing an occupied $\mathrm{Ce}$$(4f^1)$ state below the Fermi level that adds a significant contribution to MCA and magnetic moment. A half $\mathrm{Ce}$ substitution forms a metallic state, while a full substitution retains the semiconducting state of the strontium hexaferrite (host). In the latter, the band gap is reduced due to the formation of charge-transferred states in the gap region of the host. The optical absorption coefficient shows an enhanced anisotropy between light polarization in parallel and perpendicular directions. Calculated formation energies, including the analysis of probable competing phases, and elastic constants confirm that both compositions are chemically and mechanically stable. With successful synthesis, the $\mathrm{Ce}$ hexaferrite can be an alternative high-performing critical-element-free permanent magnet material adapted for use in devices such as automotive traction drive motors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

High-Throughput Screening and Accurate Prediction of Ionic Liquid Viscosities Using Interpretable Machine Learning

Ionic liquids (ILs) are a novel group of green solvents with great promise for various industrial applications, including carbon capture and lignocellulosic biomass deconstruction. However, the use of ILs at the industrial scale remains challenging due to their high viscosities at ambient temperatures. To develop ILs with lower viscosities, a systematic study of their quantitative structure–property relationship (QSPR) is desirable. Here, we developed four machine learning (ML) models to predict viscosity at various temperature and pressure ranges, trained over a wide range of ILs consisting of various cationic and anionic families. ML methods including two-factor polynomial regression (two-factor PR), support vector regression (SVR), feed-forward neural networks (FFNN), and categorical boosting (CATBoost) were developed based on features that have proven useful in previous ML studies: COSMO-RS (conductor-like screening model for real solvents)-derived surface screening charge densities (sigma profiles). FFNN and CATBoost were the most accurate in predicting IL viscosities with lower average absolute relative deviation and higher R2 values on the test set. Tanimoto similarity scores were calculated to characterize the chemical space and structural similarity of the investigated ions. Furthermore, SHapley Additive exPlanation (SHAP) analysis was employed to interpret the ML results. Temperature, the polar area of ILs, and the nonpolar regions of ions are key features that influence the viscosity predictions. Importantly, the IL viscosity prediction here is the most accurate reported to date.

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

van der Waals Epitaxy, Superlubricity, and Polarization of the 2D Ferroelectric SnS

Tin monosulfide (SnS) is a two-dimensional layered semiconductor that exhibits in-plane ferroelectric order at very small thicknesses and is of interest in highly scaled devices. Here, in this study, we report the epitaxial growth of SnS on hexagonal boron nitride (hBN) using a pulsed metal-organic chemical vapor deposition process. Lattice matching is observed between the SnS($100$) and hBN{$1\bar{10}$} planes, with no evidence of strain. Atomic force microscopy reveals superlubricity along the commensurate direction of the SnS/hBN interface, and first-principles calculations suggest that friction is controlled by the edges of the SnS islands, rather than interface interactions. Differential phase contrast imaging detects remnant polarization in SnS islands with domains that are not dictated by step-edges in the SnS. The growth of ferroelectric SnS on high quality hBN substrates is a promising step toward electrically switchable ferroelectric semiconducting devices.

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

Pareto Optimization of Oligomer Polarizability and Dipole Moment Using a Genetic Algorithm

High-performance electronic components are highly sought after in order to produce increasingly smaller and cheaper electronic devices. Drawing inspiration from inorganic dielectric materials, in which both polarizability and polarization contribute, organic materials can also maximize both. For a large set of small molecules drawn from PubChem, a Pareto-like front appears between the polarizability and dipole moment, indicating the presence of an apparent trade-off between these two properties. We tested this balance in π-conjugated materials by searching for novel conjugated hexamers with simultaneously large polar- izabilities and dipole moments with potential use for dielectric materials. Using a genetic algorithm (GA) screening technique in conjunction with an approximate density functional tight-binding method for property calculations, we were able to efficiently search chemical space for optimal hexamers. Given the scope of chemical space, using the GA technique saves considerable time and resources by speeding up molecular searches compared to a systematic search. Here, we also explored the underlying structure–function relationships, including sequence and monomer properties, that characterize large polarizability and dipole moment regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗