Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Active matter”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Multi-modal dynamic radiography using short-pulse laser-generated probe beams

Radiography is an important tool for the interrogation of dynamic experiments in the fields of dynamic properties of materials, and in condensed matter, high explosive, and high-energy-density physics. Multi-modal radiography advances the hypothesis that combining the information delivered by multiple radiographic modalities can lead to more constrained (improved) “reconstruction” of the scene than can be obtained from a single probe. We identify four modalities: multi-probe, time sequence, multi-view, and multi-messenger. Multi-probe radiography is a promising candidate for a next-generation dynamic radiographic facility. High-energy X-rays are the most frequently used probe for dynamic radiography, although recent developments show the utility of proton (pRad), electron (eRad), and neutron probe beams. Because each probing species interacts with material in the radiographic scene through quantitatively different mechanisms, each returns independent information about the scene, which can add extra constraints to the reconstruction process. How to conduct detailed, quantitative “co-analysis” of multiple data streams remains an area of active research. Multi-beam, short-pulse, laser-generated probes offer sufficient dose, an appropriate spectrum, and appropriate spatio-temporal resolution to produce high-quality dynamic radiographs. This paper reports on technology development to advance the state of the art of multi-modal/multi-probe radiography and the pursuit of both deterministic and inferential (AI/ML assisted) co-analysis methodologies to produce more constrained reconstructions from multi-modal data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of Galactic Irradiation on Thermal and Electronic Transport in Tungsten

The impact of irradiation on the thermal and electronic properties of materials is a persistent puzzle, particularly defect formation at the atomic and nanoscales. This work examines the nanoscale effects of low-energy irradiation on tungsten (W), focusing on defect-induced modifications to thermal and electronic transport. Using the site-projected thermal conductivity (SPTC) method [A. Gautam et al. PSS-RRL, 2400306, 2024], bulk and twin-grain boundary W with vacancy defects based on the Norgett–Robinson–Torrens displacements per atom model is analyzed. SPTC provides a detailed prediction of postcascade spatial thermal conductivity distribution. Electronic conductivity activity is estimated using the “N 2 method” [K. Nepal et al. Carbon, 119711, 2025] to explore the consequences of vacancies and grain boundaries, highlighting the defect-dependent nature of charge transport behavior. These findings offer high-resolution insights into irradiation-driven transport phenomena, with implications for space-exposed materials and nanoscale thermal/electronic management.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Extended defects-enhanced oxygen diffusion in ThO 2

Oxygen self-diffusion is key to understanding stoichiometry and defect structures in oxide nuclear fuels. Experimentally, low activation-barrier oxygen migration was found in ThO 2 , a candidate nuclear fuel, possibly due to short-circuit diffusion mechanisms. Here, we perform extensive molecular dynamics simulations to show that various types of extended defects can enhance oxygen self-diffusion with a much-reduced activation barrier in ThO 2 . In this work, we consider extended defects including 1D (dislocation), 2D (grain boundary), and 3D (void) defects. Due to the distinct characteristics of each type of extended defect, the modulation of oxygen diffusion varies. Further, these results provide a quantitative description of oxygen transport, which is significantly enhanced within a close distance (nanometer scale) from the extended defects. Among all these considered defects, grain boundary, particularly the low-energy $Σ$3 twin boundary, exhibits the strongest effect on increasing oxygen transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effect of bacterial cell addition on Fe(III) reduction and soil organic matter transformation in a farmland soil

The coupled biogeochemical processes of Fe(III) reduction and organic matter transformation profoundly impact terrestrial carbon cycling. However, little is known about the microbial role in soil organic matter (SOM) transformation during Fe(III) bio-reduction. Here we investigated the bio-reduction behavior of a black farmland soil and corresponding SOM transformation under circumneutral and anoxic conditions. A model dissimilatory Fe reducing bacterium, Geobacter sulfurreducens, was added in either live or dead form in order to enhance Fe(III) reduction in soil and to evaluate the accompanying transformation of SOM. The progress of Fe reduction was monitored and SOM transformation was characterized by various spectroscopy methods. Results showed that addition of either dead or live G. sulfurreducens cells increased the Fe(III) reduction rate and extent. Without cell addition, SOM transformation was insignificant within 13 days of incubation, only with some consumption of aliphatic/protein compounds, apparently due to their higher bio-degradability. Addition of dead or live cells resulted in more drastic SOM transformation, but through different mechanisms. With dead cell amendment, cell necromass and debris stimulated the activity of indigenous soil microbes by serving as extra carbon/energy sources. During Fe(III) reduction, the aliphatic/protein compounds were preferentially consumed by indigenous microbial communities, similar to the treatment without cell addition but with a greater extent of consumption. In comparison, addition of live G. sulfurreducens cells stimulated degradation of less bioavailable compounds, including more saturated and higher molecular weight SOM. It is possible that fast depletion of labile SOM by live G. sulfurreducens cells favored utilization of less bioavailable molecules (such as those with more aromatic structures) by the originally dormant species in native microbial community, suggesting an active role of live Fe(III)-reducing bacteria in affecting SOM transformation. In addition, fast assimilation of microbial related carbon, such as aromatic proteins and microbial byproducts, into SOM pools was also observed. Overall, our results suggest that addition of Fe(III)-reducing bacteria to a native soil can enhance Fe(III) reduction and accelerate the turnover of SOM through various mechanisms. The study provides new insights into coupled Fe(III) reduction and SOM transformation, as well as the “priming effect” of SOM using microbial cells as substrate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Identifying $\textit{f}$-electron symmetries of UTe 2 with O-edge resonant inelastic x-ray scattering

The recent discovery of spin-triplet superconductivity emerging from a nonmagnetic parent state in UTe 2 has stimulated great interest in the underlying mechanism of Cooper pairing. Experimental characterization of short-range electronic and magnetic correlations is vital to understanding these phenomena. Here we use resonant inelastic x-ray scattering (RIXS), x-ray absorption spectroscopy (XAS), and atomic-multiplet-based modeling to shed light on the active debate between $5f^26d^1$-based models with singlet crystal field states versus $5f^3$-based models that predict atomic Kramers doublets and much greater $5\textit{f}$ itinerancy. In this work, the XAS and RIXS data are found to agree strongly with predictions for a $5f^26d^1$-like valence electron configuration with weak intradimer magnetic correlations, and provide a context for interpreting recent investigations of the electronic structure and superconducting pairing mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The number of catalytic cycles in an enzyme’s lifetime and why it matters to metabolic engineering

Metabolic engineering uses enzymes as parts to build biosystems for specified tasks. Although a part’s working life and failure modes are key engineering performance indicators, this is not yet so in metabolic engineering because it is not known how long enzymes remain functional in vivo or whether cumulative deterioration (wear-out), sudden random failure, or other causes drive replacement. Consequently, enzymes cannot be engineered to extend life and cut the high energy costs of replacement. Guided by catalyst engineering, we adopted catalytic cycles until replacement (CCR) as a metric for enzyme functional life span in vivo. CCR is the number of catalytic cycles that an enzyme mediates in vivo before failure or replacement, i.e., metabolic flux rate/protein turnover rate. We used estimated fluxes and measured protein turnover rates to calculate CCRs for ~100–200 enzymes each from Lactococcus lactis, yeast, and Arabidopsis. CCRs in these organisms had similar ranges (<10 3 to >10 7 ) but different median values (3–4 × 10 4 in L. lactis and yeast versus 4 × 10 5 in Arabidopsis). In all organisms, enzymes whose substrates, products, or mechanisms can attack reactive amino acid residues had significantly lower median CCR values than other enzymes. Taken with literature on mechanism-based inactivation, the latter finding supports the proposal that 1) random active-site damage by reaction chemistry is an important cause of enzyme failure, and 2) reactive noncatalytic residues in the active-site region are likely contributors to damage susceptibility. Enzyme engineering to raise CCRs and lower replacement costs may thus be both beneficial and feasible.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-temperature phonon-mediated superconductivity in monolayer Mg 2 B 4 C 2

A two-dimensional material – Mg 2 B 4 C 2 , belonging to the family of the conventional superconductor MgB 2 , is theoretically predicted to exhibit superconductivity with critical temperature T c estimated in the 47–48 K range (predicted using the McMillian-Allen-Dynes formula) without any tuning of external parameters such as doping, strain, or substrate-induced effects. The origin of such a high intrinsic T c is ascribed to the presence of strong electron-phonon coupling and large density of states at the Fermi level. This system is obtained after replacing the chemically active boron-boron surface layers in a MgB 2 slab by chemically inactive boron-carbon layers. Hence, the surfaces of this material are inert. Our calculations confirm the stability of 2D Mg 2 B 4 C 2 . We also find that the key features of this material remain essentially unchanged when its thickness is increased by modestly increasing the number of inner MgB 2 layers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Vibrational properties of CuInP 2 S 6 across the ferroelectric transition

In order to explore the properties of a two-sublattice ferroelectric, we measured the infrared and Raman scattering response of CuInP 2 S 6 across the ferroelectric and glassy transitions and compared our findings to a symmetry analysis, calculations of phase stability, and lattice dynamics. In addition to uncovering displacive character and a large hysteresis region surrounding the ferroelectric transition temperature T C , we identify the vibrational modes that stabilize the polar phase and confirm the presence of two ferroelectric variants with opposite polarizations. Below T C , a poorly understood relaxational or glassy transition at T g is characterized by local structure changes in the form of subtle peak shifting and activation of low frequency out-of-plane Cu- and In-containing modes. Here, the latter are due to changes in the Cu/In coordination environments and associated order-disorder processes. Moreover, T g takes place in two steps with another large hysteresis region and significant underlying scattering. Combined with imaging of the room temperature phase separation, this effort lays the groundwork for studying CuInP 2 S 6 under external stimuli and in the ultrathin limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

An Accreting Supermassive Black Hole Buried in a Faint Dwarf Galaxy

Abstract In the last decade, there have been several discoveries of active galactic nuclei (AGN) in dwarf galaxies including an AGN in an ultracompact dwarf galaxy with a black hole mass >10 6 M ⊙ . However, finding a supermassive black hole (SMBH) in a dwarf low surface brightness (LSB) galaxy is rare. We report the discovery of a Seyfert type 2 class AGN that is associated with a nuclear SMBH of mass ∼6.5 × 10 6 M ⊙ in a dwarf LSB galaxy (μ 0,r > 23.8 mag arcsec −2 ) that we denote by MJ0818+2257. The galaxy was previously thought to be an outlying emission blob around the large spiral galaxy LEDA 1678924. In our current analysis, which includes the detection of the optical counterpart of MJ0818+2257, we study its ionized gas kinematics and find that the dynamical mass within the ionized gas disk is ∼5.3 × 10 9 M ⊙ . This is comparable to its stellar mass, which is ∼3 × 10 9 M ⊙ , and suggests that MJ0818+2257 is moderately dark matter dominated within the stellar disk. The SMBH-mass-to-galaxy-stellar-mass ratio isM BH /M(*) > 0.022, which is high compared to disk galaxies. Our detection of an SMBH in a bulgeless LSB dwarf galaxy raises questions about the growth of SMBHs in low-luminosity galaxies and suggests the possibility of detecting heavy seed black holes from early epochs in LSB dwarf galaxies in the low-redshift Universe.

Astronomy & Astrophysics↗
DOE OSTI2020

Hall-coefficient diagnostics of the surface state in pressurized SmB 6

In this study, we report results of the high-pressure Hall coefficient (R H ) measurements in the putative topological Kondo insulator SmB 6 up to 37 GPa. Below 10 GPa, our data reveal that R H (T) exhibits a prominent peak upon cooling below 20 K. Remarkably, the temperature at which surface conduction dominates coincides with the temperature of the peak in R H (T). The temperature dependent resistance and Hall coefficient can be well fitted by a two-channel model with contributions from the metallic surface and the thermally activated bulk states. When the bulk of SmB 6 becomes metallic and magnetic at ~10 GPa, both the R H (T) peak and the resistance plateau disappear simultaneously. Here, our results indicate that the R H (T) peak is a fingerprint to diagnose the presence of a metallic surface state in SmB 6 . The high-pressure magnetic state of SmB 6 is robust to 180 GPa, and no evidence of superconductivity is observed in the metallic phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

First-principles insight into all-optical spin switching in the half-metallic Heusler ferrimagnet Mn 2 RuGa

All-optical spin switching (AOS) represents a new frontier in magnetic storage technology—spin manipulation without a magnetic field—but its underlying working principle is not well understood. Many AOS ferrimagnets such as GdFeCo are amorphous and renders the high-level first-principles study unfeasible. The crystalline half-metallic Heusler Mn 2 RuGa presents an opportunity. Here we carry out hitherto the comprehensive density functional investigation into the material properties of Mn 2 RuGa, and introduce two concepts—the spin anchor site and the optical active site—as two pillars for AOS in ferrimagnets. In Mn 2 RuGa, Mn(4a) serves as the spin anchor site, whose band is below the Fermi level and has a strong spin moment, while Mn(4c) is the optical active site whose band crosses the Fermi level. Our magneto-optical Kerr spectrum and band structure calculation jointly reveal that the delicate competition between the Ru-4d and Ga-4p states is responsible for the creation of these two sites. Furthermore these two sites found here not only present a unified picture for both Mn 2 RuGa and GdFeCo, but also open the door for future applications. Specifically, we propose a Mn 2 Ru x Ga-based magnetic tunnel junction where a single laser pulse can control magnetoresistance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Excitons in periodic potentials

The energy band structure of excitons is studied in periodic potentials produced by the short-range interaction between the exciton and electrons of Wigner lattices. Treating the exciton as a pointlike dipole that interacts with the periodic potential, we can solve a simple one-body problem that provides valuable information on excitons in many-body problem settings. By employing group theory, we identify the excitonic energy bands that can couple to light and then quantify their energy shifts in response to a change in the period of the potential. This approach allows us to emulate the response of optically active exciton and trion states to a change in electron density. In agreement with experimental results, the exciton (trion) state shows a relatively strong (weak) energy shift when the electron density increases. This behavior stems from the difficulty of exciton states to remain extended, whereas the trion state remains largely the same as long as its radius is smaller than the average distance between resident electrons. Here, we identify an important relation between the energy blueshift of bright excitons and the presence of a Wigner crystal, where the blueshift magnitude is proportionate to the symmetry of the electron lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Particle Matter Index and Fuel Wall-wetting Relations on Stochastic Pre-ignition

This works explores the effect of the particle matter index (PMI) and aromatic content on fuel wall impingement and associated stochastic pre-ignition (SPI) propensity. Statically significant measurements of SPI rates are directly coupled with laser induced florescence (LIF) measurements of fuel dilution from spray-linear impingement. Literature suggests that PMI is could be correlated with the number of SPI events, but the root cause(s) of PMI and SPI are directly causational or are a predicator of SPI. Three fuels have been used in this study with 3 different PMI and two different aromatic contents. These fuels are direct injected at two different injection timings, an earlier injection timing which targets the piston crown, 310°CA bTDC, and a later injection timing that the linear, 220°CA bTDC start of injection timings (SOI) respectively. The earlier 310 SOI injection increases soot, whereas the later 220°CD SOI targets the linear and increases wall-wetting. The findings of this work highlight that the SPI activity has a weak correlation with the soot promoted by the combustion process and increased PMI. However, SPI activity shows a strong dependency with high levels of oil dilution and aromatic content, suggesting that PMI might be an indicator of fuel chemistry conducive to SPI and the soot form increased PMI fuels is not a strong source of SPI relative to fuel wall wetting.

42 ENGINEERING↗
DOE OSTI2023

Building a Computational and Experimental Rapid Response Pipeline to Counter the Coronavirus Disease 2019 Outbreak and Emerging Biothreats

The COVID-19 pandemic underscored the promise of monoclonal antibody-based prophylactic and therapeutic drugs, especially where protective candidates can be rapidly identified or developed for emerging biothreats and escape variants. Current cutting-edge technologies for this purpose still rely on pathogen-exposed convalescent volunteers and a large screening effort to find a proverbial needle in a haystack. Computational design of protective antibodies based on pre-existing templates skips those requirements and allows for greater control over the breadth and target epitope, while also co-optimizing for potency and developability or other biophysical characteristics. We approached this problem by building and expanding an in vitro experimental rapid antibody production and characterization pipeline to support development of an autonomous, closed loop, active learning software system based on structural simulation and ground truth experimental data to design and evaluate antibody antigen interactions. Starting from early in the pandemic, we targeted SARS-CoV-2, the causative agent of COVID-19, by re-purposing neutralizing antibodies against SARS-CoV-1 that had been identified in the wake of that outbreak in the early 2000’s. We successfully re-targeted three different anti-SARS-CoV-1 antibodies to neutralize SARS-CoV-2 in vitro, where the antibodies were generated externally and tested through conventional binding and neutralization assays internally or with collaborators. As antibodies were identified from the blood of humans infected with SARS-CoV-2, we shifted to collaborate with academic partners to develop improved versions of their human-derived antibodies. This work reached its most important stage in rapid response to the emergence of the Omicron variant of concern (VOC) in late 2021. In a matter of weeks, enabled by on demand innovation to our screening pipeline, we computationally designed and experimentally characterized derivative antibodies of COV2-2130, one of two antibodies from Vanderbilt that form the basis of the AstraZeneca Evusheld prophylactic drug product. This drug product suffers a serious loss of efficacy against Omicron BA.1 and BA.1.1, the first Omicron strains. Due to tight integration of computational design and experimental evaluation, we were able to identify a pair of designs with potent neutralization of the main targets Omicron BA.1 and BA.1.1; but also the earlier Delta variant, and subsequent Omicron strains including BA.2, BA.4, BA.5, and BA.2.75, demonstrating that our multi-target design process can, by its nature, produce robust antibody designs that strictly improve over the parental antibody. These results, recognized by a 2022 Director’s Science and Technology award, have enabled the follow-on GUIDE program, to commence in FY23. While earlier design campaigns were substantially outsourced, we have engineered better and faster processes internally to better compliment, calibrate, and speed computational designs. As part of the follow-on GUIDE program, we will stand up a rapid and high-throughput antibody production and characterization facility staffed with the expertise and capabilities to foster our current collaboration across PLS and ENG as well as other partnerships toward computational design of biologics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Security-by-Design: Light Water Small Modular Reactor

The growing demand for nuclear power is increasing pressure to find solutions to cost prohibitive requirements of both construction and security. Offsite response has been proposed as an option to reduce costs associated with training and maintaining an onsite response force. A previous report explored this option and revealed that security could be provided at the required level, but cost savings was not a result of this methodology. An offsite response strategy required costly active and passive delay barriers to provide sufficient time for responders to muster and deploy to a site in time to interrupt a determined and well-equipped adversary. Also, contrary to the hypothesis, the number of responders required for this strategy exceeded that needed for an onsite response force, as the adversaries could avail themselves of advantageous positions within the facility to repel arriving responders. This report builds upon the previous evaluation by using the same hypothetical light water small modular reactor (LWSMR) facility model, but this time an onsite response strategy was assessed. The goal of this analysis was to show that an onsite response strategy could be implemented effectively at a cost point that removes barriers within the industry at this critical time of growth and development. The assessment of the facility design and response strategy was completed through modeling using Scribe3D© and subsequent scenario analysis over the course of a two-day tabletop exercise. Subject matter experts in nuclear security, nuclear facility design, and response strategy and tactics contributed to the effort to ensure accurate representation of hypothetical scenarios. Several adaptations were made to the layout of the LWSMR based on lessons learned during the first day of scenario analysis. The subsequent design evaluated on the second day proved to provide a robust response posture against a large and well-trained adversary force. This report details the process of the analysis and compares the cost of the final facility design with that of the LWSMR model used for evaluation of offsite response. Ultimately, the results of this effort indicate that, when implemented correctly, an onsite response strategy is the best option from a security and cost perspective.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Fungal Community, Not Substrate Quality, Drives Soil Microbial Function in Northeastern U.S. Temperate Forests

Mycorrhizal fungi can affect soil organic matter cycling through several mechanisms including priming, nutrient competition, and direct enzyme production. Differences in nutrient foraging strategies between ectomycorrhizal (EcM) and arbuscular mycorrhizal (AM) fungi produce divergent belowground dynamics: where EcM can take up organic nitrogen and directly break down soil organic matter (SOM) by producing enzymes, AM fungi are limited to scavenging mineral N. EcM-associated tree species also have leaf litter with relatively higher ratios of carbon to nitrogen (C:N), and belowground saprotrophic communities more dominated by fungi. Consequently, free-living microbes in EcM-dominated soils should experience nitrogen limitation, with subsequent increases in enzyme production and decreased carbon use efficiency (CUE). However, the relative importance of the effects of substrate quality and fungal community composition on enzyme production and CUE are unclear. To assess this distinction, we sampled the organic horizon and 10 cm of the mineral horizon in northern temperate forest soils along a gradient of EcM dominance. We characterized fungal community composition by measuring EcM relative abundances from extracted fungal DNA and the fungal to bacterial (F:B) ratios from phospholipid fatty acid (PLFA) analysis. We assessed soil substrate quality as the soil C:N ratio. Soil microbial functions were measured as potential activities of five hydrolytic and two oxidative enzymes, and microbial CUE. We found that the fungal community, represented by either the F:B ratio, EcM relative abundance, or both, affected CUE and six measured enzyme activities, while the C:N ratio affected only oxidative and chitin-targeting extracellular enzyme activities. Our results highlight the use of EcM relative dominance as a predictor of soil microbial community composition and function independent of substrate quality.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Holographic reflective Bragg gratings of gain and their role in the operation of high-power pulsed lasers

The influence of reflective Bragg gratings that are formed in an active medium due to the spatial periodic modulation of gain on the spectral properties of lasing are investigated by an example of a pulsed Nd : YAG laser. It is shown that, at a certain choice of the formation regime of reflective gratings of gain, the lasing bandwidth can be reduced by a factor of 3 – 5, with a corresponding increase in the pulse energy and significant improvement of spatial beam quality. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗