Search NASA⌕ Search

SEARCH · Search NASA

Results for “POLAR COORDINATE”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2021

Time Refraction of Spin Waves

Here, we present an experimental study of time refraction of spin waves (SWs) propagating in microscopic waveguides under the influence of time-varying magnetic fields. Using space- and time-resolved Brillouin light scattering microscopy, we demonstrate that the broken translational symmetry along the time coordinate results in a loss of energy conservation for SWs and thus allows for a broadband and controllable shift of the SW frequency. With an integrated design of SW waveguide and microscopic current line for the generation of strong, nanosecond-long, magnetic field pulses, a conversion efficiency up to 39% of the carrier SW frequency is achieved, significantly larger compared to photonic systems. Given the strength of the magnetic field pulses and its strong impact on the SW dispersion relation, the effect of time refraction can be quantified on a length scale comparable to the SW wavelength. Furthermore, we utilize time refraction to excite SW bursts with pulse durations in the nanosecond range and a frequency shift depending on the pulse polarity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Warm Arctic, Cold Siberia Pattern: Role of Full Arctic Amplification Versus Sea Ice Loss Alone

The effect of future Arctic amplification (AA) on the extratropical atmospheric circulation remains unclear in modeling studies. Using a collection of coordinated atmospheric and coupled global climate model perturbation experiments, we find an emergent relationship between the high-latitude 1,000–500 hPa thickness response and an enhancement of the Siberian High in winter. This wave number-1-like sea level pressure anomaly pattern is linked to an equatorward shift of the eddy-driven jet and a dynamical cooling response in eastern Asia. Additional simulations, where AA is imposed directly into the model domain by nudging, demonstrate how the sea ice forcing is insufficient by itself to capture the vertical extent of the warming and by extension the amplitude of the response in the Siberian High. Additionally, this study demonstrates the importance of the vertical extent of the tropospheric warming over the polar cap in revealing the “warm Arctic, cold Siberia” anomaly pattern in future projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

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

CO and H 2 adsorption on Au-Ni bimetallic surfaces: a combined experimental and DFT theoretical study

Au-Ni bimetallic thin films were grown on refractory metal substrates. CO and H 2 adsorption on Au-Ni bimetallic surfaces have been studied by a combination of in situ polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS), temperature-programmed desorption (TPD), and density functional theory (DFT) calculations. It is found CO desorption peak shifts from 413 K on pure Ni surfaces to 293 K on the isolated Ni atoms formed by alloying with Au atoms. The sharp decrease of CO desorption temperature on Au-Ni surfaces with increasing Au coverage is caused by the change of the favored CO adsorption sites from bridge/hollow sites on pure Ni surfaces to Ni top sites on Au-Ni bimetallic surfaces. In situ PM-IRRAS shows two CO adsorption bands on Au sites at 2119 cm −1 and 2103 cm −1 on Au-Ni surfaces at 80 K, which are due to CO bound on under-coordinated Au atoms and electron negatively charged Au sites modified with nearby Ni atoms, respectively. Even with the Au-Ni surface temperature at as low as 100 K, CO adsorption induced Ni surface segregation has been observed by in situ PM-IRRAS. Furthermore, DFT calculation results discover the adsorption energy of CO on Ni top sites continues to decrease with increasing Au coverage due to the geometric ensemble effect and the lowered d-band center after Ni alloying with Au. H 2 desorption temperature decreases from 363 K on pure Ni thin films to 302 K with increasing Au coverage to 0.6 ML. A new H 2 peak appears at around 170 K on the Au-Ni surfaces with Au coverages between 0.6 ML and 0.9 ML. This new H 2 TPD peak is assigned to H 2 desorption from the totally isolated Ni sites. With Au coverage above 1.5 ML, there is no any H 2 desorption detected. Finally, the combined surface science studies and DFT calculations provide new insights into the surface structure-activity correlation of Ni-base bimetallic surface alloys.

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

Atomistic insight into the ferroelastic post-stishovite transition by high-pressure single-crystal X-ray diffraction

Abstract The post-stishovite transition is a classic pseudo-proper typed ferroelastic transition with a symmetry-breaking spontaneous strain. This transition has been studied using high-pressure spontaneous strains, optic modes, and elastic moduli (Cij) based on the Landau modeling, but its atomistic information and structural distortion remain poorly understood. Here we have conducted synchrotron single-crystal X-ray diffraction measurements on stishovite crystals up to 75.3 GPa in a diamond-anvil cell. Analysis of the data reveals atomic positions, bond lengths, bond angles, and variations of SiO6 octahedra across the transition at high pressure. Our results show that the O coordinates split at ~51.4 GPa, where the apical and equatorial Si-O bond lengths cross over, the SiO6 octahedral distortion vanishes, and the SiO6 octahedra start to rotate about the c axis. Moreover, distortion mode analysis shows that an in-plane stretching distortion (GM1+ mode) occurs in the stishovite structure at high pressure while a rotational distortion (GM2+ mode) becomes dominant in the post-stishovite structure. These results are used to correlate with elastic moduli and Landau parameters (symmetry-breaking strain e1–e2 and order parameter Q) to provide atomistic insight into the ferroelastic transition. When the bond lengths of two Si-O bonds are equal due to the contribution from the GM1+ stretching mode, C11 converges with C12, and the shear wave VS1[110] polarizing along [110] and propagating along [110] vanishes. Values of e1–e2 and Q are proportional to the SiO6 rotation angle from the occurrence of the GM1+ rotational mode in the post-stishovite structure. Our results on the pseudo-proper type transition are also compared with that for the proper type in albite and improper type in CaSiO3 perovskite. The symmetry-breaking strain, in all these types of transitions, arises as the primary effect from the structural angle (such as SiO6 rotation or lattice constant angle) and its relevant distortion mode in the low-symmetry ferroelastic phase.

Geochemistry & Geophysics↗
DOE OSTI2022

Stabilized Synthesis of 2D Verbeekite: Monoclinic PdSe 2 Crystals with High Mobility and In-Plane Optical and Electrical Anisotropy

PdSe 2 has a layered structure with an unusual, puckered Cairo pentagonal tiling. Its atomic bond configuration features planar 4-fold-coordinated Pd atoms and intralayer Se–Se bonds that enable polymorphic phases with distinct electronic and quantum properties, especially when atomically thin. PdSe 2 is conventionally orthorhombic, and direct synthesis of its metastable polymorphic phases is still a challenge. Here, we report an ambient-pressure chemical vapor deposition approach to synthesize metastable monoclinic PdSe 2 . Monoclinic PdSe 2 is shown to be synthesized selectively under Se-deficient conditions that induce Se vacancies. These defects are shown by first-principles density functional theory calculations to reduce the free energy of the metastable monoclinic phase, thereby stabilizing it during synthesis. Further, the structure and composition of the monoclinic PdSe 2 crystals are identified and characterized by scanning transmission electron microscopy imaging, convergent beam electron diffraction, and electron energy loss spectroscopy. Polarized Raman spectroscopy of the monoclinic PdSe 2 flakes reveals their strong in-plane optical anisotropy. Electrical transport measurements show that the monoclinic PdSe 2 exhibits n-type charge carrier conduction with electron mobilities up to ~298 cm 2 V –1 s –1 and a strong in-plane electron mobility anisotropy of ~1.9. The defect-mediated growth pathway identified in this work is promising for phase-selective direct synthesis of other 2D transition metal dichalcogenides.

2D material↗
DOE OSTI2021

Elucidating the Location of Cd 2+ in Post-synthetically Treated InP Quantum Dots Using Dynamic Nuclear Polarization 31 P and 113 Cd Solid-State NMR Spectroscopy

Indium phosphide quantum dots (InP QD) are a promising alternative to traditional QD materials that contain toxic heavy elements such as lead and cadmium. However, InP QD obtained from colloidal synthesis are often plagued by poor photoluminescence quantum yields (PL-QYs). In order to improve the PL-QY of InP QD, a number of post-synthetic treatments have been devised. Recently, it has been shown that InP post-synthetically treated with Lewis acid metal divalent cations (M-InP) exhibit enhanced PL-QY; however, the molecular structure and mechanism behind the improved PL-QY are not fully understood. Here, to determine the surface structure of M-InP QD, dynamic nuclear polarization surface-enhanced nuclear magnetic resonance spectroscopy (DNP SENS) experiments were employed on a series of InP magic size clusters treated with Cd ions, InP QD, cadmium phosphide (Cd 3 P 2 ) QD, and Cd-treated InP QD (Cd–InP QD). With the use of DNP SENS, we were able to obtain the 1D 31 P and 113 Cd NMR spectra, 113 Cd{ 31 P} rotational-echo double-resonance (REDOR) NMR spectra, and 31 P{ 113 Cd} dipolar heteronuclear multiple quantum correlation (D-HMQC) sequence. Changes in the phosphide 31 P chemical shifts after Cd treatment provide indirect evidence that some Cd alloys into the sub-surface regions of the particle. DNP-enhanced 113 Cd solid-state NMR spectra suggest that most Cd ions are coordinated by oxygen atoms from either carboxylate ligands or surface phosphate groups. 113 Cd{ 31 P} REDOR and 31 P{ 113 Cd} D-HMQC experiments confirm that a subset of Cd ions are located on the surface of Cd–InP QD and coordinated with phosphate groups.

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

Data for "Plasmon-driven exciton formation in a non-equilibrium Fermi liquid"

This repository contains source data for key plots presented in the manuscript "Plasmon-driven exciton formation in a non-equilibrium Fermi liquid." Experimental data that was analyzed in Igor Pro 8 are presented as the .pxp files used to generate individual sub-plots. Electronic spectral function calculations are provided as .txt files, in which consecutive rows refer to the meshgrid x coordinate, y coordinate, spectral function (and, where relevant, axis-projected local angular momentum). We additionally include the Wannier model and DFT-obtained bulk band structure on which the Wannier model was based. Files are named as the number of the figure in the manuscript to which they correspond, with additional details included where necessary. Details of file names: 2a_DOS_Lxz_Ek_KGM_40layer_xnum_800kpt_tot.txt: Density of states, xz-axis projected local orbital angular momentum, for 800 points along the K-Gamma-M path, for a 40-layer model. 2c_composite_y.pxp: ARPES (angle-resolved photoemission spectroscopy) spectra along the ky axis, including both a scan near the Fermi level and a scan at high kinetic energies. 2d_LCP_RCP_diff_Sect_20K.pxp: difference between ARPES constant energy cuts at T=20 K at E0 + 0.23 eV taken with left- and right-circularly polarized photons. The polarization-integrated intensity at the constant energy cut is also included. 2e_DOS_L45_E11pt79_m0pt25to0pt25_xnum_800kpt_tot.txt: Density of states, xz-projected local orbital angular momentum, and corresponding k-points in two dimensions from ab-initio electronic structure calculations for a constant-energy cut. 3a_[x]_[y]ps: ARPES cut under excitation at a fluence of x uJ/cm2, measured y ps after photoexcitation. Measurements were performed at 9 K. 3b_[x]: Energy distribution curves under excitation at a fluence x uJ/cm2 at selected delay times after photoexcitation. 4a_ImSigma_vs_temperature.pxp: Imaginary self energy (extracted from ARPES linewidths) at different energies above E0 for selected lattice temperatures. 4b_EELS_lowE.pxp: Electron energy loss spectrum over a low energy range 5b_diff_55m15.pxp: Difference between momentum-integrated Tr-ARPES traces at 55 uJ/cm2 and 15 uJ/cm2 photoexcitation. Time-dependent intensity at each energy level has been normalized to a maximum of 1 for each individual fluence prior to subtraction. 5d_invtau_at_EX_vs_fluence.pxp: decay rate at a specified energy EX for different excitation fluences, from single exponential fits. NOTE: Analyses based on the Wannier model presented here should cite both the associated Article and this dataset. For all other files in the repository, citing the dataset alone is sufficient.

Acharya, Rishi [University of Illinois] (ORCID:000↗
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

Volumetric imaging of the 3D orientation of cellular structures with a polarized fluorescence light-sheet microscope

Polarized fluorescence microscopy is a valuable tool for measuring molecular orientations in biological samples, but techniques for recovering three-dimensional orientations and positions of fluorescent ensembles are limited. We report a polarized dual-view light-sheet system for determining the diffraction-limited three-dimensional distribution of the orientations and positions of ensembles of fluorescent dipoles that label biological structures. We share a set of visualization, histogram, and profiling tools for interpreting these positions and orientations. We model the distributions based on the polarization-dependent efficiency of excitation and detection of emitted fluorescence, using coarse-grained representations we call orientation distribution functions (ODFs). We apply ODFs to create physics-informed models of image formation with spatio-angular point-spread and transfer functions. We use theory and experiment to conclude that light-sheet tilting is a necessary part of our design for recovering all three-dimensional orientations. We use our system to extend known two-dimensional results to three dimensions in FM1-43-labeled giant unilamellar vesicles, fast-scarlet-labeled cellulose in xylem cells, and phalloidin-labeled actin in U2OS cells. Additionally, we observe phalloidin-labeled actin in mouse fibroblasts grown on grids of labeled nanowires and identify correlations between local actin alignment and global cell-scale orientation, indicating cellular coordination across length scales.

Science & Technology - Other Topics↗
DOE OSTI2022

A Heterogeneous Palladium Catalyst for the Polymerization of Olefins Prepared by Halide Abstraction Using Surface R 3 Si + Species

The silylium-like surface species [ i Pr 3 Si][(R F O) 3 Al-OSi≡)] activates (N^N)Pd(CH 3 )Cl (N^N=Ar-N=CMeMeC=N-Ar, Ar=2,6-bis(diphenylmethyl)-4-methylbenzene) by chloride ion abstraction to form [(N^N)Pd-CH 3 ][(R F O) 3 Al-OSi≡)] ( 1 ). A combination of FTIR, solid-state NMR spectroscopy, and reactions with CO or vinyl chloride establish that 1 shows similar reactivity patterns as (N^N)Pd(CH 3 )Cl activated with Na[B(Ar F ) 4 ]. Multinuclear 13 C{ 27 Al} RESPDOR and 1 H{ 19 F} S-REDOR experiments are consistent with a weakly coordinated ion-pair between (N^N)Pd-CH 3 + and [(R F O) 3 Al-OSi≡)]. 1 catalyzes the polymerization of ethylene with similar activities as [(N^N)Pd-CH 3 ] + in solution and incorporates up to 0.4% methyl acrylate in copolymerization reactions. 1 produces polymers with significantly higher molecular weight than the solution catalyst, and generates the highest molecular weight polymers currently reported in copolymerization reactions of ethylene and methylacrylate.

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

Modeling and Rapid Prototyping of Integrated Transmission-Distribution OPF Formulations with PowerModelsITD.jl

Conventional electric power systems are composed of different unidirectional power flow stages of generation, transmission, and distribution, managed independently by transmission system and distribution system operators. However, as distribution systems increase in complexity due to the integration of distributed energy resources, coordination between transmission and distribution networks will be imperative for the optimal operation of the power grid. However, coupling models and formulations between transmission and distribution is non-trivial, in particular due to the common practice of modeling transmission systems as single-phase, and distribution systems as multi-conductor phase-unbalanced. To enable the rapid prototyping of power flow formulations, in particular in the modeling of the boundary conditions between these two seemingly incompatible data models, we introduce PowerModelsITD.jl, a free, open-source toolkit written in Julia for integrated transmission-distribution (ITD) optimization that leverages mature optimization libraries from the InfrastructureModels.jl-ecosystem. The primary objective of the proposed framework is to provide baseline implementations of steady-state ITD optimization problems, while providing a common platform for the evaluation of emerging formulations and optimization problems. In this work, we introduce the nonlinear formulations currently supported in PowerModelsITD.jl, which include AC-polar, AC-rectangular, current-voltage, and a linear network transportation model. Results are validated using combinations of IEEE transmission and distribution networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Lower-length-scale modeling of chemical additions, corrosion and fission products in select MSR base salts

Molten salts play a crucial role in advancing Generation IV nuclear energy technology, with chloride-based salts like NaCl-UCl 3 garnering significant attention due to their distinctive properties. However, the corrosive nature of molten salts causes the dissolution of chromium (Cr), leading to the formation of CrCl 2 and CrCl 3 species in molten chlorides. Moreover, the radioactive decay of nuclear fuel gives rise to fission products, including Cs, Sr, and I chlorides. The first part of this report presents a comprehensive study utilizing ab initio molecular dynamics (AIMD) simulations to investigate the properties of eutectic NaCl-UCl 3 molten salt in the presence of corrosion products (CrCl 2 and CrCl 3 ) and fission product (CsCl). The study focuses on essential structural and thermophysical properties such as density, mixing energy, coordination numbers (CN), and Radial Distribution Functions (RDF) of the salts with varying compositions of corrosion products and fission product, covering a range from 0% to 13.5%. The results offer valuable insights into the behavior of corrosion and fission products in uranium-based molten salts, providing essential data that can be used as input to the MTDB-TC and MSTDB-TP property databases being developed by the NEAMS program. Due to their favorable characteristics such as low melting points, high boiling points, and low costs, MgCl 2 +NaCl+KCl (MNK) eutectic salts have recently attracted attention as high-temperature heat transfer fluids. The incorporation of LiCl into MNK salts can further reduce their melting points and increase their specific heat capacities, which is de- sirable for high-temperature heat transfer applications. The second part of this report presents the development and validation of a new shell-model potential for the MgCl 2 +NaCl+KCl+LiCl system, which captures the polarization of Cl anions. The extensive comparison with experimental data and AIMD simulations demonstrates the accuracy and reliability of the potential. Furthermore, using this potential, we elucidate the intricate network structure in MNK eutectic salts that is formed through polyhedron sharing. This research contributes to a fundamental understanding of the atomic structures and thermophysical properties of multi-component molten salts, which is crucial for future development of heat transfer fluids for applications in molten salt reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Mesoscale Design of Magnetoelectric Heterostructures and Nanocomposites

Multiferroic materials have a coexistence of at least two ferroic orders. The coupling interaction between the different order parameters can produce new phenomena, such as the magnetoelectric (ME) effect. The ME response is the appearance of an electric polarization P upon applying a magnetic field H, and/or vice versa. Magnetoelectricity has been observed as an intrinsic effect in some natural materials at low temperature; however, such single phase materials suffer from extremely weak ME exchange. Better alternatives are two-phase composites consisting of magnetostrictive and piezoelectric phases, which have been shown to have large ME coefficients. Said composites exploit a unique product ME tensor property that depends upon the individual strictions of the piezoelectric and magnetostrictive phases, the elastic stiffness and integrity of the interphase interfaces, the phase distribution and the dimensionality of its connectivity, and any elastic constraint; amongst other things. This research program between D. Viehland PI and S.Priya at Virginia Tech investigated the mesoscale physics and chemistry of such two-phase ME heterostructures. The results of the proposed research are relevant to new generations of multifunctional materials involving heterogeneous components. How the product ME property be optimized? The answers lies with the interphase-interfaces, where the atomic layers coordinate stress transfer from one phase to another. It was the primary goal of this research program to discover such fundamental relationship for two phase ME heterostructures and nanocomposites. This research program integrated all the components required to fully investigate the ME phenomenon: starting from synthesis of nanomaterials with diversified interfacial characteristics; to full electrical, magnetic, and structural characterizations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Process-based climate change assessment for European winds using EURO-CORDEX and global models

Wind energy is an important pillar of decarbonization strategies and potentially vulnerable to climate change. Existing wind climate change assessments rely on climate models but a systematic investigation of the global-to-regional climate modeling chain is missing. In this study, I highlight key limitations, namely (a) the differing representation of land use change in global and regional climate models which compromises comparability, and (b) the consistency of large-scale features along the global-to-regional climate modeling chain. To this end, I analyze the large European Coordinated Downscaling Experiment (EURO-CORDEX) ensemble (rcp85: N = 49; rcp45: N= 18; rcp26: N = 22) along with the driving global models (rcp85: N = 7; rcp45: N = 5; rcp26: N = 7), finding evidence that climate change reduces mean wind speeds by up to -0.8 ms -1 (offshore) and -0.3 ms -1 (onshore). I provide physical explanations for these changes by identifying two key drivers. First, onshore wind speeds drop in the driving global models in regions and scenarios with strong land use change but show no drop in EURO-CORDEX where land use is held constant. Second, offshore wind reductions follow decreases in the equator-to-pole temperature gradient remarkably well with correlations reaching around 0.9 in resource-rich European countries like Ireland, the United Kingdom and Norway, implying that arctic amplification is a severe risk for European offshore wind energy. My results suggest that earlier conclusions of negligible climate change impacts on wind energy might be premature if either land use changes strongly or polar amplification is at or above the range sampled in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The multiwavelength view of shocks in the fastest nova V1674 Her

Classical novae are shock-powered multiwavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t 2 = 1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV γ-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift, and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1–100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 ± 0.4) × 10 −6 photons cm −2 s −1 . Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to k $T$ shock = 4 keV. The lack of a detectable 6.7 keV Fe Kα emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin–Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova.

79 ASTRONOMY AND ASTROPHYSICS↗