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

Whence the Interstellar Magnetic Field Shaping the Heliosphere?

Measurements of starlight polarized by aligned interstellar dust grains are used to probe the relation between the orientation of the ambient interstellar magnetic field (ISMF) and the ISMF traced by the ribbons of energetic neutral atoms discovered by the Interstellar Boundary Explorer spacecraft. We utilize polarization data, many acquired specifically for this study, to trace the configuration of the ISMF within 40 pc. A statistical analysis yields a best-fit ISMF orientation, B magpol , aligned with Galactic coordinates ℓ = 42°, b = 49°. Further analysis shows the ISMF is more orderly for "downfield" stars located over 90° from Bmagpol. The data subset of downfield stars yields an orientation for the nearby ISMF at ecliptic coordinates λ, β ≈ 219° ± 15°, 43° ± 9° (Galactic coordinates l, b ≈ 40°, 56°, ±17°). This best-fit ISMF orientation from polarization data is close to the field direction obtained from ribbon models. This agreement suggests that the ISMF shaping the heliosphere belongs to an extended ordered magnetic field. Extended filamentary structures are found throughout the sky. A previously discovered filament traversing the heliosphere nose region, "Filament A," extends over 300° of the sky, and crosses the upwind direction of interstellar dust flowing into the heliosphere. Filament A overlaps the locations of the Voyager kilohertz emissions, three quasar intraday variables, cosmic microwave background (CMB) components, and the inflow direction of interstellar grains sampled by Ulysses and Galileo. These features are likely located in the upstream outer heliosheath where ISMF drapes over the heliosphere, suggesting Filament A coincides with a dusty magnetized plasma. A filament 55° long is aligned with a possible shock interface between local interstellar clouds. A dark spot in the CMB is seen within 5° of the filament and within 10° of the downfield ISMF direction. Two large magnetic arcs are centered on the directions of the heliotail. The overlap between CMB components and the aligned dust grains forming Filament A indicates the configuration of dust entrained in the ISMF interacting with the heliosphere provides a measurable foreground to the CMB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗
DOE OSTI2021

Coordinated Modeling of Electric Grid and Natural Gas Network Operations

Presentation based upon full report of a Colorado case study and coordination framework, which is available at https://www.nrel.gov/docs/fy20osti/77096.pdf. JISEA analysts Brian Sergi, Omar Guerra, and Bri-Mathias Hodge will present on recent JISEA work on coordination between the natural gas and electricity sectors. Power and gas are becoming increasingly interdependent but weren't designed to function together. How does greater coordination impact system operations with different levels of solar and wind penetrations? This free JISEA presentation will take place in a Webinar on Tuesday, May 11, 2021 at 12 p.m.: Learn how greater coordination between the natural gas and electricity sectors impacts system operations with different levels of solar and wind penetrations. Presentation information includes: Historical and projected data for natural gas consumption and power generation in the United States as demonstrated by 2014 East Coast Polar Vortex, 2021 Texas Winter Storm Uri. Showing coupling points, FERC identified need for better coordination, different levels of coordination (decision making and optimization control), coordination framework. Case study on Colorado Front Range, ramping requirements and gas nominations, real-time dispatch (June and December scenarios), results of analysis of impacts of coordination on unserved load. Total real-time gas offtakes by node, impacts on unserved gas, effect on CO2 emissions. Discussions on clusions drawn from Colorado Case study expanded co-simulation via HELICS, objectives of the HELICS+ natural gas use case, modeling of hydrogen blending, techno economic assessment of blending, and blending impacts on energy content and pressure.

analysis↗
DOE OSTI2023

Expeditious Coordination-Driven Construction of Hierarchically Nanoporous Barium Salts

Organic moieties-derived salts with permanent porosity and polarized channels have shown unique features and attractive performance in the field of adsorption, separation, and conduction. However, state-of-the-art organic salts generally rely on ionic interaction and hydrogen bonding formation to maintain the porous channels. The synthesis of organic moiety-derived salts with permanent accessible pores even after removal of the trapped guest molecules, and without the constraint of hydrogen bonding formation still remains a great challenge. Herein, we present an expeditious construction pathway to generate hierarchically nanoporous barium salts without hydrogen bonding formation. The strong ionic interaction of the barium cation and sulfonate anions led to rapid reaction equilibrium (~2 min), affording diverse barium-derived ionic polymer (Ba-IP) with permanent porosity and highly polarized channels. The produced Ba-IP materials with abundant cations and anions displayed high CO 2 /N 2 and CO 2 /CH 4 separation performance, with the selectivities reaching up to 89.5 and 280, respectively, at 273 K, surpassing most of the organic polymers functionalized by ionic moieties.

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

Catalytic consequences of hydrogen addition events and solvent-adsorbate interactions during guaiacol-H 2 reactions at the H 2 O-Ru(0 0 0 1) interface

Catalytic reactions of biomass-derived phenolics and H 2 occur on transition metal surfaces via competitive C–O cleavage and ring saturation pathways, with both requiring multiple hydrogen addition events before forming their respective rate limiting transition states. These events are markedly affected by solvent chemical identity, with polar protic solvents ionizing hydrogen adatoms (H*) to interfacial protons (H + ) and opening up new catalytic routes. Here, we establish the reaction coordinate space for guaiacol-H 2 reactions on Ru(0 0 0 1) using density functional theory and describe the atomic-scale effect of a polar protic solvent, H 2 O. Coupled H + and H* attack leads to quasi-equilibrated enol and keto intermediates as the precursors for C–O cleavage and ring saturation, respectively. For C–O cleavage, H 2 O solvent enables a lower energy pathway via concomitant transfer of the hydroxyl H + to the methoxy oxygen during C–OCH 3 cleavage, forming a charge separated [Ru(s)–(C 6 H 5 O – )…(H + )…OCH 3 ] transition state and reducing the barrier by up to 0.8 eV as compared to unassisted C–OCH 3 cleavage. For ring saturation, H* attack onto an unsaturated meta carbon is rate limiting with no direct solvent participation, suggesting that protic polar solvents selectively promote the C–O cleavage pathway. Taken together, we show that activating guaiacol for either C–O bond cleavage or ring saturation product formation depends on the reactive hydrogen identity (H* or H + ), enol/keto isomerization equilibrium, and accessibility of the proton assisted Car–OCH 3 cleavage transition state. Here, all such factors are tunable via changes to the solvent or metal identity.

09 BIOMASS FUELS↗
DOE OSTI2020

Tunable Layer Circular Photogalvanic Effect in Twisted Bilayers

In this work, we develop a general theory of the layer circular photogalvanic effect (LCPGE) in quasi-two-dimensional chiral bilayers, which refers to the appearance of a polarization-dependent, out-of-plane static dipole moment induced by circularly polarized light. We elucidate the geometric origin of the LCPGE as two types of interlayer coordinate shift weighted by the quantum metric tensor and the Berry curvature, respectively. As a concrete example, we calculate the LCPGE in twisted bilayer graphene, and find that it exhibits a resonance peak whose frequency can be tuned from visible to infrared as the twisting angle varies. The LCPGE thus provides a promising route toward frequency-sensitive, circularly polarized light detection, particularly in the infrared range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Interactions of Oxide Surfaces with Water Revealed with Solid-State NMR Spectroscopy

Hydrous materials are ubiquitous in the natural environment and efforts have previously been made to investigate the structures and dynamics of hydrated surfaces for their key roles in various chemical and physical applications, with the help of theoretical modeling and microscopy techniques. However, an overall atomic-scale understanding of the water–solid interface, including the effect of water on surface ions, is still lacking. Herein, we employ ceria nanorods with different amounts of water as an example and demonstrate a new approach to explore the water–surface interactions by using solid-state NMR in combination with density functional theory. NMR shifts and relaxation time analysis provide detailed information on the local structure of oxygen ions and the nature of water motion on the surface: the amount of molecularly adsorbed water decreases rapidly with increasing temperature (from room temperature to 150 °C), whereas hydroxyl groups are stable up to 150 °C, and dynamic water molecules are found to instantaneously coordinate to the surface oxygen ions. Additionally, the applicability of dynamic nuclear polarization for selective detection of surface oxygen species is also compared to conventional NMR with surface selective isotopic-labeling: the optimal method depends on the feasibility of enrichment and the concentration of protons in the sample. These results provide new insight into the interfacial structure of hydrated oxide nanostructures, which is important to improve performance for various applications.

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

Quantifying Outer- and Inner-Coordination Sphere Effects Using Uranium Redox Chemistry in Molten Salt Solutions

Defining the relative influence of intramolecular and intermolecular forces is a fundamental problem in chemistry that is difficult to quantify. To address this challenge, we developed a method to evaluate the relative impact of direct chemical bonding in the inner-coordination sphere vs effects from cations in the outer-coordination sphere by comparative analysis of uranium redox reactivity in various molten salts. We observed that outer-coordination sphere cations (M 1+ ) and inner-coordination sphere anions (X 1– ) both affected uranium redox reactivity, with more polarizing M 1+ and larger X 1– favoring uranium in low oxidation states. Changing M 1+ (Li, Na, K) shifted the U IV + e 1– ⇌ U III (U IV/III ) and U III + 3e 1– → U 0 metal potentials by +330 and +240 mV, respectively. Changing X 1– (Cl, Br, I) caused larger shifts of +440 mV for the U IV/III redox potential and +1060 mV for the U 0 metal deposition potential. Using Coulomb’s Law, we correlated these potentials with electrostatic interactions between UIII and the molten salt. This model provided a facile way of predicting redox chemistry within molten salts.

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

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

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

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

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

Unprecedented lattice volume expansion on doping stereochemically active Pb 2+ into uniaxially strained structure of CaBa 1- x Pb x Zn 2 Ga 2 O 7

Lone pair cations like Pb 2+ are extensively utilized to modify and tune physical properties, such as nonlinear optical property and ferroelectricity, of some specific structures owing to their preference to adopt a local distorted coordination environment. Here we report that the incorporation of Pb 2+ into the polar “114”-type structure of CaBaZn 2 Ga 2 O 7 leads to an unexpected cell volume expansion of CaBa 1- x Pb x Zn 2 Ga 2 O 7 (0 ≤ x ≤ 1), which is a unique structural phenomenon in solid state chemistry. Structure refinements against neutron diffraction and total scattering data and theoretical calculations demonstrate that the unusual evolution of the unit cell for CaBa 1- x Pb x Zn 2 Ga 2 O 7 is due to the combination of the high stereochemical activity of Pb 2+ with the extremely strained [Zn 2 Ga 2 O 7 ] 4- framework along the c -axis. The unprecedented cell volume expansion of the CaBa 1- x Pb x Zn 2 Ga 2 O 7 solid solution in fact is a macroscopic performance of the release of uniaxial strain along c -axis when Ba 2+ is replaced with smaller Pb 2+ .

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

Electronic reconstruction at the polar (111)-oriented oxide interface

Atomically flat (111) interfaces between insulating perovskite oxides provide a landscape for new electronic phenomena. For example, the graphene-like coordination between interfacial metallic ion layer pairs can lead to topologically protected states [Xiao et al., Nat. Commun. 2, 596 (2011) and A. Rüegg and G. A. Fiete, Phys. Rev. B 84, 201103 (2011)]. The metallic ion/metal oxide bilayers that comprise the unit cell of the perovskite (111) heterostructures require the interface to be polar, generating an intrinsic polar discontinuity [Chakhalian et al., Nat. Mater. 11, 92 (2012)]. Here, we investigate epitaxial heterostructures of (111)-oriented LaAlO 3 /SrTiO 3 (LAO/STO). We find that during heterostructure growth, the LAO overlayer eliminates the structural reconstruction of the STO (111) surface with an electronic reconstruction, which determines the properties of the resulting two-dimensional conducting gas. This is confirmed by transport measurements, direct determination of the structure and atomic charge from coherent Bragg rod analysis, and theoretical calculations of electronic and structural characteristics. Interfacial behaviors of the kind discussed here may lead to new growth control parameters useful for electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation) [SWR-26-095]

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation): Multifidelity aerodynamic polar data generation for hydrofoil/tidal-turbine airfoil sections. Foilpolars ties together three pieces: *AeroSandbox supplies the baseline airfoil coordinates (UIUC database). *G2Aero parameterizes those shapes on a Grassmannian manifold (Karcher mean + PGA basis) and samples new perturbed shapes around that basis. *XFoil (panel method) and NeuralFoil (neural-network surrogate, shipped with AeroSandbox) each solve the resulting shapes for lift, drag, moment, and pressure at the swept angles of attack, Reynolds numbers, and n_crit values. Design optimization of foil shapes in a computationally efficient way requires polars data across many candidate shapes, not just a handful of baseline foils. However, high-fidelity CFD at that scale is too costly, and naive shape perturbation strays from realistic geometries. FOILPOLARS addresses this by loading baseline airfoils (via AeroSandbox) and mapping them onto a Grassmannian manifold (via G2Aero), computing a Karcher mean and principal geodesic analysis (PGA) basis. New shapes are sampled by perturbing PGA coefficients, keeping them close to the manifold of realistic foils. Each sampled shape is evaluated across a configurable sweep of angle of attack, Reynolds number, and critical amplification factor using two solvers: XFoil (panel method) and NeuralFoil (neural-network surrogate), producing a paired dataset of lift, drag, moment, pressure, convergence, and confidence, indexed alongside each shape's PGA coefficients and shared Grassmannian basis in a single xarray dataset. From this, FOILPOLARS produces convergence summaries and comparison plots per shape, Reynolds number, and n_crit. A command-line interface exposes each pipeline stage independently, supporting data-driven design, optimization, and machine-learning workflows for foils.

Sandhu, Rimple [National Laboratory of the Rockies↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Direct Observation of Circularly Polarized Nonlinear Optical Activities in Chiral Hybrid Lead Halides

Circularly polarized light emission is a crucial application in imaging, sensing, and photonics. However, utilizing low-energy photons to excite materials, as opposed to high-energy light excitation, can facilitate deep-tissue imaging and sensing applications. The challenge lies in finding materials capable of directly generating circularly polarized nonlinear optical effects. In this study, we introduce a chiral hybrid lead halide (CHLH) material system, R/S-DPEDPb 3 Br 8 ·H 2 O (DPED = 1,2-diphenylethylenediammonium), which can directly produce circularly polarized second harmonic generation (CP-SHG) through linearly polarized infrared light excitation, exhibiting a polarization efficiency as high as 37% at room temperature. To understand the spin relaxation mechanisms behind the high polarization efficiency, we utilized two models, so-called D’yakonov–Perel’ (DP) and Bir–Aronov–Pikus (BAP) mechanisms. The unique zigzag inorganic frameworks within the hybrid structure are believed to reduce the dielectric confinement and exciton binding energy, thus enhancing spin polarization, especially in regions with a high excitation pump fluence based on the DP mechanism. In the case of low excitation pump fluence, the BAP mechanism dominates, as evidenced by the observed decrease in the polarization ratio from CP-SHG measurement. Using density functional theory analysis, we elucidate how the distinctive 8-coordination environment of lead bromide building blocks effectively suppresses spin–orbit coupling at the conduction band minimum. This suppression significantly diminishes spin-splitting, thereby slowing the spin relaxation rate.

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

Metastable piezoelectric group-IV monochalcogenide monolayers with a buckled honeycomb structure

Multiple two-dimensional materials are being naïvely termed stable on the grounds of displaying phonon dispersions with no negative frequencies and of not collapsing on molecular dynamics calculations at fixed volume. But, if these phases do not possess the smallest possible structural energy, how does one understand and establish their actual meta stability? To answer this question, twelve two-dimensional group-IV monochalcogenide monolayers (SiS, SiSe, SiTe, GeS, GeSe, GeTe, SnS, SnSe, SnTe, PbS, PbSe, and PbTe) with a buckled honeycomb atomistic structure—belonging to symmetry group P3m1—displaying an out-of-plane intrinsic electric polarization are shown to be metastable by three independent methods. First, we uncover a coordination-preserving structural transformation from the low-buckled honeycomb structure onto the lower-energy Pnm2 1 (or Pmmn for PbS, PbSe, and PbTe) phase to estimate energy barriers E B that must be overcome during such structural transformation. Using the curvature of the local minima and E B as inputs to Kramers escape formula, large escape times are found, implying the structural metastability of the buckled honeycomb phase (with the exception of PbS and PbSe, these phases display escape times ranging from 700 years to multiple times the age of the universe and can be considered “stable” for practical purposes in that relative sense). The second demonstration is provided by phonon dispersion relations that include the effect of long-range Coulomb forces and display no negative vibrational modes. The third and final demonstration of structural metastability is furnished by room-temperature ab initio molecular dynamics for selected compounds. Here, the magnitude of the electronic band gap evolves with chemical composition. Different from other binary two-dimensional compounds such as transition metal dichalcogenide monolayers and hexagonal boron nitride monolayers which only develop an in-plane piezoelectric response, the twelve group-IV monochalcogenide monolayers with a buckled honeycomb structure also display out-of-plane piezoelectric properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

QM Investigation of Rare Earth Ion Interactions with First Hydration Shell Waters and Protein-Based Coordination Models

Here, conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln 3+ ) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln 3+ interaction with representatives of solvent (water, H 2 O) and protein binding sites (acetate, CH 3 COO – ), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln 3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln 3+ interactions with (H 2 O) 9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu 3+ , which bound less favorably than the preceding element, Yb 3+ . Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln 3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln 3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗