Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial evolution”

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 343 records · Page 19

Cycling Performance and Structure Evolution of Co-Free Lithium- and Manganese-Rich Layered Oxides in Lithium-Ion Batteries

Lithium- and manganese-rich (LMR) layered oxides are high-capacity cathode materials that are being considered for electric vehicle (EV) lithium-ion batteries (LIBs). Here, we investigate the electrochemical cycling behavior of cells containing a cobalt-free LMR oxide, 0.3Li 2 MnO 3 ·0.7LiMn 0.5 Ni 0.5 O 2 , paired with either Li or graphite anodes. Two- and three-electrode cells are cycled under varying conditions to examine the effects of oxide activation, upper cutoff voltage, separator type, and electrolyte composition. Post-cycling characterization of harvested electrodes, with electrochemical, X-ray absorption spectroscopy, and solid-state 6 Li nuclear magnetic resonance spectroscopy techniques, are used to correlate cell performance changes with redox state evolution in the oxide and Li inventory shifts in the cell. Our results show that capacity fade and impedance rise primarily originate at the graphite anode and LMR cathode, respectively. Notably, while Ni undergoes reversible redox, Mn shows no evidence of bulk redox activity even in highly-aged electrodes: only relithiated LMR electrodes show some reduction of Mn. Graphite electrodes degrade due to particle isolation and SEI growth, exacerbated by Mn dissolution and deposition on graphite anodes. These findings highlight the coupled structural and electrochemical degradation mechanisms in LMR/graphite cells and underscore lithium inventory retention and electrode interfacial stability as critical factors for enhancing long-term performance.

Badami, Pavan P. [Argonne National Laboratory (ANL↗

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, where droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Computational Fluid Dynamics↗

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, wherein droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Self-Pressurization↗

Development and Validation of Two-Phase Flow Models in MOOSE and Application to Molten Salt Reactors

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated a new two-phase model into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model. It provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model’s accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase model implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, the model is applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of the mixture drift-flux model for capturing the two-phase flow dynamics critical to MSR operations

42 - ENGINEERING↗

Reactive capture and electrochemical conversion of CO 2 with ionic liquids and deep eutectic solvents

Ionic liquids (ILs) and deep eutectic solvents (DESs) have tremendous potential for reactive capture and conversion (RCC) of CO 2 due to their wide electrochemical stability window, low volatility, and high CO 2 solubility. There is environmental and economic interest in the direct utilization of the captured CO 2 using electrified and modular processes that forgo the thermal- or pressure-swing regeneration steps to concentrate CO 2 , eliminating the need to compress, transport, or store the gas. The conventional electrochemical conversion of CO 2 with aqueous electrolytes presents limited CO 2 solubility and high energy requirement to achieve industrially relevant products. Additionally, aqueous systems have competitive hydrogen evolution. In the past decade, there has been significant progress toward the design of ILs and DESs, and their composites to separate CO 2 from dilute streams. In parallel, but not necessarily in synergy, there have been studies focused on a few select ILs and DESs for electrochemical reduction of CO 2 , often diluting them with aqueous or non-aqueous solvents. The resulting electrode–electrolyte interfaces present a complex speciation for RCC. In this review, we describe how the ILs and DESs are tuned for RCC and specifically address the CO 2 chemisorption and electroreduction mechanisms. Critical bulk and interfacial properties of ILs and DESs are discussed in the context of RCC, and the potential of these electrolytes are presented through a techno-economic evaluation.

36 MATERIALS SCIENCE↗

Exploration of Electrochemical Reactions at Organic–Inorganic Halide Perovskite Interfaces via Machine Learning in In Situ Time-of-Flight Secondary Ion Mass Spectrometry

The instability of hybrid organic–inorganic perovskite (HOIP) devices is one of the significant challenges preventing commercialization. Exploring these phenomena is severely limited by the complexity of the intrinsic electrochemistry of HOIPs, the presence of multiple volatile and mobile ionic species, and the possible role of environmentally induced reactions at surfaces and triple-phase junctions. Here, in situ studies of the electrochemistry of methylammonium lead bromide perovskite with the Au electrode interface are reported via light- and voltage-dependent time-of-flight secondary ion mass spectrometry (ToF-SIMS) imaging of lateral perovskite heterostructures. While ToF-SIMS allows for the visualization of the chemical composition along the surface and its evolution with light and electrical bias, the interpretation of the multidimensional data obtained is often limited due to strong correlations between chemical signatures and the need to track multiple peaks at once. Here, a machine learning workflow combining the Hough transform and non-negative matrix factorization and no-negative tensor decomposition is developed to avoid this limitation and extract salient features of associated chemical changes and to separate the light- and voltage-dependent dynamics. Finally, combining these in situ characterizations and the machine learning workflow provides comprehensive information on the chemical nature of moving species, ion accumulation, and interfacial electrochemical reactions in HOIP devices.

36 MATERIALS SCIENCE↗

Pt Nanoparticle Disintegration at Oxide Interfaces Enhances CO Oxidation Catalysis

Understanding how supported metal nanoparticles dynamically evolve under reaction conditions is critical for controlling their catalytic function. Here, the mechanism behind the dynamic disintegration of Pt nanoparticles (NPs) supported on CeOx-TiO2 (CT) during CO oxidation is elucidated, leading to the formation of single atoms (SAs) and/or sub-nanometer clusters. Density functional theory (DFT) calculations reveal that strong Pt-CO interactions weaken Pt─Pt cohesion, while electronic coupling between Pt and Ce ions stabilizes Pt-CO* intermediates at the oxide interface. Surface oxygen vacancies kinetically trap Pt-CO*, but the vacancies are replenished under oxygen-rich conditions, enabling Pt-CO* surface diffusion and subsequent structural reorganization. In situ spectroscopic analyses confirm the oxygen-driven transformation of Pt NPs, correlating with a threefold increase in mass-specific activity at 150 °C. These findings demonstrate that interfacial oxygen dynamics and metal-support interactions can be leveraged to induce nanoparticle disintegration and optimize catalytic performance, highlighting the catalytic potential of interface-engineered Pt nanostructures.

CO oxidation↗

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗

Identification of localized radiation damage in power MOSFETs using EBIC imaging

The rapidly increasing use of electronics in high-radiation environments and the continued evolution in transistor architectures and materials demand improved methods to characterize the potential damaging effects of radiation on device performance. Here, electron-beam-induced current is used to map hot-carrier transport in model metal-oxide semiconductor field-effect transistors irradiated with a 300 KeV focused He + beam as a localized line spanning across the gate and bulk Si. By correlating the damage to the electronic properties and combining these results with simulations, the contribution of spatially localized radiation damage on the device characteristics is obtained. Furthermore, this identified damage, caused by the He + beam, is attributed to localized interfacial P b centers and delocalized positive fixed-charges, as surmised from simulations. Comprehension of the long-term interaction and mobility of radiation-induced damage are key for future design of rad-hard devices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Oxidation Driven Thin‐Film Solid‐State Metal Dealloying Forming Bicontinuous Nanostructures

Abstract Solid‐state metal dealloying (SSMD) is a promising method for fabricating nanoscale metallic composites and nanoporous metals across a range of materials. Thin‐film SSMD is particularly attractive due to its ability to create fine features via solid‐state interfacial reactions within a thin‐film geometry, which can be integrated into devices for various applications. This work examines a new dealloying couple, namely the Nb–Al alloy with the dealloying agent Sc, as previously predicted in the machine‐learning (ML) models. Prior ML predictions aimed to guide the design of nanoarchitectured materials through dealloying, relying on intuition‐driven discovery within a large parameter space. However, this work reveals that at the nanoscale, the involvement of oxygen in thin film processing may instead drive the dealloying process, resulting in the formation of bicontinuous nanostructures similar to those formed by metal‐agent dealloying. The phase evolution, as well as chemical and morphological changes, are closely analyzed using a combination of X‐ray absorption spectroscopy, diffraction, and scanning transmission electron microscopy to understand the mechanisms behind nanostructure formation. The findings suggest a potential pathway for utilizing oxygen to drive the formation of bicontinuous metal–metal oxide nanocomposites, paving the way for further development of functional nanoporous materials in diverse fields.

36 MATERIALS SCIENCE↗

Advancing Photo(Electro)Chemical Water Splitting: The Promise of Atomically Dispersed Single‐, Dual‐, and Alloy‐Site Catalysts

Single-atom catalysts (SACs) have rapidly gained prominence as an emerging class of electrocatalysts for water splitting, owing to their uniform and precisely defined active sites. By maintaining uniform reaction pathways, SACs minimize the formation of unwanted byproducts, thus exhibiting extremely high selectivity and atomic efficiency. A key determinant of SAC performance lies in the interfacial interaction between the isolated metal atoms and the supporting material under strong metal–support coordination, which is vital for maintaining long-term activity. However, despite these benefits, reproducibly synthesizing SACs with high metal loadings while retaining uniform dispersion remains a significant challenge. To address the intrinsic challenges of SACs, recent research has expanded into dual-atom catalysts (DACs) and single-atom alloy catalysts (SAACs), providing synergistic active sites and combining the benefits of SACs with bimetallic systems. Furthermore, this review systematically explores the latest advancements in synthesis methods and innovations for SACs for electrochemical water splitting. Additionally, it examines the evolution of catalyst design, emphasizing the unique structural and electronic characteristics of single-site, dual-site, and alloyed SAC systems and highlighting their critical roles in accelerating water-splitting reaction kinetics as well as the prevailing challenges and outlining promising directions for advancing hydrogen production via water electrolysis.

08 HYDROGEN↗

Coherent Acoustic Interferometry during the Photodriven Oxygen Evolution Reaction Associates Strain Fields with the Reactive Oxygen Intermediate (Ti–OH*)

The oxygen evolution reaction (OER) from water requires the formation of metastable, reactive oxygen intermediates to enable oxygen–oxygen bond formation. Conversely, such reactive intermediates could also structurally modify the catalyst. Here, a descriptor for the overall catalytic activity, the first electron and proton transfer OER intermediate from water, (M–OH*), has been associated with significant distortions of the metal–oxygen bonds upon charge-trapping. Time-resolved spectroscopy of in situ, photodriven OER on transition metal oxide surfaces has characterized M–OH* for the charge trapping and the symmetry of the lattice distortions by optical and vibrational transitions, respectively, but had yet to detect an interfacial strain field arising from a surface coverage M–OH*. Here, we utilize picosecond, coherent acoustic interferometry to detect the uniaxial strain normal to the SrTiO 3 /aqueous interface directly caused by Ti–OH*. The spectral analysis applies a fairly general methodology for detecting a combination of the spatial extent, magnitude, and generation time of the interfacial strain through the coherent oscillations’ phase. For lightly n-doped SrTiO 3 , we identify the strain generation time (1.31 ps), which occurs simultaneously with Ti–OH* formation, and a tensile strain of 0.06% (upper limit 0.6%). In addition to fully characterizing this intermediate across visible, mid-infrared, and now GHz-THz probes on SrTiO 3 , we show that strain fields occur with the creation of some M–OH*, which modifies design strategies for tuning catalytic activity and provides insight into photo-induced degradation so prevalent for OER. To that end, the work put forth here provides a unique methodology to characterize intermediate-induced interfacial strain across OER catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topological control of liquid-metal-dealloyed structures

Abstract The past few years have witnessed the rapid development of liquid metal dealloying to fabricate nano-/meso-scale porous and composite structures with ultra-high interfacial area for diverse materials applications. However, this method currently has two important limitations. First, it produces bicontinuous structures with high-genus topologies for a limited range of alloy compositions. Second, structures have a large ligament size due to substantial coarsening during dealloying at high temperature. Here we demonstrate computationally and experimentally that those limitations can be overcome by adding to the metallic melt an element that promotes high-genus topologies by limiting the leakage of the immiscible element during dealloying. We further interpret this finding by showing that bulk diffusive transport of the immiscible element in the liquid melt strongly influences the evolution of the solid fraction and topology of the structure during dealloying. The results shed light on fundamental differences in liquid metal and electrochemical dealloying and establish a new approach to produce liquid-metal-dealloyed structures with desired size and topologies.

36 MATERIALS SCIENCE↗

Evolution of Metastable Phases During Mg Metal Corrosion: An In Situ Cryogenic X-ray Photoelectron Spectroscopy Study

Magnesium and its alloys are potential structural materials candidates for a wide variety of applications due to their high strength-to-weight ratio. However, ductility and poor corrosion resistance under ambient environmental conditions are the bottleneck for industrial deployment. Designing passivation layers and/or corrosion resistant alloys requires fundamental understanding of the corrosion process. The traditional ex-situ spectroscopic measurements of polycrystalline metal surface with ubiquitous surface impurities and grain boundaries only provided an indistinct view of the corrosion process. To clearly distinguish the mechanism and sequence of corrosion process, we employed in-situ cryo-based x-ray photoelectron spectroscopy (XPS) measurements on Mg single crystal surface exposed to aqueous salt solution. Clean Mg (0001) surfaces were exposed to pure D2O and NaCl aqueous solution (5 wt% NaCl+95 wt% D2O). The interfacial reactions were studied using multimodal analysis including XPS, x-ray diffraction (XRD) and scanning electron microscopy (SEM). In contrast to previous studies, our experiments demonstrated the formation of magnesium chloride hydroxide hydrate during aqueous salt corrosion processes. Evidence of metastable ClO* radicals were also found during the initial aqueous salt solution exposure.

Shutthanandan, Vaithiyalingam↗

Reconstructing parent microstructures in martensitic and pearlitic Ti-Cu

Evaluating the microstructural evolution of parent phases has long been a challenge in metals with a partial or complete solid state transformation. Most parent microstructure evolution has to be inferred from the product microstructure, with some information permanently lost. Parent microstructure reconstructions relying on orientation relationships and EBSD data are a useful tool to overcome this challenge, though much background knowledge is required to do so. Here, this work introduces a reconstruction process for a eutectoid composition Ti-Cu binary alloy to evaluate the high temperature β-Ti phase from martensitic and pearlitic product microstructures. Martensitic microstructures were accurately reconstructed, enabling elucidation of β-Ti grain size and texture, and the employed script included for scientific reference. Pearlitic microstructures did not reconstruct effectively with the current process. However, new insight into the interfacial crystallographic orientations of pearlitic Ti-Cu microstructures was gleaned via this analysis, and future investigations of interest are discussed accordingly. Ultimately, these findings demonstrate tools widely implemented on conventional titanium and ferrous alloys can also be applied to understand parent microstructures in Ti-Cu and other alloys of emerging interest. This is especially of use as new manufacturing processes for such materials require understanding how parent phase microstructures respond to new, novel material processing such as AM.

36 MATERIALS SCIENCE↗

Interface instabilities in hafnium hydride entrained iron metal matrix composites

The chemical interactions in Fe–HfH 2 metal matrix composites (MMCs) are studied across multiple length scales to elucidate the decomposition of the parent phases and corresponding reaction zone physics during direct current sintering. Fe–HfH 2 composites were synthesized with increasing as-mixed hydride contents of Fe–25% HfH 2 , Fe–40% HfH 2 , Fe–55% HfH 2 , and Fe–70% HfH 2 (all in vol. %) to demonstrate the ability to achieve sintered MMCs with target hydride contents. Samples were probed across multiple length scales through a multi-modal workflow employing x-ray diffraction, scanning electron microscopy and segmentation analysis, and synchrotron techniques including hard x-ray fluorescence mapping and nanoprobe x-ray absorption near-edge structure measurements. Under the selected sintering temperature and pressure conditions, hydrogen evolution is seen to evolve through parallel paths: thermal decomposition from during the transformation of HfH 2 to HfH x<2 and through subsequent reaction with the Fe matrix leading to intermetallic phase formation. Specifically, HfFe and HfFe 2 intermetallic formation accelerates the release of hydrogen with a subsequent HfO 2 phase forming at grain boundaries. For this MMC, the consumption or loss of hydrogen can be considerable in compacts with initial hydride loading of 25%–40% HfH 2 approaching 83% hydrogen loss for the lower volume fraction composites. Increasing the volume fraction of HfH 2 to 70% enhanced the retained hydrogen content to 53% and attributed to the reduced interfacial area intrinsic to the increased HfH 2 loading in this MMC.

36 MATERIALS SCIENCE↗

Elucidating the Discharge Behavior of Aqueous Zinc Sulfur Batteries in the Presence of Molybdenum(IV) Chalcogenide Catalyst: The Criticality of Interfacial Electrochemistry

The aqueous zinc-sulfur battery holds promise for significant capacity and energy density with low cost and safe operation based on environmentally benign materials. However, it suffers from the sluggish kinetics of the conversion reaction. Here, we highlight the efficacy of molybdenum(IV) sulfide (MoS 2 ) to reduce the overpotential of S-ZnS conversion in aqueous electrolytes and study the discharge products formed at the solid-solid and solid-liquid interfaces using experimental and theoretical approaches. Specifically, the MoS 2 -catalyzed electrochemical conversion reaction is characterized via ex situ X-ray diffraction (XRD), transmission electron microscopy (TEM) with energy dispersive spectroscopy (EDS), Raman spectroscopy, synchrotron-based Mo K-edge X-ray absorption spectroscopy (XAS), and in situ synchrotron-based X-ray computed tomography (XCT). Additionally, operando synchrotron-based S K-edge XAS and X-ray fluorescence (XRF) maps are collected to determine the spatial evolution of sulfur-based species at the electrode-electrolyte interface. Further, coupling the operando S K-edge XAS data with the simulated spectra and fitting the data suggested a possible ZnS 2 intermediate phase.

25 ENERGY STORAGE↗

Molecular dynamics studies of simple membrane-water interfaces: Structure and functions in the beginnings of cellular life

Molecular dynamics computer simulations of the structure and functions of a simple membrane are performed in order to examine whether membranes provide an environment capable of promoting protobiological evolution. Our model membrane is composed of glycerol 1-monooleate. It is found that the bilayer surface fluctuates in time and space, occasionally creating thinning defects in the membrane. These defects are essential for passive transport of simple ions across membranes because they reduce the Born barrier to this process by approximately 40%. Negative ions are transferred across the bilayer more readily than positive ions due to favorable interactions with the electric field at the membrane-water interface. Passive transport of neutral molecules is, in general, more complex than predicted by the solubility-diffusion model. In particular, molecules which exhibit sufficient hydrophilicity and lipophilicity concentrate near membrane surfaces and experience 'interfacial resistance' to transport. The membrane-water interface forms an environment suitable for heterogeneous catalysis. Several possible mechanisms leading to an increase of reaction rates at the interface are discussed. We conclude that vesicles have many properties that make them very good candidates for earliest protocells. Some potentially fruitful directions of experimental and theoretical research on this subject are proposed.

Pohorille, Andrew↗