Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “carbon dynamics”

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 631 records · Page 35

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗

Carbon Mineralization of Sulfate Wastes Containing Pb: Synchrotron Pb M3-Edge XANES Analysis of Simultaneous Heavy Metal and Carbon Sequestration

Sulfate wastes are produced in large quantities and contain toxic heavy metals such as lead (Pb), posing environmental risks. Because of favorable solubility differences, these wastes can be repurposed for engineered carbon dioxide (CO 2 ) sequestration. Understanding the fate and mobility of heavy metals during this process is important. This study focuses on Pb and the effect of zinc (Zn) on Pb in carbon mineralization. Synthesized gypsum was treated with a carbonate-rich solution at pH 11.5 to convert the sulfates to carbonates. Aqueous solutions and mineral solids were analyzed. Synchrotron-based micro-X-ray fluorescence and a novel application of Pb M3-edge X-ray absorption near-edge structure provided detailed insights into Pb distribution and mineral forms. Results showed significant reductions in aqueous Pb and Zn concentrations, indicating effective metal sequestration. Carbon mineralization transformed Pb from soluble anglesite (PbSO 4 ) into insoluble cerussite (PbCO 3 ) and hydrocerussite (Pb 3 (CO 3 ) 2 (OH) 2 ). Pb primarily precipitated onto calcium carbonate surfaces through surface-mediated precipitation reactions. While the presence of Zn modified crystallization dynamics, it did not impede Pb sequestration and potentially enhanced surface reactivity, facilitating greater Pb immobilization. These findings highlight carbon mineralization as a sustainable approach to immobilize toxic metals in sulfate wastes while advancing CO 2 sequestration efforts.

54 ENVIRONMENTAL SCIENCES↗

Infrared Fingerprint and Unimolecular Decay Dynamics of the Hydroperoxyalkyl Intermediate (•QOOH) in Cyclopentane Oxidation

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C 5 H 9 O 2 ), building on higher level reference calculations for ethane oxidation (C 2 H 5 O 2 ). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. Furthermore, the various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Chemical calculations↗

Reservoir Dynamics in Proposed Operational Scenarios Where CO 2 -EOR Fields are Transitioned From CO 2 -Flood Enhanced Oil Recovery to Dedicated Carbon Storage: A Field Case Study

This study models the transition of CO 2 -enhanced oil recovery (CO 2 -EOR) fields to dedicated carbon storage using a generalized reservoir model informed by real field data, SACROC. Through scenario-based simulations, the work evaluates how reservoir depletion levels, boundary conditions, fluid properties, and domain size affect pressure buildup, CO 2 plume migration, and long-term containment performance. The results inform practices for repurposing oil fields into reliable CO 2 storage sites.

02 PETROLEUM↗

Terrestrial photosynthesis inferred from plant carbonyl sulfide uptake

Terrestrial photosynthesis, or gross primary production (GPP), is the largest carbon flux in the biosphere, but its global magnitude and spatiotemporal dynamics remain uncertain. The global annual mean GPP is historically thought to be around 120 PgC yr −1 , which is about 30–50 PgC yr −1 lower than GPP inferred from the oxygen-18 ( 18 O) isotope and soil respiration. This disparity is a source of uncertainty in predicting climate–carbon cycle feedbacks. Here, in this study, we infer GPP from carbonyl sulfide, an innovative tracer for CO 2 diffusion from ambient air to leaf chloroplasts through stomata and mesophyll layers. We demonstrate that explicitly representing mesophyll diffusion is important for accurately quantifying the spatiotemporal dynamics of carbonyl sulfide uptake by plants. From the estimate of carbonyl sulfide uptake by plants, we infer a global contemporary GPP of 157 (±8.5) PgC yr −1 , which is consistent with estimates from 18 O (150–175 PgC yr −1 ) and soil respiration ($149^{+29}_{-23}$) PgC yr −1 ), but with an improved confidence level. Our global GPP is higher than satellite optical observation-driven estimates (120–140 PgC yr –1 ) that are used for Earth system model benchmarking. This difference predominantly occurs in the pan-tropical rainforests and is corroborated by ground measurements, suggesting a more productive tropics than satellite-based GPP products indicated. As GPP is a primary determinant of terrestrial carbon sinks and may shape climate trajectories, our findings lay a physiological foundation on which the understanding and prediction of carbon–climate feedbacks can be advanced.

Lai, Jiameng↗

Extension of Plant Dynamics Code Capabilities for Simulation of TerraPower Pascal Reactor

Argonne National Laboratory has been developing the Plant Dynamics Code (PDC) for design and transient analysis of supercritical carbon dioxide (sCO 2 ) Brayton cycles. In previous analyses with PDC, only indirect sCO 2 cycles, where the heat is being added through a heat exchanger, such as sodium-to-CO 2 HX, were analyzed. Under the U.S. Department of Energy Technology Commercialization Fund (TCF), Argonne cooperated with TerraPower to bring the Plant Dynamics Code to commercial market. The main focus of the TCF project is to extend the application base and the code usability by developing the capabilities to be able to simulate reactor systems with direct sCO 2 cycles. These new PDC capabilities have been developed in application to the TerraPower’s Pascal reactor concept. This report documents the Pascal reactor modeling with the PDC, including simulation of the Pascal split-expansion cycle, the development of the reactor module in the PDC, modeling of the Pascal’s shutdown heat removal system in PDC, and other updates to the code. The report also describes the results of the steady state and transient demonstration of the newly developed code features.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Bulk density, grain size, carbon content and isotopes, and nitrogen content of floodplain sediment along the Koyukuk River, Alaska collected in Summer 2018.

The dataset is comprised of grain size, bulk density, carbon and nitrogen content, stable carbon isotopes, and radiocarbon content of soils sampled on the floodplain of the Koyukuk River near Huslia, AK. This dataset was collected to characterize sediment properties and floodplain organic carbon stocks in the discontinuous permafrost floodplain of the Koyukuk River. The data help quantify the flux and storage of sediment and carbon between the Koyukuk River and the surrounding floodplain as the river migrates laterally. The samples were collected in June 27 - July 9, 2018. The soil samples come from a total of 24 locations and were sampled in the upper 490 cm of the river cutbanks, point bars, and floodplain. Samples were analyzed for bulk density, full particle size distribution, nitrogen content and carbon content and isotopes, though not all samples were included in each analysis. An additional 10 samples of floodplain woody debris were measured for radiocarbon content. Data, detailed methods, analyses, and interpretations were previously published at: Douglas, M. M., Lingappa, U. F., Lamb, M. P., Rowland, J. C., West, A. J., Li, G., et al. (2021). Impact of river channel lateral migration on microbial communities across a discontinuous permafrost floodplain. Applied and Environmental Microbiology, 87(20), AEM.01339-21. https://doi.org/10.1128/AEM.01339-21 Douglas, M. M., Li, G. K., Fischer, W. W., Rowland, J. C., Kemeny, P. C., West, A. J., et al. (2022). Organic carbon burial by river meandering partially offsets bank erosion carbon fluxes in a discontinuous permafrost floodplain. Earth Surface Dynamics, 10(3), 421–435. https://doi.org/10.5194/esurf-10-421-2022 The dataset is comprised of a total of three separate data files based on the sample type and analysis, including a file that lists the sample locations, permafrost presence, and geomorphic landform.

54 ENVIRONMENTAL SCIENCES↗

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials (Final Technical Report)

The goal of this project was to conduct dynamic compression experiments and predictive simulations to reveal the fundamental high-energy-density (HED) physics of amorphous carbon. These results are essential for assessing amorphous carbon as a potential ablator material for next generation of inertial confinement fusion (ICF) capsules. We made significant progress in exploring the HED properties of amorphous carbon through experiments at Omega EP Laser and the European XFEL, in addition to billion-atom, quantum-accurate molecular dynamics (MD) simulations. Through our joint experimental and simulation program, we mapped the phase diagram of amorphous carbon, uncovering its range of metastability and identifying phase transitions to diamond and liquid carbon along the Hugoniot and at higher pressures using double shock compression pathways. Our findings indicate that amorphous carbon melts at significantly lower shock pressures than high-density carbon (diamond). However, nanocrystalline diamond nucleates across a broad range of pressures and temperatures. This emergence of the nanocrystalline microstructure during compression can negatively impact the planarity of the shock front and potentially trigger ablator/fuel mixing during Inertial Fusion Energy (IFE) applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion Solvation and Transport in Narrow Carbon Nanotubes: Effects of Polarizability, Cation-π Interaction, and Confinement

Understanding ion solvation and transport under confinement is critical for a wide range of emerging technologies, including water desalination and energy storage. While molecular dynamics (MD) simulations have been widely used to study the behavior of confined ions, considerable deviations between simulation results depending on the specific treatment of intermolecular interactions remain. In the following, we present a systematic investigation of the structure and dynamics of two representative solutions, that is, KCl and LiCl, confined in narrow carbon nanotubes (CNTs) with a diameter of 1.1 and 1.5 nm, using a combination of first-principles and classical MD simulations. Our simulations show that the inclusion of both polarization and cation-π interactions is essential for the description of ion solvation under confinement, particularly for large ions with weak hydration energies. Beyond the variation in ion solvation, we find that cation-π interactions can significantly influence the transport properties of ions in CNTs, particularly for KCl, where our simulations point to a strong correlation between ion dehydration and diffusion. Finally, our study highlights the complex interplay between nanoconfinement and specific intermolecular interactions that strongly control the solvation and transport properties of ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A scalable solid-state nanoporous network with atomic-level interaction design for carbon dioxide capture

Carbon capture and sequestration reduces carbon dioxide emissions and is critical in accomplishing carbon neutrality targets. Here, we demonstrate new sustainable, solid-state, polyamine-appended, cyanuric acid–stabilized melamine nanoporous networks (MNNs) via dynamic combinatorial chemistry (DCC) at the kilogram scale toward effective and high-capacity carbon dioxide capture. Polyamine-appended MNNs reaction mechanisms with carbon dioxide were elucidated with double-level DCC where two-dimensional heteronuclear chemical shift correlation nuclear magnetic resonance spectroscopy was performed to demonstrate the interatomic interactions. We distinguished ammonium carbamate pairs and a mix of ammonium carbamate and carbamic acid during carbon dioxide chemisorption. The coordination of polyamine and cyanuric acid modification endows MNNs with high adsorption capacity (1.82 millimoles per gram at 1 bar), fast adsorption time (less than 1 minute), low price, and extraordinary stability to cycling by flue gas. This work creates a general industrialization method toward carbon dioxide capture via DCC atomic-level design strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Nature of Mineral-Organic Associations within Redox-Active Mountainous Floodplain Sediments

Floodplains are critical terrestrial-aquatic interfaces that act as hotspots of organic carbon (OC) cycling, regulating ecosystem carbon storage as well as export to riverine systems. Within floodplain sediments, regular flooding and textural gradients interact to create dynamic redox conditions. While anaerobic protection of OC upon burial is a well-recognized carbon storage mechanism in redox-active floodplain sediments, the impact of protective mineral-organic associations is relatively unknown. Here we determined the quantitative importance and chemical composition of mineral-organic associations along well-defined redox gradients emerging from textural variations and depth within meander sediments of the subalpine East River watershed (Gothic, CO). Here we characterized mineral-organic associations using a combination of sequential extractions, physical fractionation, and high-resolution mass spectrometry. Across the meander, we found that mineral-associated OC constitutes a significant fraction of total OC, and that extractable iron (Fe) and aluminum (Al) phases as well as high-density isolates were strongly correlated with total OC content, suggesting that mineral-organic associations are quantitatively important for floodplain sediment OC protection. Our mass spectrometry results showed OC associated with increasingly ordered Fe and Al phases are relatively enriched in low-molecular weight, oxidized, aromatic compounds. Surprisingly, however, total OC content showed weak or no correlation with indicators of anaerobic protection, such as relatively bioavailable OC pools (water-extractable and particulate OC) or the molecular weight and oxidation state of OC. Overall, this work highlights that protection of OC bound to reactive mineral phases-in addition to anaerobic protection-can play a quantitatively important role in controlling soil carbon storage in redox-active floodplain sediments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Trends and Drivers of Terrestrial Sources and Sinks of Carbon Dioxide: An Overview of the TRENDY Project

The terrestrial biosphere plays a major role in the global carbon cycle, and there is a recognized need for regularly updated estimates of land-atmosphere exchange at regional and global scales. An international ensemble of Dynamic Global Vegetation Models (DGVMs), known as the “Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide” (TRENDY) project, quantifies land biophysical exchange processes and biogeochemistry cycles in support of the annual Global Carbon Budget assessments and the REgional Carbon Cycle Assessment and Processes, phase 2 project. DGVMs use a common protocol and set of driving data sets. A set of factorial simulations allows attribution of spatio-temporal changes in land surface processes to three primary global change drivers: changes in atmospheric CO 2 , climate change and variability, and Land Use and Land Cover Changes (LULCC). Here, we describe the TRENDY project, benchmark DGVM performance using remote-sensing and other observational data, and present results for the contemporary period. Simulation results show a large global carbon sink in natural vegetation over 2012–2021, attributed to the CO 2 fertilization effect (3.8 ± 0.8 PgC/yr) and climate (–0.58 ± 0.54 PgC/yr). Forests and semi-arid ecosystems contribute approximately equally to the mean and trend in the natural land sink, and semi-arid ecosystems continue to dominate interannual variability. The natural sink is offset by net emissions from LULCC (–1.6 ± 0.5 PgC/yr), with a net land sink of 1.7 ± 0.6 PgC/yr. Despite the largest gross fluxes being in the tropics, the largest net land-atmosphere exchange is simulated in the extratropical regions.

58 GEOSCIENCES↗

Impact of Open Communication Networks on Load Frequency Control with Plug-In Electric Vehicles by Cyber-Physical Dynamic Co-Simulation

With the increasing electrification of the transportation sector to achieve the carbon neutrality objective, despite the challenges of charging electric vehicles (EV), there are also opportunities through smart charging EVs to improve system frequency stability; however, EV control technologies might require nontraditional communication support. This paper investigates the impacts of communication variations of EV on power system load frequency control through a cyber-physical dynamic system (CPDS) co-simulation. Here, the CPDS is built upon our previously developed transmission-and-distribution dynamic co-simulation model with the added communication variation functions (i.e., delay and packet loss). The case studies consider multiple communication variation scenarios when the system experiences an N-1 generation trip contingency. The scenarios include communication delays and packet loss using both homogeneous and heterogeneous assumptions. The outcomes of this work can help improve EV frequency regulation services and provide robust and effective tests for different load frequency control algorithms of the future power systems.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗

Spatial access and resource limitations control carbon mineralization in soils

Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale SOM mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble C (acetate) and N (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Respiration showed evidence of C-limitation, as C amendments increased the cumulative CO 2 evolved. Drought and rewetting increased soil respiration, and these soils also exhibited increased complex aromatic molecules in porewater. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We expected that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (porewater samples), compared to the baseline soils, drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules. This indicates that as soils were wet from below, the fresh inputs simulated the C-and N-limited microbial populations present deeper in the soil profile. Furthermore, our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of Fluoroethylene Carbonate Additives on the Initial Formation of the Solid Electrolyte Interphase on an Oxygen-Functionalized Graphitic Anode in Lithium-Ion Batteries

The formation of a solid electrolyte interphase (SEI) at the electrode/electrolyte interface substantially affects the stability and lifetime of lithium-ion batteries (LIBs). One of the methods to improve the lifetime of LIBs is by the inclusion of additive molecules to stabilize the SEI. To understand the effect of additive molecules on the initial stage of SEI formation, in this work we compare the decomposition and oligomerization reactions of a fluoroethylene carbonate (FEC) additive on a range of oxygen-functionalized graphitic anodes to those of an ethylene carbonate (EC) organic electrolyte. A series of density functional theory (DFT) calculations augmented by ab initio molecular dynamics (AIMD) simulations reveal that EC decomposition on an oxygen-functionalized graphitic ($11\bar{2}0$) edge facet through a nucleophilic attack on an ethylene carbon site (C E ) of an EC molecule (S2 mechanism) is spontaneous during the initial charging process of LIBs. However, decomposition of EC through a nucleophilic attack on a carbonyl carbon (C C ) site (S1 mechanism) results in alkoxide species regeneration that is responsible for continual oligomerization along the graphitic surface. In contrast, FEC prefers to decompose through an S1 pathway, which does not promote alkoxide regeneration. Including FEC as an additive is thus able to suppress alkoxide regeneration and results in a smaller and thinner SEI layer that is more flexible toward lithium intercalation during the charging/discharging process. In addition, we find that the presence of different oxygen functional groups at the surface of graphite dictates the oligomerization products and the LiF formation mechanism in the SEI.

25 ENERGY STORAGE↗

Atmosphere-Biosphere Interaction Study (Field Campaign Report)

The objective of this campaign was to collect hyperspectral reflectance measurements and complementary optical and infrared data of land surface (i.e., vegetation and soil) throughout the growing season (i.e., from March to September) at high temporal frequency (e.g., < 30 minutes). To collect the measurements, the EcoSpec system, a tower-based sensor system housing seven optical and infrared sensors that was developed at Argonne National Laboratory, was deployed to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) atmospheric observatory (within the crop, south of the eddy correlation flux measurement system [ECOR] tower). The EcoSpec system consists of a spectroradiometer having 2,151 channels (350–2,500 nm), red-greenblue (RGB) camera, diffuse radiometer, thermal infrared (TIR) sensors, albedometers, Photochemical Reflectance Index (PRI) sensors, Normalized Difference Vegetation Index (NDVI) sensors, and wetness sensor. The system is designed to collect an array of land surface measurements, including hyperspectral reflectance of plants and soils, RGB photos capturing contextual information of land surface, radiant temperature of sky and land surface, incoming and reflected shortwave radiation, and incoming radiation components (direct and diffuse) as frequent as every minute from up to 12 discrete positions around the tower. This would allow us to capture inherent heterogeneity of the terrestrial ecosystem and to obtain representative spectral signatures of the plant and soil surfaces via averaging those measurements. The SGP site provides advanced atmospheric and other physical measurements that are essential for our investigation of ecosystem functions using hyperspectral reflectance measurements of land surfaces. Located in the mid-latitudes, the SGP site exhibits dynamic seasonality and a wide range of meteorological variability and flux properties that make it ideal for studying interactions between near-surface atmosphere and biosphere across a range of environmental conditions. The atmosphere, plants, and soil control terrestrial carbon and water cycles. To make more accurate climate forecasts, researchers need to understand ecosystem dynamics at the biosphere–atmosphere interface. Even though our ability to forecast ecosystem dynamics and climate at both the regional and global scales has advanced significantly, we still do not know how local phenomena affecting water and carbon fluxes—such as diurnal variations—interrelate with large-scale atmospheric/climate phenomena, and vice versa. The lack of a more thorough understanding of ecosystem dynamics and climate interactions is a key knowledge gap for the improvement of earth system models (ESMs).

54 ENVIRONMENTAL SCIENCES↗

Inference of photosynthetic capacity parameters from chlorophyll a fluorescence is affected by redox state of PSII reaction centers

Solar-induced chlorophyll fluorescence (SIF) has been used to infer photosynthetic capacity parameters (e.g., the maximum carboxylation rate V cmax , and the maximum electron transport rate J max ). However, the precise mechanism and practical utility of such approach under dynamic environments remain unclear. In this study, we used the balance between the light and carbon reactions to derive theoretical equations relating chlorophyll a fluorescence (ChlF) emission and photosynthetic capacity parameters, and formulated testable hypotheses regarding the dynamic relationships between the true total ChlF emitted from PSII (SIF PSII ) and V cmax and J max . We employed concurrent measurements of gas exchanges and ChlF parameters for 15 species from six biomes to test the formulated hypotheses across species, temperatures, and limitation state of carboxylation. Our results revealed that SIF PSII alone is incapable of informing the variations in V cmax and J max across species, even when SIF PSII is determined under the same environmental conditions. In contrast, the product of SIF PSII and the fraction of open PSII reactions q L , which indicates the redox state of PSII, is a strong predictor of both V cmax and J max , although their precise relationships vary somewhat with environmental conditions. Our findings suggest the redox state of PSII strongly influences the relationship between SIF PSII and V cmax and J max .

59 BASIC BIOLOGICAL SCIENCES↗