Search NASA⌕ Search

SEARCH · Search NASA

Results for “Regularization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Achieving geometric accuracy in FFT-based micromechanical models using conformal grid

Owing to its efficiency, simplicity and robustness, the FFT-based method has become the standard for computation of mechanical fields in a heterogeneous periodic unit cell. One of the main disadvantages of the FFT-based method is the inaccurate representation of the initial microstructure on a regular grid of voxels, which can be alleviated through the use of distorted initial grids. Here, in this paper, a method for generation of distorted initial grids conforming to the microstructural features (e.g. straight/curved boundaries) is proposed. The method determines the positions of the grid nodes in the initial configuration by solving a system of springs connecting the nodes. Microstructures consisting of layers, Voronoi tessellation and circular/spherical inclusions are considered, and mechanical fields simulated using the FFT-based method. It is found that distorted initial grids, conforming to the microstructural features, lead to more accurate mechanical fields in comparison to the corresponding non-distorted initial grid solution. The effect of initial grid distortion on the convergence of the FFT-based method is analyzed and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Jensen–Shannon divergence based novel loss functions for Bayesian neural networks

Bayesian neural networks (BNNs) are state-of-the-art machine learning methods that can naturally regularize and systematically quantify uncertainties using their stochastic parameters. Kullback–Leibler (KL) divergence-based variational inference used in BNNs suffer from unstable optimization and challenges in approximating light-tailed posteriors due to the unbounded nature of the KL divergence. To resolve these issues, we formulate a novel loss function for BNNs based on a new modification to the generalized Jensen–Shannon (JS) divergence, which is bounded. In addition, we propose a Geometric JS divergence-based loss, which is computationally efficient since it can be evaluated analytically. We found that the JS divergence-based variational inference is intractable, and hence employed a constrained optimization framework to formulate these losses. Our theoretical analysis and empirical experiments on multiple regression and classification data sets suggest that the proposed losses perform better than the KL divergence-based loss, especially when the data sets are noisy or biased. Specifically, there are approximately 5% and 8% improvements in accuracy for a noise-added CIFAR-10 dataset and a regression dataset, respectively. There is about 13% reduction in false negative predictions of a biased histopathology dataset. Additionally, we quantify and compare the uncertainty metrics for the regression and classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Large radial shift experiments in RHIC and their implications for EIC design

The Hadron Storage Ring (HSR) in the future Electron-Ion Collider (EIC) must operate over a broad range of design circumferences. In 2018 preliminary beam studies on the circumference adjustment capabilities of the Relativistic Heavy Ion Collider (RHIC) were performed by applying a small momentum offset to the circulating bunches without adjusting any bending magnets. The off-momentum linear optics were corrected back to on-momentum conditions. Applying a similarly small deviation to the dipole fields of a select set of bending magnets provides a large radial shift over much of the RHIC (or HSR) circumference while leaving the design trajectory unchanged in the insertion regions. Here, this paper presents the design of the different lattice configurations foreseen as the most viable options for the required HSR circumference changes, and highlights the modifications necessary for regular operations and to allow for testing these new settings in RHIC. Experimental results from 2021 and 2022 are reviewed and compared to model predictions obtained from both MAD-X and Bmad. The implications of these results for HSR design are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

The porosity surrounding carbides and second phase stringers in monolithic U-10Mo fuel plate after irradiation

Post-irradiation microstructure characterization plays an important role in qualifying the low-enriched uranium (LEU) monolithic U-10 wt%Mo plate-type fuel for United States high-performance research reactors (USHPRRs) program. Inhomogeneous features resulting from manufacturing and irradiation processes, including carbides, second phase stringers, and extensive void spaces caused by the combining of small porosities, may increase the risk of heat concentration in local regions of the fuel plate over the operating conditions. In this study, characteristics of carbides, stringers, and porosity after multiple levels of irradiation at varying fission densities were studied by electron microscopes to decipher the morphology of pores and the porosity evolution in U-10 wt%Mo. For carbides, the result shows that porosities start forming on UMo grain boundaries, then on UMo/carbides interfaces as the burn-up going higher. However, the porosities surrounding carbides grow larger than the ones on UMo grain boundaries. The porosities around the uranium carbides could interconnect to form larger void space. The study revealed that the void spaces larger than 5 µm were found around uranium carbides after high burnup, while no evidence was observed to support the similar voids formed near second phase stringers even though the size of the stringers (> 50 µm) was much larger than uranium carbides (< 20 µm). The evolution of porosities suggests that the formation of second phase stringers may not create more significant porosities compared to regular uranium carbides regions during fuel operating conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Unraveling plant phenotype to genotype associations with daily hyperspectral traits in Populus trichocarpa

Hyperspectral remote sensing is a powerful, high-throughput phenotyping tool that quantifies physiologically and structurally relevant wavelengths across diverse genotypes and over varying temporal scales. In this study, we combined tower-based continuous hyperspectral sensing with genome-wide association studies to analyze 1423 wavebands (400-900 nm) and derivative vegetation indices across 505 genotypes and the genetic architecture of hyperspectral phenotypes over time in Populus trichocarpa Torr. & Gray grown under field conditions. Wavelengths related to chlorophyll and carotenoid absorption spectra exhibited the strongest genetic variation resulting in 98 significant SNP associations. Notably, we found substantial overlap in genetic association between the blue and red spectral regions, indicative of carotenoids and chlorophyll, respectively, and identified more than 10 candidate genes associated with chloroplast function, underpinning photosynthetic activity. Furthermore, fluctuations in associations for vegetative indices, such as the chlorophyll:carotenoid index (CCI), across the growing season reveal a temporally dynamic genetic architecture of physiological traits associated with fall senescence of this temperate tree species. Finally, we also observed correlations (spearman rho = 0.3, p < 1x10 −8 ) between individual wavebands or vegetative indices and growth rate, assessed as the relative change of tree height over the growing season. The growth rate prediction was substantially improved by a regularization multivariate model (spearman rho>0.5, p < 1x10 −16 ), reinforcing the value of hyperspectral measurements for predicting traits linked to tree productivity. These findings highlight the potential of high-throughput, rapid, hyperspectral genome wide association studies GWAS to uncover physiologically meaningful genetic variation and offer promising insights for future acceleration for plant breeding.

09 BIOMASS FUELS↗
DOE OSTI2026

PyHydroGeophysX: An extensible open-source platform for integrating hydrological models with geophysical measurements

Hydrological models and geophysical measurements are widely used tools for understanding subsurface hydrological processes relevant to water resource management, yet they typically remain disconnected due to technical barriers. We present PyHydroGeophysX, an open-source Python platform bridging this gap by providing standardized interfaces between hydrological modeling software (MODFLOW, ParFlow) and geophysical simulation tools (PyGIMLi, SimPEG). The platform implements bidirectional workflows: translating hydrological outputs into simulated geophysical responses through petrophysical models, and extracting hydrological information from geophysical inversions. Key features include bidirectional workflow modules, configurable petrophysical models, time-lapse inversion with temporal regularization, parallel computing, and mesh utilities for property transfer between geophysical and hydrological grids. The modular architecture of PyHydroGeophysX enables researchers to incorporate additional models and methods, fostering broader adoption of integrated hydrogeophysical approaches. The software is freely available on GitHub and is intended for researchers and practitioners working at the intersection of hydrology and geophysics.

Hydrogeophysics↗
DOE OSTI2024

Calcium-organic matter fouling in nanofiltration: Synchrotron-based X-ray fluorescence and absorption near-edge structure spectroscopy for speciation

Calcium (Ca)-enhanced organic matter (OM) fouling of nanofiltration (NF) membranes leads to reduced flux during desalination and requires frequent cleaning. Fouling mechanisms are not fully understood, which limits the development of targeted fouling control methods. This study employed synchrotron-based X-ray fluorescence (XRF) and X-ray absorption near-edge structure (XANES) spectroscopy to quantify the spatial distribution and mass of Ca deposition as well as changes in the Ca coordination environment characteristic of specific fouling mechanisms, respectively. Bench-scale filtration experiments were performed using feed solutions containing Ca and ten different types of organic matter (OM), as well as the common scalants, calcium carbonate (CaCO 3 ) and calcium sulfate (CaSO 4 ). Osmotic backwash (OB) was performed at regular intervals for fouling control. Ca-OM aggregation resulted in greater flux decline and lower flux recovery during OB than Ca conditioning of membranes followed by filtration of feed solution with OM. Linear combination fitting (LCF) of XANES absorption spectra from fouled membranes indicated that Ca-OM aggregation preferentially occurred for OM types that exhibited both high carboxylic group and negative charge density. Consequently, these OM types exhibited greater deposition of Ca and TOC on the membrane surface when compared to other OM types. For the coexistence of scalants and OM, Ca speciation within the fouling layer was characteristic of both Ca bound to the membrane (i.e. potential bridging, charge screening) as well as Ca-OM aggregation and deposition mechanisms, while a range of crystal polymorphs were observed to occur simultaneously. XRF and XANES represent powerful tools for the elucidation of NF fouling mechanisms by quantification of Ca deposition as well as Ca speciation. Fouling control methods should target OM types with high carboxyl group density and negative charge to neutralize or eliminate interactions with Ca.

42 ENGINEERING↗
DOE OSTI2024

Tribological behavior of hybrid Aluminum-TiB 2 metal matrix composites for brake rotor applications

Here, this research investigates the feasibility of hybrid aluminum metal matrix composites (MMC) incorporating TiB 2 particles for brake rotor applications. The composites were produced by incorporating both in-situ, submicron-sized TiB 2 particles and regular micron-sized TiB 2 powders via stir and squeeze casting techniques into A206 aluminum alloy matrix. Systematic adjustments in the fractions of in-situ and ex-situ TiB 2 particles were conducted to evaluate their impact on wear behavior and mechanisms. Combination of both particle types allowed composites with up to 10 vol% of reinforcements. Composites with higher proportions of ex-situ particles demonstrated increased wear resistance compared to those solely composed of in-situ particles, control specimens without TiB 2 , and conventional cast iron counterparts. The lowest wear rate for the hybrid composites sliding against phenolic brake pads was 1.1 × 10 -5 mm 3 /Nm, signifying a 3-fold reduction relative to cast iron sliding against the same pads. Wear analysis elucidated distinctive mechanisms within the hybrid composites, characterized by mild fragmented abrasive wear, adhesive wear, and plastic deformation, accompanied by the formation of an intermixed tribo-oxide layer.

42 ENGINEERING↗
DOE OSTI2025

Structural response reconstruction using a system-equivalent singular vector basis

Here, this paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.

42 ENGINEERING↗
DOE OSTI2024

Venturing into Unexplored Phase Space: Synthesis, Structure, and Properties of MgCo 3 B 2 Featuring a Rumpled Kagomé Network

MgCo 3 B 2 , a novel ternary boride in a previously unexplored phase space, was synthesized using the hydride route. In situ powder X-ray diffraction and DFT calculations aided in the discovery of this compound, whose structure was then determined by single-crystal X-ray diffraction. Like the closely related CeCo 3 B 2 , MgCo 3 B 2 crystallizes in centrosymmetric space group P6/mmm (a = 4.883(2) Å, c = 2.926(2) Å at 210 K, Z = 1). Unlike CeCo 3 B 2 , however, it adopts a disordered structure that features a rumpled Kagomé network of Co atoms, and Mg atoms fill the channels of a Co–B framework. Although the structural disorder leads to motifs that are similar to those observed in MgNi 3 B 2 and other related ternary borides, no evidence of an ordered superstructure was found by single-crystal X-ray diffraction or high-resolution powder X-ray diffraction. In the case of CeCo 3 B 2 , boron atoms occupy the center of regular Co 6 trigonal prisms; in MgCo 3 B 2 , boron atoms are shifted from the center of the prism to form B–B dimers with roughly the same length as those found in MgNi 3 B 2 . Magnetic susceptibility data exhibit an unusual temperature dependence that cannot be convincingly modeled by the modified Curie–Weiss equation, consistent with DFT calculations predicting a nonmagnetic ground state. Intrinsic susceptibility at 300 K is 1.42 × 10 –3 emu/mol Oe, which is comparable to that of paramagnetic YCo 3 B 2 and CeCo 3 B 2 with a similar structure and composition. Here, this study showcases the efficacy of combining several methodologies to discover new solids in unexplored phase spaces. This approach includes in situ PXRD data to monitor reactions of precursors upon heating, a diffusion-enhanced synthesis method, and DFT assessment of compound stability.

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

Decorated Honeycomb Lattice and Successive Magnetic Transitions of the Delafossite Derivative K 4 Ni 5 Te 3 O 16

Here, we report the synthesis and characterization of the delafossite derivative K 4 Ni 5 Te 3 O 16 , which hosts a decorated honeycomb Ni 2+ lattice. Single-crystal X-ray diffraction, electron microscopy, and neutron diffraction establish a noncentrosymmetric Pmm2 crystal structure that combines regular honeycomb units with elongated motifs formed by Ni 2+ trimers along extended edges. Magnetic susceptibility reveals strong antiferromagnetic interactions, while heat capacity measurements identify two successive magnetic transitions at 30 and 12 K. Neutron diffraction shows that these transitions correspond to a progression from partial to full long-range magnetic order on distinct Ni sites, reflecting competition among the various exchange interactions. The resulting magnetic point group, mm2.1′, inherits the lattice polarity and permits high-order magnetoelectric coupling. Consistent with this symmetry, a quadratic magnetic-field-induced polarization is experimentally observed below the magnetic ordering temperature, likely arising from exchange-driven magnetostriction coupled to the polar lattice. These results establish K 4 Ni 5 Te 3 O 16 as an interesting platform for engineering B-site ordering in delafossite derivatives to realize complex quantum magnetism and functional responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗
DOE OSTI2025

An Experimental and Theoretical Study of the Valence Shell Electronic Structure of Nitromethane

Vibrationally resolved CH 3 NO 2 and CD 3 NO 2 photoelectron spectra and angular distribution parameters, β, have been measured in the photon energy range 20–80 eV. This allows a comprehensive investigation of conflicting literature interpretations of the two outermost photoelectron bands, complemented by high-level calculations characterizing the three O lone-pair-based orbital ionizations expected in this region. Franck–Condon simulations allow the regular vibrational progressions observed at the beginning of the first and second bands to be assigned as predominantly modes v 5 and v 6 of the D 0 and D 2 cation states, respectively. Irregular structuring observed midband in both bands is attributed to avoided crossings between the adiabatic D 0 /D 1 and D 1 /D 2 states as identified in potential surface cuts taken along the v 5 and v 6 normal mode coordinates. The associated D 0 /D 1 and D 1 /D 2 conical intersections were located and characterized, allowing possible vibronic coupling and its impact on the photoelectron spectra of the D 0 , D 1 , and D 2 states to be discussed. This helps rationalize the apparently missing D 1 ionization in this region, although some contribution by this state may be identifiable in the β-parameter spectra. Theoretical analysis of the complete valence region spectrum identifies the early breakdown of the independent electron model beyond these first two bands.

Powis, Ivan [University of Nottingham (United King↗
DOE OSTI2025

Nonorthogonal Configuration Interaction for Singlet Fission: Beyond the Dimer

Non-orthogonal configuration interaction with fragment calculations are presented for a number of compounds that show singlet fission properties: (i) four perylene-diimide derivatives, (ii) crystalline pentacene and its (B,N)-substituted variant, and (iii) a regular and a distorted stack of three indolonaphthyridine molecules. The electronic couplings between the singlet excitonic states (S 1 ) and the singlet-coupled double triplet (T 1 T 1 ), the so-called singlet fission coupling, were computed from ensembles with two and three molecules, and except for some small deviations when charge transfer states were included, results are virtually the same. Ensembles of three molecules were used to study the mechanisms of triplet separation, double triplet diffusion, and singlet and triplet exciton diffusion. The calculations show that apart from the standard mechanism for the generation of two uncoupled triplet states (S 1 → T 1 T 1 → T 1 ...T 1 ), there are two other possible pathways: the direct generation from the singlet excitonic state (S 1 → T 1 ...T 1 ) and the process in which the excitonic state evolves in a superposition of T 1 T 1 and T 1 ...T 1 states. Furthermore, the electronic coupling for triplet diffusion is in general much smaller than for singlet diffusion.

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

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗
DOE OSTI2026

Domain-Dependent Electronic Properties of 2D Copper Boride Synthesized from Reversible Subsurface Diffusion on Cu(111)

Emerging boron-based 2D materials hold properties with potential applications ranging from electronic devices to catalysis and energy storage. Using physical vapor deposition, 2D copper borides are formed on Cu(111) at elevated temperatures. However, the as-prepared surfaces often contain boron clusters that compromise structural homogeneity, and the dynamics of their formation is not yet fully understood. Through in situ low-energy electron microscopy (LEEM), Auger electron spectroscopy (AES), and scanning tunneling microscopy/spectroscopy (STM/STS), we show that cluster-free 2D CuB x is synthesized via reversible subsurface diffusion of boron atoms. The resolved regular atomic structure aligns with the previously calculated Cu 8 B 14 model. In a rare case, possibly due to surface strain and a change of orientation, a second domain of CuB x was observed with a pronounced change in electronic properties, yet electronic states remain delocalized in both domains. In conclusion, these findings demonstrate the sensitivity of electronic properties to domain structure and highlight a potential opportunity to tune 2D CuB x for diverse applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Ion Transport in Charged Membranes: Linking Electric-Field-Driven Mechanisms to Pore Size via Perturbation Analysis

Ion-exchange membranes are a critical component in electrochemical systems. Nevertheless, the understanding and modeling of ion transport within these porous structures have been limited by particular complexity reductions, either ignoring the dimensionality of their porous network architecture or imposing geometric assumptions (i.e., overlapping double layers). Before addressing this morphology-transport gap, a framework that relates the driving forces of transport to the geometry of a single pore is required. In this work, our modeling domain consists of a two-dimensional single pore with charged walls, connecting two identical electrolyte reservoirs. Using the Poisson-Nernst-Planck equations and regular perturbation theory, we decouple the electric fields and analyze the driving forces of ion transport, specifically electromigration and induced electroosmosis within the pore. These processes are described as analytical functions of the interaction aspect ratio, ?, defined as the ratio of the pore radius to the Debye length. Using this parameter, our study (i) describes the interplay between electromigrative and electroosmotic mechanisms that set ionic conductivity, (ii) identifies a dimensionless group of intrinsic electrolyte properties that indicates the predominant driving force, and (iii) provides a qualitative, confinement-dependent perspective on selectivity in ion-conducting membranes.

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

Identifying Topological Defects in Lamellar Phases through Contour Analysis of Complex Wave Fields

Lamellar phases frequently contain structural imperfections that significantly affect their behaviors and properties. Our previous research successfully reconstructed real-space configurations of defective lamellar phases from diffuse scattering patterns, indicating the presence of phase vortices as a potential method for identifying topological defects disrupting the smectic ordering. Here, this report presents a mathematical framework using regularized wave fields to represent defective lamellar structures in real space. Phase singularities, resulting from the interference of random waves and indicating lamellar order disruption, are identified through a contour integral. These wave fields, derived from coherent scattering in reciprocal space, were validated via computational benchmarks analyzing small-angle neutron scattering data from AOT surfactant solutions, facilitating further statistical analysis of the defects. Our study highlights the potential to extract meaningful information about topological defects in lyotropic phases by inversely analyzing experimentally measured two-point static correlations. Our method allows for detailed structural analysis of various lyotropic phases, both particulate and nonparticulate, in their quiescent states and facilitates quantitative investigation of defects’ role in phase transitions. By integrating small-angle scattering, deep learning, and vortex tangle analysis, our comprehensive approach shows promise in addressing complex challenges in the structural analysis of soft matter systems.

36 MATERIALS SCIENCE↗