Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic structural analysis”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Elucidating the structure of the magnesium aluminum chloride complex electrolyte for magnesium-ion batteries

Non-aqueous Mg-ion batteries offer a promising way to overcome safety, costs, and energy density limitations of state-of-the-art Li-ion battery technology. Here, we present a rigorous analysis of the magnesium aluminum chloride complex (MACC) in tetrahydrofuran (THF), one of the few electrolytes that can reversibly plate and strip Mg. We use ab initio calculations and classical molecular dynamics simulations to interrogate the MACC electrolyte composition with the goal of addressing two urgent questions that have puzzled battery researchers: (i) the functional species of the electrolyte, and (ii) the complex equilibria regulating the MACC speciation after prolonged electrochemical cycling, a process termed as conditioning, and after prolonged inactivity, a process called aging. A general computational strategy to untangle the complex structure of electrolytes, ionic liquids and other liquid media is presented. The analysis of formation energies and grand-potential phase diagrams of Mg–Al–Cl–THF suggests that the MACC electrolyte bears a simple chemical structure with few simple constituents, namely the electro-active species MgCl + and AlCl 4 - in equilibrium with MgCl 2 and AlCl 3 . Knowledge of the stable species of the MACC electrolyte allows us to determine the most important equilibria occurring during electrochemical cycling. We observe that Al deposition is always preferred to Mg deposition, explaining why freshly synthesized MACC cannot operate and needs to undergo preparatory conditioning. Similarly, we suggest that aluminum displacement and depletion from the solution upon electrolyte resting (along with continuous MgCl 2 regeneration) represents one of the causes of electrolyte aging. Finally, we compute the NMR shifts from shielding tensors of selected molecules and ions providing fingerprints to guide future experimental investigations.

25 ENERGY STORAGE↗
DOE OSTI2016

Forward-backward multiplicity fluctuation and longitudinal harmonics in high-energy nuclear collisions

In this paper, an analysis method is proposed to study the forward-backward (FB) multiplicity fluctuation in high-energy nuclear collisions, built on the earlier work of Bzdak and Teaney [Phys. Rev. C 87, 024906 (2013)]. The method allows the decomposition of the centrality dependence of average multiplicity from the dynamical event-by-event (EbyE) fluctuation of multiplicity in pseudorapidity. Application of the method to AMPT (A Multi-Phase Transport model) and HIJING (Heavy Ion Jet INteraction Generator) models shows that the long-range component of the FB correlation is captured by a few longitudinal harmonics, with the first component driven by the asymmetry in the number of participating nucleons in the two colliding nuclei. The higher-order longitudinal harmonics are found to be strongly damped in AMPT compared to HIJING, due to weaker short-range correlations as well as the final-state effects present in the AMPT model. Two-particle pseudorapidity correlation reveals interesting charge-dependent short-range structures that are absent in HIJING model. Lastly, the proposed method opens an avenue to elucidate the particle production mechanism and early time dynamics in heavy-ion collisions. Future analysis directions and prospects of using the pseudorapidity correlation function to understand the centrality bias in p + p, p + A, and A + A collisions are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Experimental characterization and atomistic modeling of interfacial void formation and detachment in short pulse laser processing of metal surfaces covered by solid transparent overlayers

The short pulse laser interaction with metal surfaces covered by solid transparent overlayers is investigated in experiments and atomistic simulations, with a particular aim of revealing the mechanisms responsible for structural modification of the metal–overlayer interfacial regions. Experimental characterization of Al–silica targets modified by single-pulse laser irradiation with the pulse duration of 10 ps reveals the transitions from the generation of extended interfacial voids with internal nanoscale surface roughness to the partial detachment of the overlayer from the metal substrate, and to the cracking/chipping or complete removal of the overlayer as the laser fluence increases. The mechanisms responsible for the appearance, growth, and percolation of the interfacial voids leading to the detachment of the overlayer from the metal substrate are investigated in a large-scale atomistic simulation. The results of the simulation demonstrate that the processes of nucleation and growth of the interfacial voids are driven by the dynamic relaxation of laser-induced stresses proceeding simultaneously with rapid phase transformations and temperature variation in the interfacial region. The growth and coalescence of the interfacial voids results in the formation of liquid bridges connecting the overlayer and the metal substrate, whereas solidification of the transient liquid structures produced by the breakup of the liquid bridges may be responsible for the formation of the nanoscale roughness of the interfacial voids observed in experiments. Computational analysis of the effect of preexisting interfacial voids reveals a complex dynamic picture of the initial expansion and subsequent compaction of the surface region of the metal substrate and suggests a possible scenario for the formation of voids below the metal–overlayer interface.

Materials Science↗
DOE OSTI2021

How Do Ionic Liquids “Fold” Ionenes? Computational and Experimental Analysis of Imidazolium Polymers Based on Ether and Alkyl Chain Variations Dissolved in an Ionic Liquid

Molecular dynamics (MD) simulations and complementary experiments are used here to understand the inter- and intramolecular structure and conformational properties of ionenes (i.e., cationic polymers formed from condensation reactions) dissolved in an ionic liquid (IL), 1-ethyl-3-methylimidazolium bistriflimide [C 2 mim + ][Tf 2 N - ]. The simulated structural properties are benchmarked against experimental analyses of these same polymers, mainly using dynamic light scattering (DLS) experiments. Four different imidazolium ionenes are considered, corresponding to variations in the chemistry and structure of the repeat units including: (1) poly (decylimidazolium) (PD 10 ); (2) poly (tetraethyleneglycolimidazolium) (PE 10 ); (3) alternating copolymer (P(ED) 5 ); and (4) a block copolymer (PE 5 D 5 ). Detailed computational analyses of the polymer structure and conformational properties were performed, including the radius of gyration, end-to-end distance, torsional distributions, and site-site and spatial distribution functions. Overall, there is a competition between intramolecular associations between the imidazolium groups and the ether sites in the polymers versus the intermolecular interactions of the polymer imidazolium groups with the surrounding anion molecules. The polymer with only ether linkages (PE 10 ) results in the most intramolecular interactions, leading to significant coiling behavior and chain contraction. These strong interactions reduce chain flexibility, but they also result in much more linear chain configurations and alignment of the imidazolium groups. The polymer structure is strongly affected as the concentration of alkyl groups is increased, and we find that the polymer architecture (alternating versus block copolymer) also has an important influence. Furthermore, these detailed observations are important for understanding the connection between the molecular design of imidazolium polymers and their emergent structural properties, which could lead to unique opportunities for creating polymer composites.

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

Designing the Protocols for Programmable Ammonia Catalysis

Programmable catalysis can provide a more energy-efficient and cost-effective route to enhancing commercial ammonia production, a key process in the advancement of renewable energy technologies and the manufacture of fertilizers and basic chemicals. This work explores the computational discovery of optimal forcing protocols to drive such dynamic catalysis models. By employing matrix-free time-stepper methods, coupled with an optimization approach, that integrates Bayesian optimization with a Bayesian continuation strategy to efficiently discover the periodic steady states of such periodically forced systems, we enable the discovery of complex optimal catalyst strain waveforms, while ensuring robust solver convergence. We demonstrate the flexibility of our approach to discover optimized forcing protocols under varying physical constraints on strain modulation or other catalyst operating parameters. We show that these can have a temporal structure more complex than simple step functions. In order to detect undesirable catalytic loops that may correlate with overall reduced performance, we perform a study using graph-theoretical analysis to investigate the dynamics of catalytic kinetic networks formed.

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

Molecular dynamics investigation of threshold displacement energies in CaF 2

Understanding the propagation of radiation damage in a material is paramount to predicting the material damage effects. To date, no current literature has investigated the Threshold Displacement Energy (TDE) of Ca and F atoms in CaF 2 through molecular dynamics and simulated statistical analysis. A set of interatomic potentials between Ca-Ca, F-F, and F-Ca were splined, fully characterizing a pure CaF 2 simulation cell, by using published Born-Mayer-Huggins, standard ZBL, and Coulomb potentials, with a resulting structure within 1% of standard density and published lattice constants. Using this simulation cell, molecular dynamics simulations were performed with LAMMPS using a simulation that randomly generated 500 Ca and F PKA directions for each incremental set of energies, and a simulation in each of the [1 0 0], [1 1 0], and [1 1 1] directions with 500 trials for each incremental energy. MD simulations of radiation damage in CaF 2 are carried out using F and Ca PKAs, with energies ranging from 2 to 200 eV. Probabilistic determinations of the TDE and Threshold Vacancy Energy (TVE) of Ca and F atoms in CaF 2 were performed, as well as examining vacancy, interstitial, and antisite production rates over the range of PKA energies. Many more F atoms were displaced from both PKA species, and though F recombination appears more probable than Ca recombination, F vacancy numbers are higher. In conclusion, the higher number of F vacancies than Ca vacancies suggests F Frenkel pairs dominate CaF 2 damage.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Characterization of Complex Interactions at the Gas–Solid Interface with in Situ Spectroscopy: The Case of Nitrogen-Functionalized Carbon

Interactions at the gas-solid interface drive physicochemical processes in many energy and environmental applications; yet, the challenges associated with characterization and development of these dynamic interactions in complex systems limit progress in developing effective materials. Thus, structure-property-performance correlations greatly depend on the development of advanced techniques and analysis methods for the investigation of gas-solid interactions. In this work, adsorption behavior of O 2 and humidified O 2 on nitrogen-functionalized carbon (N-C) materials was investigated to provide a better understanding of the role of nitrogen species in the oxygen reduction reaction (ORR). N-C materials were produced by solvothermal synthesis and N-ion implantation, resulting in a set of materials with varied nitrogen amount and speciation in carbon matrices with different morphologies. Adsorption behavior of the N-C samples was characterized by in situ diffuse reflectance infrared Fourier-transform spectroscopy (DRIFTS) and ambient pressure X-ray photoelectron spectroscopy (AP-XPS) experiments. A novel analysis approach for the interpretation of AP-XPS data was developed, allowing both the determination of overall adsorption behavior of each N-C material and identification of which nitrogen species were responsible for adsorption. The complementary information provided by in situ DRIFTS and AP-XPS suggests that O 2 adsorption primarily takes place on either electron-rich nitrogen species like pyridine, hydrogenated nitrogen species, or graphitic nitrogen. Adsorption of O 2 and H 2 O occurs competitively on solvothermally prepared N-Cs, whereas adsorption of H 2 O and O 2 occurs at different sites on N-ion implanted N-Cs, highlighting the importance of tuning the composition of N-C materials to promote the most efficient ORR pathway.

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

Time-resolved scanning electron microscopy with polarization analysis

We demonstrate the feasibility of investigating periodically driven magnetization dynamics in a scanning electron microscope with polarization analysis based on spin-polarized low-energy electron diffraction. With the present setup, analyzing the time structure of the scattering events, we obtain a temporal resolution of 700 ps, which is demonstrated by means of imaging the field-driven 100 MHz gyration of the vortex in a soft-magnetic FeCoSiB square. Owing to the efficient intrinsic timing scheme, high-quality movies, giving two components of the magnetization simultaneously, can be recorded on the time scale of hours.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Metatranscriptomic reconstruction reveals RNA viruses with the potential to shape carbon cycling in soil

Viruses impact nearly all organisms on Earth, with ripples of influence in agriculture, health, and biogeochemical processes. However, very little is known about RNA viruses in an environmental context, and even less is known about their diversity and ecology in soil, 1 of the most complex microbial systems. Here, we assembled 48 individual metatranscriptomes from 4 habitats within a planted soil sampled over a 22-d time series: Rhizosphere alone, detritosphere alone, rhizosphere with added root detritus, and unamended soil (4 time points and 3 biological replicates). We resolved the RNA viral community, uncovering a high diversity of viral sequences. We also investigated possible host organisms by analyzing metatranscriptome marker genes. Based on viral phylogeny, much of the diversity was Narnaviridae that may parasitize fungi or Leviviridae , which may infect Proteobacteria. Both host and viral communities appear to be highly dynamic, and rapidly diverged depending on experimental conditions. The viral and host communities were structured based on the presence of root litter. Clear temporal dynamics by Leviviridae and their hosts indicated that viruses were replicating. With this time-resolved analysis, we show that RNA viruses are diverse, abundant, and active in soil. When viral infection causes host cell death, it may mobilize cell carbon in a process that may represent an overlooked component of soil carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Inverse Methods - Users Manual 5.6

The inverse methods team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verfication/validation analysis of weapons systems. This document provides a user's guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.

42 ENGINEERING↗
DOE OSTI2022

Inverse Methods - Users Manual 5.8

The inverse methods team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verification/validation analysis of weapons systems. This document provides a user’s guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

SISGR: The regulation of carbon fixation in plant and green algae: Rubisco activase and the origin of heat inactivation of CO 2 assimilation

Rubisco activase (Rca) is a critical AAA+ ATPase protein complex that remodels and promotes the Rubisco enzyme, a key player in photosynthetic performance and carbon fixation. The assembly and function of the Rca protein complex are regulated by a range of factors, including subunit concentration, nucleotide-binding states, thermal conditions, metal-ion coordination, and post-translational modifications, such as phosphorylation. Despite its importance in photosynthesis, the detailed molecular mechanisms underlying the regulation of plant Rca and how it activates Rubisco remain elusive. This project aims to bridge this knowledge gap by integrating sophisticated enzymology tools with single-molecule methods and high-resolution electron microscopy to elucidate the structure and function of plant Rca. Through these multiple approaches, we have systematically investigated how the activity of plant Rca is impacted by various factors, such as phosphorylation and metal-ion coordination. The Rca complex assembly/disassembly dynamics were captured using anti-Brownian electrokinetic (ABEL) trap-based measurements, providing unprecedented insight into its structural flexibility and diverse assembly states. Furthermore, the structural analysis of Rca through electron crystallography and single-particle cryogenic electron microscopy (cryo-EM) reveals novel assembly states of the spinach Rca, providing insight into the mechanistic action for Rubisco remodeling. By combining cutting-edge tools and approaches, this work uncovers critical aspects of Rca’s regulation and assembly, paving the way for a deeper understanding of its role in photosynthetic efficiency and the potential for enhancing carbon fixation in crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

FuSED Users Manual, 5.24

The Fusion of Simulation, Experiment, and Data (FuSED) team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verification/validation analysis of systems. This document provides a user’s guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Enhancing Multi-Step Reservoir Inflow Forecasting: A Time-Variant Encoder–Decoder Approach

Accurate reservoir inflow forecasting is vital for effective water resource management. Reliable forecasts enable operators to optimize storage and release strategies to meet competing sectoral demands—such as water supply, irrigation, and hydropower scheduling—while also mitigating flood and drought risks. To address this need, in this study, we propose a novel time-variant encoder–decoder (ED) model designed specifically to improve multi-step reservoir inflow forecasting, enabling accurate predictions of reservoir inflows up to seven days ahead. Unlike conventional ED-LSTM and recursive ED-LSTM models, which use fixed encoder parameters or recursively propagate predictions, our model incorporates an adaptive encoder structure that dynamically adjusts to evolving conditions at each forecast horizon. Additionally, we introduce the Expected Baseline Integrated Gradients (EB-IGs) method for variable importance analysis, enhancing interpretability of inflow by incorporating multiple baselines to capture a broader range of hydrometeorological conditions. The proposed methods are demonstrated at several diverse reservoirs across the United States. Our results show that they outperform traditional methods, particularly at longer lead times, while also offering insights into the key drivers of inflow forecasting. These advancements contribute to enhanced reservoir management through improved forecasting accuracy and practical decision-making insights under complex hydroclimatic conditions.

58 GEOSCIENCES↗
DOE OSTI2022

Advanced modeling and simulation of research reactors using dynamic mode decomposition

Full text of publication follows. Due to the ever-increasing safety requirements, the current trend of nuclear reactor analysis is shifting towards high-fidelity multi-physics models, which have a very high computational cost and modelling complexity. As the cost of even a single model run makes it impossible to analyse the behaviour and performance of these models on large-scale commercial plants, it has become even more significant to provide suitable benchmarks to validate and test them extensively. In this sense, research reactors offer a promising solution for the initial validation of high-fidelity models, as they are significantly smaller than commercial reactors and their characteristics are well known. In particular, the reactors of the TRIGA family have been used to assess and validate models and methods for Generation-IV designs, as they have some similar features (such as the dominance of natural convection as cooling mechanism and the difficulties in performing sub-channel analysis using standard codes). Still, the computational requirements of high-fidelity models make them unsuitable for real-time analysis, even following their assessment on research reactors. In this sense, Model Order Reduction (MOR) techniques give an additional strategy to reduce the computational cost of high-fidelity models (whilst preserving sufficient accuracy). In particular, this work focuses on Dynamic Mode Decomposition (DMD), a non-intrusive MOR technique that aims at representing models with explicit temporal dynamics by extracting the time-varying characteristics and the governing structures based only on a set of available data, thus without needing any underlying knowledge of the governing equations. In addition, DMD also computes a low-dimensional surrogate of the dynamic matrix of the system, making it suited for stability analysis and real-time evaluations. This work focuses on the application and validation of the DMD method on the Computational Fluid-Dynamics (CFD) model TRIGA Mark II reactor, also discussing in detail the potentiality of this algorithm as an advanced modelling tool for nuclear reactor analysis. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A thin-walled composite beam model for light-weighted structures interacting with fluids

A thin-walled beam model is proposed for structures of variable cross-section, which can be either open or closed and includes multicellular cross-sections with either isotropic or orthotropic materials. The proposed model does not require any priori definition of cross-sectional warping which instead results from the solution of the problem. To achieve that a special deformation pattern is superimposed on the bending deformation described by Euler–Bernoulli beam theory. All sectional properties are automatically incorporated in the analysis as a result of the usual variational formulation of the system of equations. The proposed model is specifically designed to simulate the dynamics of wind/hydrokinetic turbine blade with low computational cost, especially in fluid–structure interaction (FSI) simulation. A number of test cases have been carried out to validate the proposed structural model which show good agreement between the results obtained her e and the solutions available in literature. Finally, FSI simulation of a hydrokinetic blade under field condition is carried out to illustrate the capability of the current thin-walled beam model in practice.

42 ENGINEERING↗
DOE OSTI2023

Structure of a 28.5 kDa duplex-embedded G-quadruplex system resolved to 7.4 Å resolution with cryo-EM

Genomic regions with high guanine content can fold into non-B form DNA four-stranded structures known as G-quadruplexes (G4s). Extensive in vivo investigations have revealed that promoter G4s are transcriptional regulators. Little structural information exists for these G4s embedded within duplexes, their presumed genomic environment. Here, we report the 7.4 Å resolution structure and dynamics of a 28.5 kDa duplex-G4-duplex (DGD) model system using cryo-EM, molecular dynamics, and small-angle X-ray scattering (SAXS) studies. The DGD cryo-EM refined model features a 53° bend induced by a stacked duplex-G4 interaction at the 5’ G-tetrad interface with a persistently unstacked 3’ duplex. The surrogate complement poly dT loop preferably stacks onto the 3’ G-tetrad interface resulting in occlusion of both 5’ and 3’ tetrad interfaces. Structural analysis shows that the DGD model is quantifiably more druggable than the monomeric G4 structure alone and represents a new structural drug target. Our results illustrate how the integration of cryo-EM, MD, and SAXS can reveal complementary detailed static and dynamic structural information on DNA G4 systems.

59 BASIC BIOLOGICAL SCIENCES↗