SEARCH · Engineering Papers

Results for “Fractionation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Search for η c ( 2 S ) → p p ¯ K + K − and measurement of χ c J → p p ¯ K + K − in ψ ( 3686 ) radiative decays

A search for η c ( 2 S ) → p p ¯ K + K − , together with measurement of branching fractions of χ c J ( J = 0 , 1 , 2 ) → p p ¯ K + K − in the ψ ( 3686 ) → γ η c ( 2 S ) and the ψ ( 3686 ) → γ χ c J radiative decays, is performed with ( 2712.4 ± 14.3 ) × 10 6 ψ ( 3686 ) events collected with the BESIII detector at the BEPCII collider. An evidence for η c ( 2 S ) → p p ¯ K + K − is found, with a significance of 3.3 σ . The product branching fraction of B [ ψ ( 3686 ) → γ η c ( 2 S ) ] · B [ η c ( 2 S ) → p p ¯ K + K − ] is determined to be ( 1.98 ± 0.4 1 stat ± 0.9 9 syst ) × 10 − 7 . The product branching fractions of B [ ψ ( 3686 ) → γ χ c J ] · B [ χ c J → p p ¯ K + K − ] are measured to be ( 2.49 ± 0.0 3 stat ± 0.1 5 syst ) × 10 − 5 , ( 1.83 ± 0.0 2 stat ± 0.1 1 syst ) × 10 − 5 , and ( 2.43 ± 0.0 2 stat ± 0.1 5 syst ) × 10 − 5 , for J = 0 , 1, and 2, respectively. Published by the American Physical Society 2025

Ablikim, M.↗

DOE OSTI2025

Study of η → π + π − l + l −

Using a sample of ( 10087 ± 44 ) × 10 6 J / ψ events accumulated with the BESIII detector, we analyze the decays η → π + π − l + l − ( l = e or μ ) via the process J / ψ → γ η . The branching fraction of η → π + π − e + e − is measured to be B ( η → π + π − e + e − ) = ( 3.07 ± 0.1 2 stat ± 0.1 9 syst ) × 10 − 4 . No signal events are observed for the η → π + π − μ + μ − decay, leading to an upper limit on the branching fraction of B ( η → π + π − μ + μ − ) < 4.0 × 10 − 7 at the 90% confidence level. Furthermore, the C P -violation asymmetry parameter is found to be A C P ( η → π + π − e + e − ) = ( − 4.04 ± 4.6 9 stat ± 0.1 4 syst ) % , showing no evidence of C P -violation with current statistics. Additionally, we extract the transition form factor from the decay amplitude of η → π + π − e + e − . Finally, axionlike particles are searched for via the decay η → π + π − a , a → e + e − , and upper limits on this branching fraction relative to that of η → π + π − e + e − are presented as a function of the axionlike particle mass in the range 5 – 200 MeV / c 2 . Published by the American Physical Society 2025

Ablikim, M.↗

DOE OSTI2025

Search for D 0 meson decays to π + π − e + e − and K + K − e + e − final states

A search for D 0 meson decays to the π + π − e + e − and K + K − e + e − final states is reported using a sample of proton-proton collisions collected by the LHCb experiment at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 6 fb − 1 . The decay D 0 → π + π − e + e − is observed for the first time when requiring that the two electrons are consistent with coming from the decay of a ϕ or ρ 0 / ω meson. The corresponding branching fractions are measured relative to the D 0 → K − π − [ e + e − ] ρ 0 / ω decay, where the two electrons are consistent with coming from the decay of a ρ 0 or ω meson. No evidence is found for the D 0 → K + K − e + e − decay and world-best limits are set on its branching fraction. The results are compared to, and found to be consistent with, the branching fractions of the D 0 → π + π − μ + μ − and D 0 → K + K − μ + μ − decays recently measured by LHCb and confirm lepton universality at the current precision. © 2025 CERN, for the LHCb Collaboration 2025 CERN

Aaij, R. (ORCID:0000000305331952)↗

DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

DOE OSTI2020

Unstable Invasion of Sedimenting Granular Suspensions

Here we investigate the development of mobility inversion and fingering when a granular suspension is injected radially between horizontal parallel plates of a cell filled with a miscible fluid. While the suspension spreads uniformly when the suspension and the displaced fluid densities are exactly matched, even a small density difference is found to result in a dense granular front which develops fingers with angular spacing that increase with granular volume fraction and decrease with injection rate. We show that the timescale over which the instability develops is given by the volume fraction dependent settling timescale of the grains in the cell. We then show that the mobility inversion and the nonequilibrium Korteweg surface tension due to granular volume fraction gradients determine the number of fingers at the onset of the instability in these miscible suspensions.

02 PETROLEUM↗

DOE OSTI2023

Z 3 Parafermion in the Double Charge Kondo Model

Quantum impurity models with frustrated Kondo interactions can support quantum critical points with fractionalized excitations. Recent experiments [arXiv:2108.12691] on a circuit containing two coupled metal-semiconductor islands exhibit transport signatures of such a critical point. Here we show using bosonization that the double charge-Kondo model describing the device can be mapped in the Toulouse limit to a sine-Gordon model. Its Bethe-ansatz solution shows that a Z 3 parafermion emerges at the critical point, characterized by a fractional $\frac{1}{2}$ln(3) residual entropy, and scattering fractional charges e/3. Finally, we also present full numerical renormalization group calculations for the model and show that the predicted behavior of conductance is consistent with experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Observation of the Rare Decay of the 𝜂 Meson to Four Muons

A search for the rare η → μ + μ − μ + μ − double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb −1 . A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the η → μ + μ − decay as normalization, the branching fraction B(η → μ + μ − μ + μ − ) = [5.0 ± 0.8 (stat) ± 0.7 (syst) ± 0.7 (B 2μ )] × 10 −9 is measured, where the last term is the uncertainty in the normalization channel branching fraction. This work achieves an improved precision of over five orders of magnitude compared to previous results, leading to the first measurement of this branching fraction, which is found to agree with theoretical predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Here, coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to t 1/4 for surface diffusion and t 1/3 for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Reverse segregation and self-organization in inclined chute flows of bidisperse granular mixtures

In the usual segregation scenario for stable inclined chute flows of bidisperse mixtures of fine and coarse spherical particles, coarse particles rise toward the free surface, forming a coarse-rich region atop the flowing pile. Beyond a threshold coarse-to-fine diameter ratio of approximately 4, conversely, the weight of the coarse particles exceeds the segregation driving forces, causing individual coarse particles to sink within the pile and producing a reversed segregation state. However, an understanding of the collective evolution of the pile structure is still lacking when the particle diameter ratio exceeds 4 and the coarse-particle mass fraction is appreciable. To explore this broadly bidisperse limit, we perform discrete element method simulations considering mean particle diameter ratios of up to 8 and coarse-particle mass fractions spanning 0.1 to 0.9. The steady-state flow profiles reveal several intriguing behaviors that depend on the diameter ratio and mass fraction. These include a previously identified transition from usual to reverse segregation and a newfound tendency to self-organize into alternating coarse- and fine-rich particle layers stacked along the shear gradient direction, with layer thickness dictated by the coarse-particle diameter. A fuller understanding of segregation at this scale could pave the way for enhanced mixing or demixing techniques at the commercial scale.

granular flow↗

DOE OSTI2023

A New Proposal Generalized Predictive Control Algorithm With Polynomial Reference Tracking Applied for Sodium Fast Reactors

This paper proposes a generalized predictive control (GPC) with constraints and orthonormal Laguerre functions using the simplified model of the primary system (reactor core and intermediate heat exchanger (IHX)) of a prototypical sodium fast reactor (SFR). This paper develops a multiple-input multiple-output (MIMO) GPC with input constraints able to track polynomial references of any degree applied in coolant temperature difference across the core and fractional power. The manipulated variables of the GPC-SFR are the reactivity and the sodium flow rate of the primary and secondary pipes. Moreover, orthonormal Laguerre functions and step down condition number techniques were also applied to avoid the numerical ill-conditioning issue in quadratic programming of large systems. Thus, a GPC type-2 was designed to control fractional power, coolant temperature difference across the core and sodium tank temperature of the SFR primary system when temperature references change according to a linear ramp after reaching their steady-state operation, sustaining 100% power operation on the reactor. In order to analyze the load tracking capability of the GPC-SFR type-2, the load following from 100% fractional power (FP) to 60% FP at 0.8% FP/min rate is simulated. Constraints on the rate of coolant temperature difference across the core and reactivity were applied for the design safety. For comparison criteria, this paper compares the GPC-SFR type-2 with the GPC-SFR type-1, i.e., standard model predictive control (MPC), to verify the viability and superior performance of the proposal regarding: (a) ramp-tracking capability of temperature and load; (b) the rejections of a reactivity disturbance of -1 cent and a secondary sodium inlet temperature disturbance of +10°F; and (c) a simulation with uncertainty in reactor design. The simulations show that the GPC-SFR type-2 overcome the GPC-SFR type-1 robustness and performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Temperature-Dependent Magnetic Properties of Magnetorheological Elastomers

We report on an investigation of the temperature-dependent magnetic properties of magnetorheological elastomers (MREs). These are a class of composites that consist of magnetically permeable particles dispersed in a nonmagnetic polymeric matrix. Under the application of an external magnetic field, a large deformation occurs altering the mechanical properties of these materials. Due to their magnetoelastic coupling response, these materials are finding an increasing interest among the scientific community. These polymer-based composites’ performance depends on many factors, which temperature is one of the biggest influencing factors requiring further investigation. In this work, the magnetic properties of isotropic and anisotropic polyurethane-based MRE with different iron (Fe) particle loading fractions (50, 60, 70, and 80% by weight) were investigated under different temperatures. From the analysis, the magnetization curves of these materials are observed to overlap for the different measured temperature values. The variation of various magnetic properties including saturation magnetization and differential susceptibility with temperature was also determined. Furthermore, the results show ~ 2.5% decrement in the saturation magnetization for each of the loading fractions, between the lowest (300K) and the highest (400K) measured temperature. On the other hand, the initial differential susceptibility exhibits different trends with increasing temperature. Generally, the magnetization response of these materials is seen to be only slightly sensitive to temperature changes. Additionally, the magnetization response is observed to be highly dependent on particle loading fractions and particle orientation within the elastomer.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Carbon flow through energycane agroecosystems established post-intensive agriculture

As part of an integrated energy and climate system, biomass production for bioenergy based on the tropical perennial C4 grass energycane can both offset fossil fuels and store soil carbon (C). We measured energycane yields, root biomass, soil C pools, and soil C stocks in a 4 year field trial and modeled C flow from plants to soils in the surface layer of no-till energycane planted after more than a century of intensive sugarcane agriculture. Aboveground yields ranged from 16.7 to 19.0 Mg C/ha over the 4 year trial. Although total C stocks did not significantly differ in the surface layer (approx. 0–20 cm) during the study, C in free and occluded light fractions decreased, whereas C in the mineral-rich dense fraction increased over 4 years. Belowground system inputs, estimated from measurements and informed by convergence in the final soil fraction model, were set to 2.5 Mg C ha -1 year -1 . With this input value, we estimated that surface soils retained photosynthetically fixed C predominantly within the mineral-associated organic matter pool for a mean and median transit time of 177 and 110 years, respectively. Although we did not model C flow to deep soil layers (approx. 0–100 cm), observed C accumulation (11.4 Mg C ha -1 year -1 ) and root growth down to 120 cm suggest that soil processes and resulting C sequestration at the surface are likely to persist deeper into the soil profile. Energycane, as a strong candidate for climate change mitigation and land degradation remediation, showed high biomass yields and allocation of resources to roots, with sequestered soil C expected to persist for over a century.

09 BIOMASS FUELS↗

DOE OSTI2021

Functional Adaptation of LPS-affected Dentoalveolar Fibrous Joints in Rats

We report the functional interplay between cementum of the root and alveolar bone of the socket is tuned by a uniquely positioned 70–80 µm wide fibrous and lubricious ligament in a dentoalveolar joint (DAJ). In this study, structural and biomechanical properties of the DAJ, periodontal ligament space (PDL-space also known as the joint space), alveolar bone of the socket, and cementum of the tooth root that govern the biomechanics of a lipopolysaccharide (LPS)-affected DAJ were mapped both in space and time. The hemi-maxillae from 20 rats (4 control at 6 weeks of age, 4 control and 4 LPS-affected at 12 weeks of age, 4 control and 4 LPS-affected at 16 weeks of age) were investigated using a hybrid technique; micro-X-ray computed tomography (5 µm resolution) in combination with biomechanical testing in situ. Temporal variations in bone and cementum volume fractions were evaluated. Trends in mineral apposition rates (MAR) in additional six Sprague Dawley rats (3 controls, 3 LPS-affected) were revealed by transforming spatial fluorochrome signals to functional growth rates (linearity factor - RW) of bone, dentin, and cementum using a fast Fourier transform on fluorochrome signals from 100-µm hemi-maxillae sections. An overall change in LPS-affected DAJ biomechanics (a 2.5–4.5X increase in tooth displacement and 2X tooth rotation at 6 weeks, no increase in displacement and a 7X increase in rotation at 12 weeks; 27% increase in bone effective strain at 6 weeks and 11% at 12 weeks relative to control) was associated with structural changes in the coronal regions of the DAJ (15% increase in PDL-space from 0 to 6 weeks but only 5% from 6 to 12 weeks compared to control). A significant increase (p < 0.05) in PDL-space between ligated and age-matched control was observed. The bone fraction of ligated at 12 weeks was significantly lower than its age-matched control, and no significant differences (p > 0.05) between groups were observed at 6 weeks. Cementum in the apical regions grew faster but nonlinearly (11% and 20% increase in cementum fraction (CF) at 6 and 12 weeks) compared to control. Alveolar bone revealed site-specific nonlinear growth with an overall increase in MAR (108.5 µm/week to 126.7 µm/week after LPS treatment) compared to dentin (28.3 µm/week in control vs. 26.1 µm/week in LPS-affected) and cementum (126.5 µm/week in control vs. 119.9 µm/week in LPS-affected). A significant increase in CF (p < 0.05) in ligated specimens was observed at 6 weeks of age. Anatomy-specific responses of cementum and bone to the mechano-chemo stimuli, and their collective temporal contribution to observed changes in PDL-space were perpetuated by altered tooth movement. Data highlight the “resilience” of DAJ function through the predominance of nonlinear growth response of cementum, changes in PDL-space, and bone architecture. Despite the significant differences in bone and cementum architectures, data provided insights into the reactionary effects of cementum as a built-in compensatory mechanism to reestablish functional competence of the DAJ. The spatial shifts in architectures of alveolar bone and cementum, and consequently ligament space, highlight adaptations farther away from the site of insult, which also is another novel insight from this study. These adaptations when correlated within the context of joint function (biomechanics) illustrate that they are indeed necessary to sustain DAJ function albeit being pathological.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

A Double-Wiebe Function for Reactivity Controlled Compression Ignition Combustion Using Reformate Diesel

Abstract Reactivity controlled compression ignition (RCCI) combustion has previously been proposed as a method to achieve high fuel conversion efficiency and reduce engine emissions. A single-fuel RCCI combustion strategy can have decreased fuel system complexity by using a reformate fuel for port fuel injection and the parent fuel (diesel) for direct injection. This paper presents a one-dimensional computational model of a compression ignition engine with single-fuel RCCI. A Wiebe function is used to predict the combustion process by representing the mass fraction burned (MFB) on a crank angle resolved basis. One single-Wiebe function (SWF) and two double-Wiebe functions(DWFs) were fitted to experimentally derive MFB data using the least-square method. The fitted results were compared with MFBs calculated from experimental data to verify the accuracy. The SWF did not fully capture the MFB curve with high fidelity while the detailed DWF captured the MFB curve within a root mean square error of 1.4%. The reduced double-Wiebe function (RDWF) also resulted in a predicted combustion profile with similar accuracy. Hence, the RDWF was used in a GT-power thermodynamic study to understand the effects of the low-temperature heat release (LTHR) fraction and combustion phasing on combustion characteristics. At optimum phasing of 5–10 crank angle degree after the top dead center, increasing the LTHR fraction from 20% to 60% resulted in the fuel conversion efficiency increasing from 39.5% to 41.1%, thus suggesting that the reformate fuel-based RCCI strategy is viable to unlock improved combustion performance.

Energy & Fuels↗

DOE OSTI2023

MACHINE LEARNING-ENABLED PREDICTION OF TRANSIENT INJECTION MAP IN AUTOMOTIVE INJECTORS WITH UNCERTAINTY QUANTIFICATION

Accurate prediction of injection profiles is a critical aspect of linking injector operation with engine performance and emissions. However, highly resolved injector simulations can take one to two weeks of wall-clock time, which is incompatible with engine design cycles with desired turnaround times of less than a day. Hence, it is important to reduce the time-to-solution of the internal flow simulations by several orders of magnitude to make it compatible with engine simulations. This work demonstrates a data-driven approach for tackling the computational overhead of injector simulations, whereby the transient injection profiles are emulated for a side-oriented, single-hole diesel injector using a Bayesian machine-learning framework. First, an interpretable Bayesian learning strategy was employed to understand the effect of design parameters on the total void fraction field. Then, autoencoders are utilized for efficient dimensionality reduction of the flowfields. Gaussian process models are finally used to predict the spatiotemporal void fraction field at the injector exit for unknown operating conditions. The Gaussian process models produce principled uncertainty estimates associated with the emulated flowfields, which provide the engine designer with valuable information of where the data-driven predictions can be trusted in the design space. The Bayesian flowfield predictions are compared with the corresponding predictions from a deep neural network, which has been transfer-learned from static needle simulations from a previous work by the authors. The emulation framework can predict the void fraction field at the exit of the orifice within a few seconds, thus achieving a speed-up factor of up to 38 x 10(6) over the traditional simulation-based approach of generating transient injection maps.

machine learning↗

DOE OSTI2023

A Monte-Carlo Method for Evaluating the Economic Performance of Plastics Recycling Systems Using Historical Pricing

The management of plastic waste is a considerable emerging global concern. Conventional plastics recycling rates range from 5–10% in the US, which results in considerable material being sent to landfill. New sorting technologies using artificial intelligence have the potential to dramatically increase the ability of waste management companies to sort out valuable fractions of plastic waste and create high-purity streams for secondary markets. However, several fractions of the plastics recycling stream, i.e., low-density polyethylene, polypropylene, and polystyrene, do not have well-developed markets. These materials may be sent to pyrolysis facilities for thermal recycling, but the economics of the pyrolysis industry are currently uncertain. This work aims to identify the breakeven price of the plastic fractions that would be sent to pyrolysis in order to determine the economic viability of the sorting facility. The work will use Monte Carlo analysis to evaluate several scenarios and understand how compositional variation, price variability, and facility attributes affect the breakeven price. Here, the presorting conditions in Materials Recovery Facilities strongly affects the profitability of the sorting facility. High residual high density polyethylene (HDPE) and polyethylene terephthalate (PET) content in presorted bales are important revenue drivers for the facility, even in presorted 3-7 bales. Key variables that lead to profitability include a high sale price of polyethylene terephthalate, a high sale price of high-density polyethylene, and if polypropylene can be sold to recycling markets (as opposed to pyrolysis facilities).

Plastics↗

DOE OSTI2020

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗

DOE OSTI2021

Carbon isotope evidence for the global physiology of Proterozoic cyanobacteria

Ancestral cyanobacteria are assumed to be prominent primary producers after the Great Oxidation Event [≈2.4 to 2.0 billion years (Ga) ago], but carbon isotope fractionation by extant marine cyanobacteria (α-cyanobacteria) is inconsistent with isotopic records of carbon fixation by primary producers in the mid-Proterozoic eon (1.8 to 1.0 Ga ago). To resolve this disagreement, we quantified carbon isotope fractionation by a wild-type planktic β-cyanobacterium (Synechococcus sp. PCC 7002), an engineered Proterozoic analog lacking a CO 2 -concentrating mechanism, and cyanobacterial mats. At mid-Proterozoic pH and pCO 2 values, carbon isotope fractionation by the wild-type β-cyanobacterium is fully consistent with the Proterozoic carbon isotope record, suggesting that cyanobacteria with CO 2 -concentrating mechanisms were apparently the major primary producers in the pelagic Proterozoic ocean, despite atmospheric CO 2 levels up to 100 times modern. The selectively permeable microcompartments central to cyanobacterial CO 2 -concentrating mechanisms (“carboxysomes”) likely emerged to shield rubisco from O 2 during the Great Oxidation Event.

59 BASIC BIOLOGICAL SCIENCES↗