Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structure modelling”

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 127 records · Page 7

Resolving Atomic-Scale Structure and Chemical Coordination in High-Entropy Alloy Electrocatalysts for Structure–Function Relationship Elucidation

The recent breakthrough in confining five or more atomic species in nanocatalysts, referred to as high-entropy alloy nanocatalysts (HEAs), has revealed the possibilities of multielemental interactions that can surpass the limitations of binary and ternary electrocatalysts. The wide range of potential surface configurations in HEAs, however, presents a significant challenge in resolving active structural motifs, preventing the establishment of structure-function relationships for rational catalyst design and optimization. Here, we present a methodology for creating sub-5 nm HEAs using an aqueous-based peptide-directed route. Using a combination of pair distribution function and X-ray absorption spectroscopy, HEA structure models are constructed from reverse Monte Carlo modeling of experimental data sets and showcase a clear peptide-induced influence on atomic-structure and chemical miscibility. Coordination analysis of our structure models facilitated the construction of structure-function correlations applied to electrochemical methanol oxidation reactions, revealing the complex interplay between multiple metals that leads to improved catalytic properties. Our results showcase a viable strategy for elucidating structure-function relationships in HEAs, prospectively providing a pathway for future materials design.

36 MATERIALS SCIENCE↗

Oxygen Defect Configurations in Single-Phase UO 2.15

Hyperstoichiometric UO 2.15 was characterized by neutron total scattering at high temperature in the single-phase UO 2+x region of the U/O phase diagram. The diffraction data confirmed a single-phase fluorite structure at high temperature. Analysis of the short-range data showed that the same structural model does not fit the pair distribution functions well. Instead, structural models containing specific configurations of oxygen defect clusters best represent the local atomic arrangement. In conclusion, prevalent defect clusters previously proposed were fit to the experimental data, and moderately distorted oxygen cuboctahedra hypothesized by recent molecular dynamics simulations fit the data most accurately.

neutron pair-distribution function analysis↗

High throughput, accurate gene annotation through AI and HPC-enabled structural analysis

With the advances in next generation sequencing technologies, the number of sequenced genomes is growing exponentially, resulting in a technology bottleneck for the translation of sequence information into usable hypotheses about the function of each gene. We have proposed leveraging our leadership high-performance computing (HPC) resources to help break this annotation bottleneck. Here we design an HPC-based framework to infer gene function from gene sequence by incorporating information about protein structure and interactions predicted by deep learning approaches. Accurate functional prediction and gene annotation using computational methods will facilitate breakthroughs in the genomic sciences essential to understanding and harnessing life processes in bacteria, fungi and plants. The development and applications of the state-of-the-art deep neural networks to protein structural modeling, interaction prediction, sequence comparison, and quality assessment of protein structural models will be made possible by leadership computational resources. These HPC-enabled bioinformatics and molecular modeling tools will provide powerful insights into molecular functions of genes.

59 BASIC BIOLOGICAL SCIENCES↗

A Review on Lignin Liquefaction: Advanced Characterization of Structure and Microkinetic Modeling

Lignin liquefaction microkinetics is a move toward a more first principles (i.e., ab initio) based understanding at the molecular level in reaction engineering. While the microkinetic modelling of reactions to obtain kinetic rate parameters of chemical reactions have been widely used in the field of gas phase combustion and heterogeneous catalysis, this approach has not been as thoroughly developed in the area of biomass thermochemical reactions (e.g., lignin pyrolysis, hydrothermal liquefaction). The difficulties in establishing the structure of complex heterogeneous materials, like lignin, is perhaps the main challenge in developing rational microkinetic descriptions of biomass thermochemical reactions. In this manuscript we review the current state of the art and the challenges to develop microkinetic models for lignin liquefaction technologies (e.g., pyrolysis, hydrothermal liquefaction, solvolysis). A general strategy for the development of microkinetic models for lignin liquefaction technologies is discussed. Here, the first hurdle is to obtain sufficiently rich experimental data of lignin underlying polymeric structure and methodologies to use this data to build realistic lignin structural representations. Some analytical techniques for lignin structural characterization and their associated data, as well as a correlation for calculating the degree of macromolecular lignin branching are discussed. The presence of small lignin oligomeric structures and the role of these structures in lignin pyrolysis is also addressed. The ways in which elementary deconstruction and repolymerization reactions occur within this structure to form a liquid intermediate, and how these deconstruction products continue to interact with each other until they are removed from the liquid intermediate is thoroughly discussed. Further, experimental work with model compounds and the effect of reaction parameters (e.g., temperature, pressure, vapor residence time) are reviewed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics-based halo model for large scale structure

Accurate modeling of the one-to-two halo transition has long been difficult to achieve. Here, we demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to (1) particles orbiting the halo and (2) particles infalling into the halo for the first time. Motivated by this decomposition, García et al. [Mon. Not. R. Astron. Soc. 521, 2464 (2023)] advocated for defining halos as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo-mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterized by a single physical scale 𝑟 h , which we refer to as the halo radius. The two-halo term is described using a simple yet accurate empirical model based on the Zel’dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo-mass correlation function is accurate at the ≈ 2% level for 𝑟∈ [0.1, 50] ℎ −1 Mpc. The Fourier transform of our best-fit model describes the halo-mass power spectrum with comparable accuracy for 𝑘 ∈ [0.06, 6.0] ℎ Mpc −1 .

79 ASTRONOMY AND ASTROPHYSICS↗

Shaker-structure interaction modeling and analysis for nonlinear force appropriation testing

Nonlinear force appropriation is an extension of its linear counterpart where sinusoidal excitation is applied to a structure with a modal shaker and phase quadrature is achieved between the excitation and response. While a standard practice in modal testing, modal shaker excitation has the potential to alter the dynamics of the structure under test. Previous studies have been conducted to address several concerns, but this work specifically focuses on a shaker-structure interaction phenomenon which arises during the force appropriation testing of a nonlinear structure. Under pure-tone sinusoidal forcing, a nonlinear structure may respond not only at the fundamental harmonic but also potentially at sub- or superharmonics, or it can even produce aperiodic and chaotic motion in certain cases. Shaker-structure interaction occurs when the response physically pushes back against the shaker attachment, producing non-fundamental harmonic content in the force measured by the load cell, even for pure tone voltage input to the shaker. This work develops a model to replicate these physics and investigates their influence on the response of a nonlinear normal mode of the structure. Experimental evidence is first provided that demonstrates the generation of harmonic content in the measured load cell force during a force appropriation test. This interaction is replicated by developing an electromechanical model of a modal shaker attached to a nonlinear, three-mass dynamical system. Several simulated experiments are conducted both with and without the shaker model in order to identify which effects are specifically due to the presence of the shaker. Finally, the results of these simulations are then compared to the undamped nonlinear normal modes of the structure under test to evaluate the influence of shaker-structure interaction on the identified system’s dynamics.

42 ENGINEERING↗

Manifestations of metastable criticality in the long-range structure of model water glasses

Abstract Much attention has been devoted to water’s metastable phase behavior, including polyamorphism (multiple amorphous solid phases), and the hypothesized liquid-liquid transition and associated critical point. However, the possible relationship between these phenomena remains incompletely understood. Using molecular dynamics simulations of the realistic TIP4P/2005 model, we found a striking signature of the liquid-liquid critical point in the structure of water glasses, manifested as a pronounced increase in long-range density fluctuations at pressures proximate to the critical pressure. By contrast, these signatures were absent in glasses of two model systems that lack a critical point. We also characterized the departure from equilibrium upon vitrification via the non-equilibrium index; water-like systems exhibited a strong pressure dependence in this metric, whereas simple liquids did not. These results reflect a surprising relationship between the metastable equilibrium phenomenon of liquid-liquid criticality and the non-equilibrium structure of glassy water, with implications for our understanding of water phase behavior and glass physics. Our calculations suggest a possible experimental route to probing the existence of the liquid-liquid transition in water and other fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

RNA target highlights in CASP15 : Evaluation of predicted models by structure providers

Abstract The first RNA category of the Critical Assessment of Techniques for Structure Prediction competition was only made possible because of the scientists who provided experimental structures to challenge the predictors. In this article, these scientists offer a unique and valuable analysis of both the successes and areas for improvement in the predicted models. All 10 RNA‐only targets yielded predictions topologically similar to experimentally determined structures. For one target, experimentalists were able to phase their x‐ray diffraction data by molecular replacement, showing a potential application of structure predictions for RNA structural biologists. Recommended areas for improvement include: enhancing the accuracy in local interaction predictions and increased consideration of the experimental conditions such as multimerization, structure determination method, and time along folding pathways. The prediction of RNA–protein complexes remains the most significant challenge. Finally, given the intrinsic flexibility of many RNAs, we propose the consideration of ensemble models.

59 BASIC BIOLOGICAL SCIENCES↗

Partitioning uncertainty in projections of Arctic sea ice

Abstract Improved knowledge of the contributing sources of uncertainty in projections of Arctic sea ice over the 21st century is essential for evaluating impacts of a changing Arctic environment. Here, we consider the role of internal variability, model structure and emissions scenario in projections of Arctic sea-ice area (SIA) by using six single model initial-condition large ensembles and a suite of models participating in Phase 5 of the Coupled Model Intercomparison Project. For projections of September Arctic SIA change, internal variability accounts for as much as 40%–60% of the total uncertainty in the next decade, while emissions scenario dominates uncertainty toward the end of the century. Model structure accounts for 60%–70% of the total uncertainty by mid-century and declines to 30% at the end of the 21st century in the summer months. For projections of wintertime Arctic SIA change, internal variability contributes as much as 50%–60% of the total uncertainty in the next decade and impacts total uncertainty at longer lead times when compared to the summertime. In winter, there exists a considerable scenario dependence of model uncertainty with relatively larger model uncertainty under strong forcing compared to weak forcing. At regional scales, the contribution of internal variability can vary widely and strongly depends on the calendar month and region. For wintertime SIA change in the Greenland-Iceland-Norwegian and Barents Seas, internal variability contributes 60%–70% to the total uncertainty over the coming decades and remains important much longer than in other regions. We further find that the relative contribution of internal variability to total uncertainty is state-dependent and increases as sea ice volume declines. These results demonstrate that internal variability is a significant source of uncertainty in projections of Arctic sea ice.

Bonan, David B. (ORCID:0000000338676009)↗

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗

Improving Mathematical Exposition of an Industrial-Scale Linear Program

Industrial-scale models require considerable setup time; hence, once built, they are used in myriad ways to consider closely related cases. In practice, the code for these models frequently evolves without appropriate notational choices, largely as a result of the lengthy development time of, and the number of individuals contributing to, their formulation. This leads to inefficiencies and obfuscates model structures that might be leveraged to expedite solutions. In this paper, we advocate for an emerging literature on model formulation “best practices” and present the reformulation of a widely used industrial-scale linear program. The efficient mathematical expression of this linear program, used to plan capacity expansion in the energy sector, allows for greater transparency of model structures and enhanced ability to identify computational performance improvements, as well as a lucid interpretation of its solutions. This type of formulation is employed in several mathematical programming courses at our university as an example of the advantages of best practices; the model more broadly is used widely to inform policy in the U.S. energy sector.

97 MATHEMATICS AND COMPUTING↗

A measurement of the equation of state of carbon envelopes of white dwarfs

White dwarfs represent the final state of evolution for most stars. Certain classes of white dwarfs pulsate, leading to observable brightness variations, and analysis of these variations with theoretical stellar models probes their internal structure. Modelling of these pulsating stars provides stringent tests of white dwarf models and a detailed picture of the outcome of the late stages of stellar evolution. However, the high-energy-density states that exist in white dwarfs are extremely difficult to reach and to measure in the laboratory, so theoretical predictions are largely untested at these conditions. In this paper we report measurements of the relationship between pressure and density along the principal shock Hugoniot (equations describing the state of the sample material before and after the passage of the shock derived from conservation laws) of hydrocarbon to within five per cent. The observed maximum compressibility is consistent with theoretical models that include detailed electronic structure. This is relevant for the equation of state of matter at pressures ranging from 100 million to 450 million atmospheres, where the understanding of white dwarf physics is sensitive to the equation of state and where models differ considerably. The measurements test these equation-of-state relations that are used in the modelling of white dwarfs and inertial confinement fusion experiments, and we predict an increase in compressibility due to ionization of the inner-core orbitals of carbon. We also find that a detailed treatment of the electronic structure and the electron degeneracy pressure is required to capture the measured shape of the pressure–density evolution for hydrocarbon before peak compression. Our results illuminate the equation of state of the white dwarf envelope (the region surrounding the stellar core that contains partially ionized and partially degenerate non-ideal plasmas), which is a weak link in the constitutive physics informing the structure and evolution of white dwarf stars.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantification of hydraulic trait control on plant hydrodynamics and risk of hydraulic failure within a demographic structured vegetation model in a tropical forest (FATES–HYDRO V1.0)

Abstract. Vegetation plays a key role in the global carbon cycle and thus is an important component within Earth system models (ESMs) that project future climate. Many ESMs are adopting methods to resolve plant size and ecosystem disturbance history, using vegetation demographic models. These models make it feasible to conduct more realistic simulation of processes that control vegetation dynamics. Meanwhile, increasing understanding of the processes governing plant water use, and ecosystem responses to drought in particular, has led to the adoption of dynamic plant water transport (i.e., hydrodynamic) schemes within ESMs. However, the extent to which variations in plant hydraulic traits affect both plant water stress and the risk of mortality in trait-diverse tropical forests is understudied. In this study, we report on a sensitivity analysis of an existing hydrodynamic scheme (HYDRO) model that is updated and incorporated into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (FATES–HYDRO V1.0). The size- and canopy-structured representation within FATES is able to simulate how plant size and hydraulic traits affect vegetation dynamics and carbon–water fluxes. To better understand this new model system, and its functionality in tropical forest systems in particular, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We assembled 942 observations of plant hydraulic traits on 306 tropical plant species for stomata, leaves, stems, and roots and determined the best-fit statistical distribution for each trait, which was used in model parameter sampling to assess the parametric sensitivity. We showed that, for simulated leaf water potential and loss of hydraulic conductivity across different plant organs, the four most important traits were associated with xylem conduit taper (buffers increasing hydraulic resistance with tree height), stomatal sensitivity to leaf water potential, maximum stem hydraulic conductivity, and the partitioning of total hydraulic resistance above vs. belowground. Our analysis of individual ensemble members revealed that trees at a high risk of hydraulic failure and potential tree mortality generally have a lower conduit taper, lower maximum xylem conductivity, lower stomatal sensitivity to leaf water potential, and lower resistance to xylem embolism for stem and transporting roots. We expect that our results will provide guidance on future modeling studies using plant hydrodynamic models to predict the forest responses to droughts and future field campaigns that aim to better parameterize plant hydrodynamic models.

54 ENVIRONMENTAL SCIENCES↗

Liquid state theory of the structure of model polymerized ionic liquids

We employ polymer integral equation theory to study a simplified model of semiflexible polymerized ionic liquids (PolyILs) that interact via hard core repulsions and short range screened Coulomb interactions. The multi-scale structure in real and Fourier space of PolyILs (ions chosen to mimic Li, Na, K, Br, PF 6 , and TFSI) are determined as a function of melt density, Coulomb interaction strength, and ion size. Comparisons with a homopolymer melt, a neutral polymer–solvent-like athermal mixture, and an atomic ionic liquid are carried out to elucidate the distinct manner that ions mediate changes of polymer packing, the role of excluded volume effects, and the influence of chain connectivity, respectively. The effect of Coulomb strength depends in a rich manner on ion size and density, reflecting the interplay of steric packing, ion adsorption, and charge layering. Ion-mediated bridging of monomers is found, which intensifies for larger ions. Intermediate range charge layering correlations are characterized by a many-body screening length that grows with PolyIL density, cooling, and Coulomb strength, in disagreement with Debye–Hückel theory, but in accord with experiments. Qualitative differences in the collective structure, including an ion-size-dependent bifurcation of the polymer structure factor peak and pair correlation function, are predicted. The monomer cage order parameter increases significantly, but its collective ion counterpart decreases, as ions become smaller. Such behaviors allow one to categorize PolyILs into two broad classes of small and large ions. Furthermore, dynamical implications of the predicted structural results are qualitatively discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Work Hardening Model of Structure in Hall-Petch Strengthening

The microstructural length scale of metals changes by orders of magnitude under extreme processing conditions producing a concurrent wide range of mechanical strength and plasticity behaviors. A unified stress-strain σε σ(ε) model is formulated that’s based on superposing the components of asymptotic-curvilinear work hardening Θσ Θ(σ) to qualify and quantify these mechanical behaviors. This approach accounts for the rapid strengthening of metals beyond the initial yield point, through extended steady-state deformation, to the structural instability. The relationship between the softening coefficients cbi c bi of the work hardening formulation Θσ Θ(σ) and strength are found to reveal the microstructural scale in the material. Specifically, the rapid decrease in the slope of the Θσ Θ(σ) curve provides a measure for microstructural size consistent with a functional Hall-Petch relationship of strength. A successful application is shown for the tensile behavior of pure aluminum microstructures that result from extreme plastic deformation by equal-channel angle pressing.

Hall-Petch strengthening↗

KGML-ag: a modeling framework of knowledge-guided machine learning to simulate agroecosystems: a case study of estimating N<sub>2</sub>O emission using data from mesocosm experiments

Abstract. Agricultural nitrous oxide (N2O) emission accounts for a non-trivial fraction of global greenhouse gas (GHG) budget. To date, estimating N2O fluxes from cropland remains a challenging task because the related microbial processes (e.g., nitrification and denitrification) are controlled by complex interactions among climate, soil, plant and human activities. Existing approaches such as process-based (PB) models have well-known limitations due to insufficient representations of the processes or uncertainties of model parameters, and due to leverage recent advances in machine learning (ML) a new method is needed to unlock the “black box” to overcome its limitations such as low interpretability, out-of-sample failure and massive data demand. In this study, we developed a first-of-its-kind knowledge-guided machine learning model for agroecosystems (KGML-ag) by incorporating biogeophysical and chemical domain knowledge from an advanced PB model, ecosys, and tested it by comparing simulating daily N2O fluxes with real observed data from mesocosm experiments. The gated recurrent unit (GRU) was used as the basis to build the model structure. To optimize the model performance, we have investigated a range of ideas, including (1) using initial values of intermediate variables (IMVs) instead of time series as model input to reduce data demand; (2) building hierarchical structures to explicitly estimate IMVs for further N2O prediction; (3) using multi-task learning to balance the simultaneous training on multiple variables; and (4) pre-training with millions of synthetic data generated from ecosys and fine-tuning with mesocosm observations. Six other pure ML models were developed using the same mesocosm data to serve as the benchmark for the KGML-ag model. Results show that KGML-ag did an excellent job in reproducing the mesocosm N2O fluxes (overall r2=0.81, and RMSE=3.6 mgNm-2d-1 from cross validation). Importantly, KGML-ag always outperforms the PB model and ML models in predicting N2O fluxes, especially for complex temporal dynamics and emission peaks. Besides, KGML-ag goes beyond the pure ML models by providing more interpretable predictions as well as pinpointing desired new knowledge and data to further empower the current KGML-ag. We believe the KGML-ag development in this study will stimulate a new body of research on interpretable ML for biogeochemistry and other related geoscience processes.

54 ENVIRONMENTAL SCIENCES↗

Unbiased particle conformation extraction from scattering spectra using orthonormal basis expansions

A strategy is outlined for quantitatively evaluating the particle density profiles from small-angle scattering spectra of dilute solutions. The approach employs an orthonormal basis function expansion method, enabling the determination of characteristic mass distributions in self-assembled structures without the need for a specific structural model. Through computational benchmarking, the efficacy of this approach is validated by effectively reconstructing the density profile of soft-ball systems with varying fuzziness from their scattering signatures. Further, the feasibility of the method is demonstrated by fitting small-angle neutron scattering data obtained from Pluronic L64 micelles at different temperatures. This proposed approach is both simple and analytical, eliminating the requirement for a presumptive structural model in scattering analysis. The new method could therefore facilitate quantitative descriptions of complex nanoscopic structures inherent to numerous soft-matter systems using small-angle scattering techniques.

36 MATERIALS SCIENCE↗

Long-Range Allosteric Communication Modulated by Active Site Mn(II) Coordination Drives Catalysis in Xanthobacter autotrophicus Acetone Carboxylase

Acetone carboxylase (AC) from Xanthobacter autotrophicus is a 360 KDa α2β2γ2 heterohexamer that catalyzes the ATP-dependent formation of phosphorylated acetone and bicarbonate intermediates that react at Mn(II) metal active sites to form acetoacetate. Structural models of X. autotrophicus AC (XaAC) with and without nucleotides reveal that the binding and phosphorylation of the two substrates occurs ~40 Å from the Mn(II) active sites where acetoacetate is formed. Based on the crystal structures, a significant conformational change was proposed to open and close a tunnel that facilitates the passage of reaction intermediates between the sites for nucleotide binding and phosphorylation of substrates and Mn(II) sites of acetoacetate formation. We have employed electron paramagnetic resonance (EPR), kinetic assays, and hydrogen/deuterium exchange mass spectrometry (HDX-MS) of poised ligand-bound states and site-specific amino acid variants to complete an in-depth analysis of Mn(II) coordination and allosteric communication throughout the catalytic cycle. In contrast with the established paradigms for carboxylation, our analyses of XaAC suggested a carboxylate shift that couples both local and long-range structural transitions. Shifts in the coordination mode of a single carboxylic acid residue (αE89) mediate both catalysis proximal to a Mn(II) center and communication with an ATP active site in a separate subunit of a 180 kDa α2β2γ2 complex at a distance of 40 Å. This work demonstrates the power of combining structural models from X-ray crystallography with solution-phase spectroscopy and biophysical techniques to elucidate functional aspects of a multi-subunit enzyme.

Biochemistry & Molecular Biology↗