SEARCH · Search NASA

Results for “Phase field crystal 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1993

The stability of hibonite and other aluminous phases in silicate melts: Implications for the origin of hibonite-bearing inclusions

Phase fields in which hibonite (Hib) and silicate melt coexist with spinel (Sp), CaAl4O7 (CA2), gehlenitic melilite (Mel), anorthite (An), or corundum (Cor) in the system CaO-MgO-Al203-SiO2-TiO2 (CMAST) were determined and activity models developed for Mel and Hib solid solutions. Experimentally determined partition coefficients for Ti between Hib and coexisting melt, D sub t, vary from 0.8 to 2.1 and generally decrease with increasing TiO2 content in the liquid (L). Based on Ti partioning between Hib and melt, bulk inclusion compositions and Hib-saturated liquid use phase diagrams, the Hib in Fluffy Type A inclusions (FTA's) from Allende and at least some of the Hib from Hib-rich inclusions is relict; much of the Hib from Hib-glass spherules probably crystallized from a melt under nonequilibrium conditions. Bulk compositions for all of these Ca-Al-rich inclusions (CAI's) are consistent with an origin as Mel + Hib + Sp + perovskite (Pv) proto-inclusions in which Mel was partially altered. In some cases, the proto-inclusion was partially or completely melted with vaporization occurring over a period of time sufficient to remove any Na introduced by the alteration process but frequently insufficient to dissolve all of the original hibonite. If equilibration temperatures based on Hib-bearing CAI's reflect condensation in a cooling gas of solar composition, then Hib + Cor condensed at approximately 1260 C (referenced to 10 exp -3 atm) and Hib + Sp + Mel at approximately 1215 +/- 10 C. Simple thermochemical models for the substitution of trace elements into the Ca-site of meteoritic Hib suggest that virtually all Eu is divalent in early condensate Hibs but that Eu(2+)/Eu(3+) decreases by a factor of 20 or more during the course of condensation, primarily because the ratio is proportional to the partial pressure of Al, which decreases dramatically as aluminous phases condense. The relative sizes of Eu and Yb anomalies in meteoritic Hibs and CAI's may be influenced by this effect.

Beckett, J. R.↗

DOE OSTI2026

Excitation spectrum and low-temperature magnetism in the disordered defect-fluorite Ho2⁢Zr2⁢O7

In this work, we report on the thermomagnetic characterization and crystalline-electric field (CEF) energy scheme of the disordered defect-fluorite Ho2⁢Zr2⁢O7. This structural phase is distinguished by the coexistence of magnetic frustration and extensive disorder, with Ho3+ and Zr4+ sharing randomly the same 4⁢𝑎 site with even 50% occupancy, and an average 1/8 oxygen vacancy per unit cell. AC magnetic susceptibility measurements performed on powder samples down to 0.5 K revealed signs of slowing spin dynamics without glassy behavior, including a frequency dependent peak at ∼1K. Yet, no evidence for long-range magnetic order is found down to 200 mK in the specific heat. Inelastic neutron scattering measurements show a weak, low-lying CEF excitation around 2 meV, accompanied by a broad level centered at 60 meV. To fit our observations, we propose an approach to account for structural disorder in the crystal-field splitting of the non-Kramers Ho3+. Our model provides an explanation to the broadening of the high-energy, single-ion excitations and suggests that the zirconate ground-state wavefunction has zero magnetic moment. Through the breaking of local symmetry, structural disorder enables the magnetic response observed in Ho2⁢Zr2⁢O7, allowing the mixing of low-lying states at finite temperatures. Finally, we show that this scenario is in good agreement with the bulk properties reported in this work.

Gardner, Jason [ORNL] (ORCID:0000000278234072)↗

DOE OSTI2025

Observation of Unprecedented Fractional Magnetization Plateaus in a New Shastry-Sutherland Ising Compound

Geometrically frustrated magnetic systems, such as those based on the Shastry-Sutherland lattice (SSL), offer a rich playground for exploring unconventional magnetic states. The delicate balance between competing interactions in these systems leads to the emergence of novel phases. We present the characterization of Er 2 ⁢Be 2⁢ GeO 7 , an SSL compound with Er 3+ ions forming orthogonal dimers separated by nonmagnetic layers whose structure is invariant under the 𝑃⁢$\bar{4}21$⁢𝑚 space group. Neutron scattering reveals an antiferromagnetic dimer structure at zero field, typical of Ising spins on that lattice and consistent with the anisotropic magnetization observed. However, magnetization measurements exhibit fractional plateaus at 1/4 and 1/2 of saturation, in contrast to the expected 1/3 plateau of the SSL Ising model. By comparing the energy of candidate states with ground-state lower bounds we show that this behavior requires spatially anisotropic interactions, leading to an anisotropic Shastry-Sutherland Ising model symmetric under the 𝐶⁢𝑚⁢𝑚⁢2 space group. This anisotropy is consistent with the small orthorhombic distortion observed with single-crystal neutron diffraction. The other properties, including thermodynamics, which have been investigated theoretically using tensor networks, point to small residual interactions, potentially due to further couplings and quantum fluctuations. This study highlights Er 2 ⁢Be 2 ⁢GeO 7 as a promising platform for investigating exotic magnetic phenomena.

Yadav, Lalit [Duke University, Durham, NC (United ↗

DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗

NASA NTRS2020

Polarimetric Radar Characteristics of Simulated and Observed Intense Convective Cores for a Midlatitude Continental and Tropical Maritime Environment

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Z(dr) and larger specific differential phase K(dp) above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Z(dr) and K(dp) observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Zdr and Kdp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

Toshihisa Matsui↗

DOE OSTI2025

Collinear antiferromagnetic order in a quasi-two-dimensional triangular lattice compound DyNiAl 4 Ge 2

The two-dimensional magnetic triangular lattice, with rare-earth local moments coupled through long-range interactions and hosting strong spin-orbit coupling, serves as an ideal platform for exploring novel quantum phenomena. Here, we report the physical properties of DyNiAl 4 Ge 2 single crystals through comprehensive structural, magnetic, heat capacity, electrical transport, and neutron powder diffraction measurements. DyNiAl 4 Ge 2 crystallizes in a trigonal lattice with the space group $R\bar{3}m$, where Dy atoms form two-dimensional triangular lattice layers within the 𝑎𝑏 plane and stack along the 𝑐 axis. Magnetic susceptibility measurement reveals a magnetic phase transition at 𝑇 𝑁 = 9.0 K, further confirmed by the 𝜆-shaped peak in the heat capacity curve. The obtained temperature-magnetic field phase diagram is relatively simple compared to other isostructural compounds, suggesting the absence of magnetic frustration. This is further evidenced by neutron powder diffraction, which shows a collinear antiferromagnetic ground state with a propagation vector 𝒌 = (0, 0, 1.5), implying the absence of magnetic frustration in DyNiAl 4 Ge 2 . In this magnetic structure, Dy moments are aligned ferromagnetically within each triangular layer and arranged antiferromagnetically between adjacent layers along the 𝑐 axis. In conclusion, our results establish DyNiAl 4 Ge 2 as a model system for understanding the interplay between RKKY interaction and geometrical frustration in this family of compounds.

Wu, HengHeng [Chinese Academy of Sciences (CAS), S↗

DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

NASA NTRS1989

New measurements of the far ultraviolet scattering properties of interstellar dust

An analysis is presented for the spectra of the diffuse ultraviolet background taken during shuttle flight STS-61C (January 1986). Eight regions of the sky were observed for approximately 20 minutes each, using a spectrograph designed specifically to perform measurements of the UV background. The field of view was 3.8 deg x 8 min, with imaging along the slit to confine stellar contamination. The instrument featured a shutter mechanism to measure internal background during flight, a low-scatter holographically ruled diffraction grating, photon counting microchannel plate detectors, through baffling, and a crystal window to further attenuate stray light. The spectra covered the range 1400 to 1850 A and was binned in 50 A bands. The procedure for substracting the contribution of stars too faint to be detected as discrete sources during the observations (in general this represents a small fraction of the total intensity detected except at the longest wavelengths). A radiative transfer model used to interpret the data and set confidence intervals on the relevant parameters is described. It was found that the continuum component of the diffuse ultraviolet background arises primarily from two sources. One source is scattering of starlight by interstellar dust with an albedo of about 12 percent and a relatively isotropic phase function. A second source consists of about 150 photons/cm/sec/ster/A of extragalactic light which is attenuated by the dust in the galaxy. Although emission features possibly associated with molecular H2 are detected in one look direction, fluorescence of H2 is not a major contributor to the diffuse UV background, at least at galactic latitudes greater than -10 degrees.

Hurwitz, Mark↗

NASA NTRS1999

Development of Optical Diagnostic Techniques for Microgravity Materials Processing

Materials processing including crystal growth, either under a gravity environment on ground or a microgravity environment in space, involves complicated phenomena of fluid motions in gas or liquid phases as well as interaction of various species. To obtain important physical insight, it is very necessary to provide gross-field optical diagnostics for monitoring various physical properties. Materials processing inhibits easy access by ordinary instruments and thus characterizing gross-field physical properties is very challenging. Typical properties of importance can be fluid velocity, temperature, and species concentration for fluids, and surface topology and defects for solids. Observing surface grow rate during crystal growth is also important. Material microstructures, i.e., integrity of crystal structures, is strongly influenced by the existence of thermally-induced flow as well as local nucleation of particles during solidification, which may act in many detrimental ways. In both ground-based and microgravity experiments, the nature of product property changes resulting from three-dimensional fluid or particle motions need be characterized. Gross-field diagnostics is thus required to identify their effects on product defects and process deficiencies. The quantitative visualization techniques can also be used for validation of numerical modeling. For optical nonintrusive gross-field diagnostic techniques, two approaches were developed as summer projects. One optical approach allows us to provide information of species concentration and temperature for monitoring in real time. The other approach, that is, the concept which is formulated for detection of surface topography measurement can provide unprecedented spatial resolution during crystal growth.

Cha, Soyoung Stephen↗

NASA NTRS2000

The Life Cycle Application of Intelligent Software Modeling for the First Materials Science Research Rack

Marshall Space Flight Center (MSFC) has been funding development of intelligent software models to benefit payload ground operations for nearly a decade. Experience gained from simulator development and real-time monitoring and control is being applied to engineering design, testing, and operation of the First Material Science Research Rack (MSRR-1). MSRR-1 is the first rack in a suite of three racks comprising the Materials Science Research Facility (MSRF) which will operate on the International Space Station (ISS). The MSRF will accommodate advanced microgravity investigations in areas such as the fields of solidification of metals and alloys, thermo-physical properties of polymers, crystal growth studies of semiconductor materials, and research in ceramics and glasses. The MSRR-1 is a joint venture between NASA and the European Space Agency (ESA) to study the behavior of different materials during high temperature processing in a low gravity environment. The planned MSRR-1 mission duration is five (5) years on-orbit and the total design life is ten (IO) years. The MSRR-1 launch is scheduled on the third Utilization Flight (UF-3) to ISS, currently in February of 2003). The objective of MSRR-1 is to provide an early capability on the ISS to conduct material science, materials technology, and space product research investigations in microgravity. It will provide a modular, multi-user facility for microgravity research in materials crystal growth and solidification. An intelligent software model of MSRR-1 is under development and will serve multiple purposes to support the engineering analysis, testing, training, and operational phases of the MSRR-1 life cycle development. The G2 real-time expert system software environment developed by Gensym Corporation was selected as the intelligent system shell for this development work based on past experience gained and the effectiveness of the programming environment. Our approach of multi- uses of the simulation model and its intuitive graphics capabilities is providing a concurrent engineering environment for rapid prototyping and development. Operational schematics of the MSRR-1 electrical, thermal control, vacuum access, and gas supply systems, and furnace inserts are represented graphically in the environment. Logic to represent first order engineering calculations is coded into the knowledge base to simulate the operational behavior of the MSRR-1 systems. An example of engineering data provided includes electrical currents, voltages, operational power, temperatures, thermal fluid flow rates. pressures, and component status indications. These type of data are calculated and displayed at appropriate instrumentation points, and the schematics are animated to reflect the simulated operational status of the MSRR-1. The software control functions are also simulated to represent appropriate operational behavior based on automated control and response to commands received by the crew or ground controllers. The first benefit of this simulation environment is being realized in the high fidelity engineering analysis results from the electrical power system G2 model. Secondly, the MSRR-1 simulation model will be embedded with a hardware mock-up of the MSRR-1 to provide crew training on MSRR-1 integrated payload operations. G2 gateway code will output the simulated instrumentation values, termed as telemetry, in a flight-like data stream so that the crew has realistic and accurate simulated MSRR-1 data on the flight displays which will be designed for crew use. The simulation will also respond appropriately to crew or ground initiated commands, which will be part of normal facility operations. A third use of the G2 model is being planned; the MSRR-1 simulation will be integrated with additional software code as part of the test configuration of the primary onboard computer, or Master Controller, for MSRR-1. We will take advantage of the G2 capability to simulate the flight like data stream to test flight software responses and behavior. A fourth use of the G2 model will be to train the Ground Support Personnel that will monitor the MSRR-1 systems and payloads while they are operating aboard the ISS. The intuitive, schematic based environment will provide an excellent foundation for personnel to understand the integrated configuration and operation of the MSRR-1, and the anticipated telemetry feedback based on operational modes of the equipment. Expert monitoring features will be enhanced to provide a smart monitoring environment for the operators. These features include: (1) Animated, intuitive schematic-based displays which reflect telemetry values, (1) Real-time plotting of simulated or incoming sensor values, (3) High/Low exception monitoring for analog data, (4) Expected state monitoring for discrete data, (5) Data trending, (6) Automated malfunction procedure execution to diagnose problems, (7) Look ahead capability to planned MSRR-1 activities in the onboard timeline. And finally, the logic to calculate telemetry values will be deactivated, and the same environment will interface to the incoming data for the real-time telemetry stream to schematically represent the onboard hardware configuration. G2 will be the foundation for the real-time monitoring and control environment. In summary, our MSRR-1 simulation model spans many elements of the life cycle development of this project: Engineering Analysis, Test and Checkout, Training of Crew and Ground Personnel, and Real-time monitoring and control. By utilizing the unique features afforded by an expert system development environment, we have been able to synergize a powerful tool capable of addressing our project needs at every phase of project development.

Rice, Amanda↗

NASA NTRSDate not supplied

Thermodynamic and Diffusion Model Estimates on Metamorphic Temperatures and Timescales for Basaltic Eucrite GRA 98098

Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.

J S Gorce↗

DOE OSTI2025

Emergent Inductance from Chiral Orbital Currents in a Bulk Ferrimagnet

Here, we report the discovery of a new form of inductance in the bulk ferrimagnet Mn 3 ⁢Si 2 ⁢Te 6 , which features strong spin-orbit coupling, large magnetic anisotropy, and pronounced magnetoelastic interactions. Below its Curie temperature (𝑇 𝐶 ≈ 78 K), Mn 3 ⁢Si 2 ⁢Te 6 hosts chiral orbital currents (COC) that circulate within the crystal lattice and give rise to collective electronic behavior [Control of chiral orbital currents in a colossal magnetoresistance material, Nature (London) 611, 467 (2022).]. By applying a magnetic field along the hard 𝑐 axis and driving the system with low-frequency currents, we uncover a giant inductive response up to millhenry scale, originating from first-order reconfigurations of COC domains. These domains act as coherent mesoscopic inductive elements that resist reversal upon current reduction, producing a large electromotive force and sharply increasing voltage. This emergent inductance defies classical models, occurs without superconductivity or engineered nanostructures, and opens a new frontier in orbital-based quantum functionality and device concepts.

complex materials↗

DOE OSTI2025

Quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) of zwitterionic organic HCl salts

In this work, we benchmark quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) for determining the crystal structures of two zwitterionic organic HCl salts, L-ornithine HCl ( Orn ) and L-histidine HCl·H 2 O ( Hist ). These salts present an interesting challenge for QNMRX-CSP, as gas-phase geometry optimizations used to generate starting structures for the organic zwitterionic fragments fail to capture their correct solid-state geometries. To overcome this limitation, geometry optimizations using the COSMO water-solvation model are employed to generate initial structural models. Using this approach, QNMRX-CSP yields structural models of the two zwitterionic organic HCl salts that closely match experimentally determined crystal structures. In addition, the application of QNMRX-CSP to Hist represents a further step toward the de novo structural determination of solvated organic HCl salts, as Hist is the first benchmark system of this type to include a water molecule as a component of its crystal structure. This work is significant for its potential application to the structural determination of active pharmaceutical ingredients, which often feature complex organic components and solvated solid forms.

Fleischer, Carl H. [Florida State Univ., Tallahass↗

DOE OSTI2024

Comparison of Machine Learning Approaches for Prediction of the Equivalent Alkane Carbon Number for Microemulsions Based on Molecular Properties

The chemical properties of oils are vital in the design of microemulsion systems. The hydrophilic–lipophilic difference equation used to predict microemulsions’ phase behavior expresses the oils’ physiochemical properties as the equivalent alkane carbon number (EACN). The experimental determination of EACN requires knowledge of the temperature dependence of the microemulsion system and the effects of different surfactant concentrations. Thus, the experimental determination is time-intensive and tedious, requiring days to months for proper separations. Furthermore, the experiments require high purity of chemicals because microemulsions are sensitive to impurities. Our work focuses on the quick and reliable predictions of the EACN with machine learning (ML) models. Due to the immaturity of ML chemical predictions, we compare three graph neural networks (GNNs) and a gradient-boosted tree algorithm, known as XGBoost. The GNNs use the molecular structures represented as simplified molecular-input line-entry system (SMILES) codes for the initial input, which allows us to assess whether geometry optimization is necessary for reliable results. The XGBoost model also begins with the SMILES representations of the molecules but uses molecular descriptors instead of geometry optimizations. As a result, the best model tested (crystal graph convolutional neural network with Merck molecular force field-94) has an error of 1.15 EACN units of the true EACN for unknown data with the errors skewed toward zero and an R² score of 0.9

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Frustration-driven magnetic correlations in the spin-$\frac{5}{2}$ triangular lattice antiferromagnet RbFe⁢(HPO 3 ) 2

Here, a detailed study of the structural and magnetic properties of a spin-5/2 triangular lattice antiferromagnet RbFe⁢(HPO 3 ) 2 is presented using x-ray diffraction, magnetization, heat capacity, and 31 P nuclear magnetic resonance (NMR) experiments on a polycrystalline sample. The crystal structure features an equilateral triangular lattice of Fe 3+ ions. The thermodynamic measurements reveal the onset of a magnetic long-range order at 𝑇 N⁢1 ≃7.8K in zero field, followed by another low-temperature field-induced ordering at 𝑇 N⁢2 in higher fields. The transition at 𝑇 N⁢1 is further confirmed from the NMR spin lattice relaxation measurements. The value of the frustration ratio (𝑓≃7) implies moderate spin frustration in the compound. The 31 P NMR spectra exhibit two distinct spectral lines corresponding to two inequivalent phosphorus sites (P1 and P2), consistent with the crystal structure. The P1 site is strongly coupled with an isotropic hyperfine coupling of 𝐴$^{iso}_{hf}$ = 0.55⁢(2)⁢ T/𝜇 B while the P2 site is weakly coupled with 𝐴$^{iso}_{hf}$ = 0.25⁢(3)⁢ T/𝜇 B with the Fe 3+ ions. The magnetic susceptibility and NMR shift data are described well assuming a spin-5/2 isotropic triangular lattice antiferromagnetic model with an average exchange coupling of 𝐽/𝑘 B = 2.8⁢(2)⁢K. Below 𝑇 N⁢1 , the spectra evolve into a nearly rectangular powder pattern, indicating a commensurate antiferromagnetic type order. The 31 P spin-lattice relaxation rate well below 𝑇N⁢1 follows a 𝑇 3 temperature dependence, implying a two-magnon Raman scattering mechanism in the ordered state. Three well-defined phase regimes are clearly ascertained in the 𝐻−𝑇 phase diagram, reflecting a weak magnetic anisotropy in the compound.

Nagpal, V. [Indian Institute of Science Education ↗

NASA NTRS2002

Application of Powder Diffraction Methods to the Analysis of the Atomic Structure of Nanocrystals: The Concept of the Apparent Lattice Parameter (ALP)

The applicability of standard methods of elaboration of powder diffraction data for determination of the structure of nano-size crystallites is analysed. Based on our theoretical calculations of powder diffraction data we show, that the assumption of the infinite crystal lattice for nanocrystals smaller than 20 nm in size is not justified. Application of conventional tools developed for elaboration of powder diffraction data, like the Rietveld method, may lead to erroneous interpretation of the experimental results. An alternate evaluation of diffraction data of nanoparticles, based on the so-called 'apparent lattice parameter' (alp) is introduced. We assume a model of nanocrystal having a grain core with well-defined crystal structure, surrounded by a surface shell with the atomic structure similar to that of the core but being under a strain (compressive or tensile). The two structural components, the core and the shell, form essentially a composite crystal with interfering, inseparable diffraction properties. Because the structure of such a nanocrystal is not uniform, it defies the basic definitions of an unambiguous crystallographic phase. Consequently, a set of lattice parameters used for characterization of simple crystal phases is insufficient for a proper description of the complex structure of nanocrystals. We developed a method of evaluation of powder diffraction data of nanocrystals, which refers to a core-shell model and is based on the 'apparent lattice parameter' methodology. For a given diffraction pattem, the alp values are calculated for every individual Bragg reflection. For nanocrystals the alp values depend on the diffraction vector Q. By modeling different a0tomic structures of nanocrystals and calculating theoretically corresponding diffraction patterns using the Debye functions we showed, that alp-Q plots show characteristic shapes which can be used for evaluation of the atomic structure of the core-shell system. We show, that using a simple model of a nanocrystal with spherical shape and centro-symmetric strain at the surface shell we obtain theoretical alp-Q values which match very well the alp-Q plots determined experimentally for Sic, GaN, and diamond nanopowders. The theoretical models are defined by the lattice parameter of the grain core, thickness of the surface shell, and the magnitude and distribution of the strain field in the surface shell. According to our calculations, the part of the diffraction pattern measured at relatively low diffraction vectors Q (below 10/angstrom) provides information on the surface strain, whle determination of the lattice parameters in the grain core requires measurements at large Q-values (above 15 - 20/angstrom).

Palosz, B.↗

NASA NTRS1992

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗

NASA NTRS1994

The stability of hibonite, melilite and other aluminous phases in silicate melts: Implications for the origin of hibonite-bearing inclusions from carbonaceous chondrites

Phase fields in which hibonite and silicate melt coexist with spinel CaAl4O7, gehlenitic melilite, anorthite or corundum at 1 bar in the system CaO-MgO-Al2O3-SiO2-TiO2 were determined. The hibonites contain up to 1.7 wt% SiO2. For TiO2, the experimentally determined partition coefficients between hibonite and coexisting melt D(sub i)(sup Hib/L), vary from 0.8 to 2.1 and generally decrease with increasing TiO2 in the liquid. Based on Ti partitioning between hibonite and melt, bulk inclusion compositions and hibonite-saturated liquidus phase diagrams, the hibonite in hibonite-poor fluffy Type A inclusions from Allende and at least some hibonite from hibonite-rich inclusions is relict, although much of the hibonite from hibonite-glass spherules probably crystallized metasably from a melt. Bulk compositions for all of these CAIs are consistent with an origin as melite + hibonite + spinel + perovskite phase assembalges that were partially altered and in some cases partially or completely melted. The duration of the melting event was sufficient to remove any Na introduced by the alteration process but frequently insufficient to dissolve all of the original hibonite. Simple thermochemical models developed for meteoritic melilite and hibonite solid solutions were used to obtain equilibration temperatures of hibonite-bearing phase assemblages with vapor. Referenced to 10(exp -3) atm, hibonite + corundum + vapor equilibrated at approximately 1260 C and hibonite + spinel +/- melilite + vapor at 1215 +/- 10 C. If these temperatures reflect condensation in a cooling gas of solar composition, then hibonite +/- corundum condensed first, followed by spinel and then melilite. The position of perovskite within this sequence is uncertain, but it probably began to condense before spinel. This sequence of phase appearances and relative temperatures is generally consistent with observed textures but differs from expectations based on classical condensation calculations in that equilibration temperatures are generally lower than predicted and melilite initially condenses with or even after spinel. Simple thermochemical modes for the substitution of trace elements into the Ca site of meteoritic hibonites suggest that virtually all Eu is divalent in early condensate hibonites but that Eu(2+)/Eu(#+) decreases by a factor of 20 or more during the course of condensation primarily because the ratio is proportional to the partial pressure of Al, which decreases dramatically as aluminous phase condense. The relative sizes of Eu and Yb anomalies in meteoritic hibonites and inclusions may be partly due to this effect.

Beckett, J. R.↗