Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pairwise interaction model”

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 73 records · Page 4

Chemical order transitions within extended interfacial segregation zones in NbMoTaW

Interfacial segregation and chemical short-range ordering influence the behavior of grain boundaries in complex concentrated alloys. In this study, we use atomistic modeling of a NbMoTaW refractory complex concentrated alloy to provide insight into the interplay between these two phenomena. Hybrid Monte Carlo and molecular dynamics simulations are performed on columnar grain models to identify equilibrium grain boundary structures. Here our results reveal extended near-boundary segregation zones that are much larger than traditional segregation regions, which also exhibit chemical patterning that bridges the interfacial and grain interior regions. Furthermore, structural transitions pertaining to an A2-to-B2 transformation are observed within these extended segregation zones. Both grain size and temperature are found to significantly alter the widths of these regions. An analysis of chemical short-range order indicates that not all pairwise elemental interactions are affected by the presence of a grain boundary equally, as only a subset of elemental clustering types are more likely to reside near certain boundaries. The results emphasize the increased chemical complexity that is associated with near-boundary segregation zones and demonstrate the unique nature of interfacial segregation in complex concentrated alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sensitivity of Wind Turbine Array Downstream Effects to the Parameterization Used in WRF

The Weather Research and Forecasting (WRF) Model has been extensively used for wind energy applications, and current releases include a scheme that can be applied to examine the effects of wind turbine arrays on the atmospheric flow and electricity generation from wind turbines. Herein we present a high-resolution simulation using two different wind farm parameterizations: 1) the “Fitch” parameterization that is included in WRF releases and 2) the recently developed Explicit Wake Parameterization (EWP) scheme. We compare the schemes using a single yearlong simulation for a domain centered on the highest density of current turbine deployments in the contiguous United States (Iowa). Pairwise analyses are applied to diagnose the downstream wake effects and impact of wind turbine arrays on near-surface climate conditions. On average, use of the EWP scheme results in small-magnitude wake effects within wind farm arrays and faster recovery of full WT array wakes. This in turn leads to smaller impacts on near-surface climate variables and reduced array–array interactions, which at a systemwide scale lead to summertime capacity factors (i.e., the electrical power produced relative to nameplate installed capacity) that are 2%–3% higher than those from the more commonly applied Fitch parameterization. It is currently not possible to make recommendations with regard to which wind farm parameterization exhibits higher fidelity or to draw inferences with regard to whether the relative performance may vary with prevailing climate conditions and/or wind turbine deployment configuration. However, the sensitivities documented herein to the wind farm parameterization are of sufficient magnitude to potentially influence wind turbine array siting decisions. Thus, our research findings imply high value in undertaking combined long-term high-fidelity observational studies in support of model validation and verification.

17 WIND ENERGY↗

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗

Four-dimensional chromosome reconstruction elucidates the spatiotemporal reorganization of the mammalian X chromosome

Significance Breakthrough technological advances such as Hi-C have enabled large-scale sequencing of pairwise genomic contacts. However, translating the 2D contacts to 3D spatial models remains challenging. Here we develop the “4DHiC” modeling method to infer 3D information from 2D datasets across time and construct a 3D model of the inactive X-chromosome (Xi). Our model not only demonstrates a spatial phase-separation between A/B and S1/S2 compartments to form “hemispheres,” but also reveals persistent smaller-scale structures hidden under the Xi superstructure. We also simulate Xist RNA spreading dynamics and deduce that the RNA spreads between hemispheres via both surface and core to establish the Xi. Our 4DHiC will serve a growing demand for data-driven methodologies that extract spatial information from biochemical interaction data.

Lappala, Anna↗

Planet Formation and the Characteristics of Extrasolar Planets

An overview of current theories of planetary growth, emphasizing the formation of extrasolar planets, is presented. Models of planet formation are based upon observations of the Solar System, extrasolar planets, and young stars and their environments. Terrestrial planets are believed to grow via pairwise accretion until the spacing of planetary orbits becomes large enough that the configuration is stable for the age of the system. Giant planets begin their growth like terrestrial planets, but if they become massive enough before the protoplanetary disk dissipates, then they are able to accumulate substantial amounts of gas. These models predict that rocky planets should form in orbit about most single stars. It is uncertain whether or not gas giant planet formation is common, because most protoplanetary disks may dissipate before solid planetary cores can grow large enough to gravitationally trap substantial quantities of gas. A potential hazard to planetary systems is radial decay of planetary orbits resulting from interactions with material within the disk. Planets more massive than Earth have the potential to decay the fastest, and may be able to sweep up smaller planets in their path. The implications of the giant planets found in recent radial velocity searches for the abundances of habitable planets are discussed.

Lissauer, Jack J.↗

Formation and Detection of Planetary Systems

Modern theories of star and planet formation and of the orbital stability of planetary systems are described and used to discuss possible characteristics of undiscovered planetary systems. The most detailed models of planetary growth are based upon observations of planets and smaller bodies within our own Solar System and of young stars and their environments. Terrestrial planets are believed to grow via pairwise accretion until the spacing of planetary orbits becomes large enough that the configuration is stable for the age of the system. Giant planets begin their growth as do terrestrial planets, but they become massive enough that they are able to accumulate substantial amounts of gas before the protoplanetary disk dissipates. These models predict that rocky planets should form in orbit about most single stars. It is uncertain whether or not gas giant planet formation is common, because most protoplanetary disks may dissipate before solid planetary cores can grow large enough to gravitationally trap substantial quantities of gas. A potential hazard to planetary systems is radial decay of planetary orbits resulting from interactions with material within the disk. Planets more massive than Earth have the potential to decay the fastest, and may be able to sweep up smaller planets in their path. The implications of the giant planets found in recent radial velocity searches for the abundances of habitable planets are discussed, and the methods that are being used and planned for detecting and characterizing extrasolar planets are reviewed.

Lissauer, Jack J.↗

Planet Formation

Modern theories of star and planet formation and of the orbital stability of planetary systems are described and used to discuss possible characteristics of undiscovered planetary systems. The most detailed models of planetary growth are based upon observations of planets and smaller bodies within our own Solar System and of young stars and their environments. Terrestrial planets are believed to grow via pairwise accretion until the spacing of planetary orbits becomes large enough that the configuration is stable for the age of the system. Giant planets begin their growth as do terrestrial planets, but they become massive enough that they are able to accumulate substantial amounts of gas before the protoplanetary disk dissipates. These models predict that rocky planets should form in orbit about most single stars. It is uncertain whether or not gas giant planet formation is common, because most protoplanetary disks may dissipate before solid planetary cores can grow large enough to gravitationally trap substantial quantities of gas. A potential hazard to planetary systems is radial decay of planetary orbits resulting from interactions with material within the disk. Planets more massive than Earth have the potential to decay the fastest, and may be able to sweep up smaller planets in their path. The implications of the giant planets found in recent radial velocity searches for the abundances of habitable planets are discussed, and the methods that are being used and planned for detecting and characterizing extrasolar planets are reviewed.

Lissauer, Jack J.↗

Planet Formation - Overview

Modern theories of star and planet formation are based upon observations of planets and smaller bodies within our own Solar System, exoplanets &round normal stars and of young stars and their environments. Terrestrial planets are believed to grow via pairwise accretion until the spacing of planetary orbits becomes large enough that the configuration is stable for the age of the system. Giant planets begin their growth as do terrestrial planets, but they become massive enough that they are able to accumulate substantial amounts of gas before the protoplanetary disk dissipates. These models predict that rocky planets should form in orbit about most single stars. It is uncertain whether or not gas giant planet formation is common, because most protoplanetary disks may dissipate before solid planetary cores can grow large enough to gravitationally trap substantial quantities of gas. A potential hazard to planetary systems is radial decay of planetary orbits resulting from interactions with material within the disk. Planets more massive than Earth have the potential to decay the fastest, and may be able to sweep up smaller planets in their path.

Lissauer, Jack J.↗

Designers' models of the human-computer interface

Understanding design models of the human-computer interface (HCI) may produce two types of benefits. First, interface development often requires input from two different types of experts: human factors specialists and software developers. Given the differences in their backgrounds and roles, human factors specialists and software developers may have different cognitive models of the HCI. Yet, they have to communicate about the interface as part of the design process. If they have different models, their interactions are likely to involve a certain amount of miscommunication. Second, the design process in general is likely to be guided by designers' cognitive models of the HCI, as well as by their knowledge of the user, tasks, and system. Designers do not start with a blank slate; rather they begin with a general model of the object they are designing. The author's approach to a design model of the HCI was to have three groups make judgments of categorical similarity about the components of an interface: human factors specialists with HCI design experience, software developers with HCI design experience, and a baseline group of computer users with no experience in HCI design. The components of the user interface included both display components such as windows, text, and graphics, and user interaction concepts, such as command language, editing, and help. The judgments of the three groups were analyzed using hierarchical cluster analysis and Pathfinder. These methods indicated, respectively, how the groups categorized the concepts, and network representations of the concepts for each group. The Pathfinder analysis provides greater information about local, pairwise relations among concepts, whereas the cluster analysis shows global, categorical relations to a greater extent.

Gillan, Douglas J.↗

Moving beyond the constraints of chemistry via crystal structure discovery with isotropic multiwell pair potentials

The rigid constraints of chemistry-dictated by quantum mechanics and the discrete nature of the atom-limit the set of observable atomic crystal structures. What structures are possible in the absence of these constraints? Here, we systematically crystallize one-component systems of particles interacting with isotropic multiwell pair potentials. We investigate two tunable families of pairwise interaction potentials. Our simulations self-assemble a multitude of crystal structures ranging from basic lattices to complex networks. Sixteen of the structures have natural analogs spanning all coordination numbers found in inorganic chemistry. Fifteen more are hitherto unknown and occupy the space between covalent and metallic coordination environments. The discovered crystal structures constitute targets for self-assembly and expand our understanding of what a crystal structure can look like.

crystal structures↗

Interaction–Deletion: A Composite Energy Method for the Optimization of Molecular Systems Selectively Removing Specific Nonbonded Interactions

The complex interactions between different portions of a large molecule can be challenging to analyze through traditional electronic structure calculations. Moreover, standard methods cannot easily quantify the physical consequences of individual pairwise interactions inside a molecule. By creating a set of molecular fragments, we propose a composite energy method to explore changes in a molecule caused by removing selected nonbonded interactions between different molecular portions. Energies and forces are easily obtained with this composite approach, allowing geometry optimizations that lead to chemically meaningful structures that describe how the omitted interactions contribute to changes in the local geometrical minima. We illustrate the application of our new hybrid scheme by computing the influence of intramolecular hydrogen-bonding interactions in two small molecules: 1,6-(tG + G + TG + G + g – )-hexanediol and a cyclic analogue, cis-1,4-cyclohexanediol. The resulting structural and energetic changes are interpreted to yield key physical insights and quantify concepts such as “preparation energy” or “reorganization energy”. We demonstrate that the composite method can be extended to larger molecular systems by showing its application on a Si(100) surface model containing interactions between dissociated ammonia molecules on adjacent surface dimers. The scheme’s efficacy is also tested by applying it to systems having multiple intramolecular interactions, viz., 310-polyglycine and H+GPGG. Furthermore, the cooperative nature of intramolecular hydrogen bonds is explored by using interaction–deletion in 2-nitrobenzene-1,3-diol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗

Bimetallic Copper/Ruthenium/Osmium Complexes: Observation of Conformational Differences Between the Solution Phase and Solid State by Atomic Pair Distribution Function Analysis

High-energy X-ray scattering and pair distribution function analysis (HEXS/PDF) is a powerful method to reveal the structure of materials lacking long-range order, but is underutilized for molecular complexes in solution. Here we demonstrate the application of HEXS/PDF with 0.26 Å resolution to uncover the solution structure of five bimetallic Cu I /Ru II /Os II complexes. HEXS/PDF of each complex in acetonitrile solution confirms the pairwise distances in the local coordination sphere of each metal center as well as the metal...metal distances separated by over 12 Å. The metal…metal distance detected in solution is compared with that from the crystal structure and molecular models to confirm that distortions to the metal bridging ligand are unique to the solid state. This work presents the first example of observing sub-Ångström conformational differences by direct comparison of solution phase and solid-state structures and shows the potential for HEXS/PDF in the determination of solution structure of single molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Origin of Ferromagnetic Exchange Coupling in Donor–Acceptor Biradical Analogues of Charge-Separated Excited States

A new donor–acceptor biradical complex, Tp Cum,Me Zn(SQ-VD) (Tp Cum,Me Zn + = zinc(II) hydro-tris(3-cumenyl-5-methylpyrazolyl)borate complex cation; SQ = orthosemiquinone; VD = oxoverdazyl), which is a ground-state analogue of a charge-separated excited state, has been synthesized and structurally characterized. Here, the magnetic exchange interaction between the S = 1/2 SQ and the S = 1/2 VD within the SQ-VD biradical ligand is observed to be ferromagnetic, with $J_{SQ-VD}$ = +77 cm –1 (H = -2$J_{SQ-VD}\bar{S}_{SQ}·\bar{S}_{VD}$) determined from an analysis of the variable-temperature magnetic susceptibility data. The pairwise biradical exchange interaction in Tp Cum,Me Zn(SQ-VD) can be compared with that of the related donor–acceptor biradical complex Tp Cum,Me Zn(SQ-NN) (NN = nitronyl nitroxide, S = 1/2), where $J_{SQ-NN}$ ≅ +550 cm –1 . This represents a dramatic reduction in the biradical exchange by a factor of ~7, despite the isolobal nature of the VD and NN acceptor radical SOMOs. Computations assessing the magnitude of the exchange were performed using a broken-symmetry density functional theory (DFT) approach. These computations are in good agreement with those computed at the CASSCF NEVPT2 level, which also reveals an S = 1 triplet ground state as observed in the magnetic susceptibility measurements. A combination of electronic absorption spectroscopy and CASSCF computations has been used to elucidate the electronic origin of the large difference in the magnitude of the biradical exchange coupling between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN). A Valence Bond Configuration Interaction (VBCI) model was previously employed to highlight the importance of mixing an SQ SOMO → NNLUMO charge transfer configuration into the electronic ground state to facilitate the stabilization of the high-spin triplet (S = 1) ground state in Tp Cum,Me Zn(SQ-NN). Here, CASSCF computations confirm the importance of mixing the pendant radical (e.g., VD, NN) LUMO (VD LUMO and NN LUMO ) with the SOMO of the SQ radical (SQ SOMO ) for stabilizing the triplet, in addition to spin polarization and charge transfer contributions to the exchange. An important electronic structure difference between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN), which leads to their different exchange couplings, is the reduced admixture of excited states that promote ferromagnetic exchange into the Tp Cum,Me Zn(SQ-VD) ground state, and the intrinsically weaker mixing between the VD LUMO and the SQ SOMO compared to that observed for Tp Cum,Me Zn(SQ-NN), where this orbital mixing is significant. The results of this comparative study contribute to a greater understanding of biradical exchange interactions, which are important to our understanding of excited-state singlet–triplet energy gaps, electron delocalization, and the generation of electron spin polarization in both the ground and excited states of (bpy)Pt(CAT-radical) complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗

Where, when, and why: Occurrence of fast-pairwise collective neutrino oscillation in three-dimensional core-collapse supernova models

Fast-pairwise collective neutrino oscillation represents a key uncertainty in the theory of core-collapse supernova (CCSN). Despite the potentially significant impact on CCSN dynamics, it is usually neglected in numerical models of CCSN because of the formidable technical difficulties of self-consistently incorporating this physics. In this paper, we investigate the prospects for the occurrence of fast flavor conversion by diagnosing electron neutrino lepton number (ELN) crossing in more than a dozen state-of-the-art three-dimensional CCSN models. ELN crossings provide a necessary condition for triggering flavor conversion. Although only zeroth and first angular moments are available from the simulations, our new method enables us to look into the angular distributions of neutrinos in momentum space and provide accurate insight into ELN crossings. Our analysis suggests that fast flavor conversion generally occurs in the postshock region of CCSNe, and that explosive models provide more favorable conditions for the flavor conversion than failed CCSNe. Here we also find that there are both common and progenitor-dependent characteristics. Classifying ELN crossings into two types, we analyze the generation mechanism of each case by scrutinizing the neutrino radiation field and matter interactions. We find key ingredients of CCSN dynamics driving the ELN crossings: proto-neutron star convection, asymmetric neutrino emission, neutrino absorptions, and scatterings. This study suggests that we need to accommodate fast flavor conversions in realistic CCSN models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗

Computational Optimization of Room Temperature Usable Capacity for Hydrogen Storage in MFU-4-Type Metal–Organic Frameworks via Pairwise Metal Substitutions

The efficient storage of hydrogen is a critical challenge in the quest for sustainable energy solutions. Current adsorbent-based methods achieve satisfactory storage densities predominantly under cryogenic temperatures and/or high pressures, which imposes problems with cost-efficient and safe implementation of this technology. Materials that can bind hydrogen gas reversibly at ambient temperatures and more moderate pressures could play a pivotal role in enabling hydrogen-powered technologies. In this study, we use reliable computational modeling to investigate two synthetically feasible paths for tuning the enthalpy of H2 binding in MFU-4-type metal–organic frameworks (MOFs), aiming to maximize usable capacity. This study examines MIM4 IICl3(bta)6 (bta– = benzotriazolate) Kuratowski-type clusters as a model for strong binding sites in MFU-4l frameworks. We systematically evaluate the impact of separately tuning the central MII metal ion (which plays a structural role) and the peripheral MI metal ion (which binds the substrate) on the energetics of H2 binding. Our computational study reveals that H2 binding at an MI site mostly follows the trend AgI < CuI < NiI < CoI < AuI while a larger central MII site generally weakens the H2 binding at a MI site. Importantly, we have identified three new combinations of MI and MII to achieve high fractional usable capacities of the total H2 adsorbed under a pressure swing from 5 to 100 bar at room temperature. Additionally, we examine the nature of the binding interaction between the peripheral metal atom and the hydrogen molecule. While charge transfer predominantly induces this interaction, for several atom combinations, a change in the polarization (associated with variations in the ionic radius of the MI binding atom) is another important factor for adjusting the strength of the interaction. We suggest that the proposed compositions of Kuratowski-type clusters are highly desirable synthetic targets for future laboratory study.

Tkachenko, Nikolay V↗