Search NASA⌕ Search

SEARCH · Search NASA

Results for “addition reactions”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Isotope effects at the origin of life: Fingerprints of the Strecker synthesis

We report strecker synthesis creates α-amino acids from prebiotically plausible substrates (cyanide, ammonia, and aldehydes) and is widely hypothesized to be a key mechanism in the chemistry that led to life on Earth and on other planets. To better constrain the synthetic environments and precursors of abiotic α-amino acids, and to determine unique signatures of abiogenic amino acids, we measured the molecular-averaged and site-specific carbon and nitrogen isotope effects for the Strecker synthesis of alanine. The reaction steps of the Strecker synthesis can be divided into two groups: an initial series of reversible amination and nitrile-addition reactions (‘equilibration’) and a second series of irreversible hydrolysis reactions (‘hydrolysis’). The equilibration of cyanide, acetaldehyde, and ammonia with the intermediate, α-aminopropionitrile (α-APN), has a measured 55.1‰ equilibrium nitrogen isotope effect between the 15 N-rich amine nitrogen in α-aminopropionitrile and the 15 N-poor ammonia and a 20.0‰ equilibrium carbon isotope effect between the 13 C-poor C-2 site in α-aminopropionitrile and the 13 C-rich carbonyl carbon in acetaldehyde. The first irreversible hydrolysis step is inferred to have an up to 10‰ normal carbon fractionation (i.e., faster for 12 C, slower for 13 C) for the whole molecule, but it also has one or more side reactions that deplete the reactive α-APN reservoir by up to 15‰. The second hydrolysis step has a 15.4‰ normal kinetic isotope effect on the amide (C-1) site of alaninamide, which becomes the carboxyl site of alanine. Other α-amino acids will likely experience similar nitrogen isotope fractionations between ammonia and their amine sites, and similar carbon isotope fractionations between the carbonyl carbon in reactant aldehydes or ketones and the intermediate α-aminonitrile, and between cyanide and the carboxyl site. Therefore, these isotope effects allow us to predict the carbon and nitrogen isotopic contents and intramolecular structures of α-amino acids formed by Strecker synthesis based on their substrates’ isotopic compositions, or to infer the isotopic compositions of substrates from which amino acids formed, for example in the case of the amino-acid-rich carbonaceous chondrites. The site-specific C and N isotopic compositions of amino acids formed by Strecker chemistry contrast with those typical of terrestrial biosynthetic amino acids, so these data also provide a means of discriminating between biogenic and abiogenic α-amino acids.

58 GEOSCIENCES↗
NASA NTRS2020

Stability of Actinolite on Venus

Venus currently has a hostile surface environment with temperatures of ~460 ºC, pres-sures near 92 bars, and an atmosphere composed of super critical CO2 hosting a myriad of other potentially reactive gases (e.g., SO2, HCl, HF). However, it has been proposed that its surface may not have always been so harsh. Models suggest there may have been billions of years of clement conditions allowing an Earth-like environment with liquid water oceans. If such conditions existed, it is possible Venus formed a similar array of hydrous or aqueous minerals as seen on other planets with liquid surface water (e.g., Mars, Earth). Based on thermodynamic modeling, many of these phases would not be stable under the current atmospheric conditions on Venus, dehydrating due to the high temperatures and low concentration of H2O in the atmosphere. However, the rate of decomposition of these phases may allow them to remain present on the surface over geologic time. For example, experiments on the reaction rate of tremolite (Ca2Mg5Si8O22(OH)2) show a 50% decomposition time of 2.7 Gyr for micrometer sized grains in unreactive atmospheres (i.e., without SO2) at 740 K, and a 50% decomposition time of 70 Gyr for crystals several millimeters to centimeters in size. If hydrous minerals can remain on the surface of Venus over geologic time, it has implications for our detection of evidence of these past environments, and also for the overall water budget of the planet. If after surficial dehydration the planet was able to still store water in its crust, possible processes such as subduction or metamorphism could still have operated using stored water long after liquid surface water evaporated. Several previous studies have focused on experimental investigations of mineral stability on Venus. In particular, the works of studied the decomposition rate of tremolite under conditions relevant to Venus. As their focus was on decomposition of the mineral due to lack of water in the atmosphere, their experiments were undertaken using only CO2 or N2 gas at atmospheric pressure. Re-cent experiments have examined reactivity of other minerals with the Venusian atmosphere using more complex gas compositions at similar pressures to those seen on Venus. These studies show reaction of silicate minerals with atmospheric components on relatively short timescales (i.e., on the order of days). The reported reactions of silicate materials in both studies produced iron oxides, Ca sulfates, and Na sulfates. These ions are present in many amphiboles, and Ca was proposed by Johnson and Fegley to potentially have an important role in the decomposition mechanism for tremolite, with the Ca-O bond being the first to break during decomposition. The potential involvement of Ca in both processes raises the question of whether or not the reaction to form a secondary mineral phase will influence the rate of amphibole break-down (e.g., discussion in for tremolite). Additionally, reaction of Ca with atmospheric gases may result in a different secondary mineral assemblage than simple amphibole decomposition, which will need to be recognized when searching for evidence of past hydrated minerals on the Venusian surface. In order to understand the effect of this reaction on the overall preservation potential of amphibole on the surface of Venus, we are conducting experiments in both reactive and nonreactive atmospheres using the mineral actinolite (Ca2(Mg,Fe)5Si8O22(OH)2), an amphibole with similar crystal structure to tremolite that contains both Ca and Fe.

Santos, A. R.↗
DOE OSTI2023

Sulfur ( 3 P) Reaction with Conjugated Dienes Gives Cyclization to Thiophenes under Single Collision Conditions

Here, we combine crossed-beam velocity map imaging with high-level ab initio/transition state theory modeling of the reaction of S( 3 P) with 1,3-butadiene and isoprene under single collision conditions. For the butadiene reaction, we detect both H and H 2 loss from the initial adduct, and from reaction with isoprene, we see both H loss and methyl loss. Theoretical calculations confirm these arise following intersystem crossing to the singlet surface forming long-lived intermediates. For the butadiene reaction, these lose H 2 to form thiophene as the dominant channel, H to form the detected 2H-thiophenyl radical, or ethene, giving thioketene. For isoprene, additional reaction products are suggested by theory, including the observed H and methyl loss radicals, but also methyl thiophene, thioformaldehyde, and thioketene. The results for S( 3 P) + 1,3-butadiene, showing direct cyclization to the aromatic product and yielding few bimolecular product channels, are in striking contrast to those for the analogous O( 3 P) reaction.

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

Data‐Driven Insights into Rare Earth Mineralization: Machine Learning Applications Using Functional Material Synthesis Data

Understanding rare‐earth element (REE) mineralization mechanisms is essential for developing efficient separation strategies. Although the geochemical pathways that generate REE deposits are qualitatively known, quantitative links between specific conditions and mineralization outcomes remain limited. Herein, the repurpose laboratory REE hydrothermal synthesis data—originally collected for functional‐materials fabrication—as a surrogate for studying mineralization with data‐driven methods. The compiled 1,200+ hydrothermal reaction records and trained three machine‐learning models—K‐nearest neighbors (KNN), random forest (RF), and extreme gradient boosting (XGB)—to predict product elements and phases from precursors, additives, reaction conditions, and engineered features. Validation shows XGB achieves the highest accuracy. Feature importance indicates thermodynamic properties of cations and anions dominate model decisions. Correlations reveal positive relationships among precursor concentration, reaction time, pH, and temperature, consistent with classical crystallization behavior. XGB‐based regressors are built to predict crystallization temperature and pH from precursor/product attributes. Performance is strongest when similar training examples exist, while accuracy declines for underrepresented reactions, notably REE carbonates and heavy‐REE systems. Overall, the study shows that functional‐materials datasets can illuminate REE mineralization and provide priors for exploration and processing. Expanding datasets with less‐studied chemistries and conditions will improve generality and support deposit discovery and more efficient REE recovery.

feature importance analysis↗
NASA NTRS1996

Detection of abundant ethane and methane, along with carbon monoxide and water, in comet C/1996 B2 Hyakutake: evidence for interstellar origin

The saturated hydrocarbons ethane (C2H6) and methane (CH4) along with carbon monoxide (CO) and water (H2O) were detected in comet C/1996 B2 Hyakutake with the use of high-resolution infrared spectroscopy at the NASA Infrared Telescope Facility on Mauna Kea, Hawaii. The inferred production rates of molecular gases from the icy, cometary nucleus (in molecules per second) are 6.4 X 10(26) for C2H6, 1.2 X 10(27) for CH4, 9.8 X 10(27) for CO, and 1.7 X 10(29) for H2O. An abundance of C2H6 comparable to that of CH4 implies that ices in C/1996 B2 Hyakutake did not originate in a thermochemically equilibrated region of the solar nebula. The abundances are consistent with a kinetically controlled production process, but production of C2H6 by gas-phase ion molecule reactions in the natal cloud core is energetically forbidden. The high C2H6/CH4 ratio is consistent with production of C2H6 in icy grain mantles in the natal cloud, either by photolysis of CH4-rich ice or by hydrogen-addition reactions to acetylene condensed from the gas phase.

Non-programmatic↗
DOE OSTI2023

Green's function knockout formalism

Knockout nuclear reactions, in which a nucleon is removed from a nucleus as a result of the collision with another nucleus, have been widely used as an experimental tool, both to populate isotopes further removed from stability and to obtain information about the single-particle nature of the nuclear spectrum. In order to fully exploit the experimental information, theory is needed for the description of both the structure of the nuclei involved, and the dynamics associated with the nucleon removal mechanisms. Here, the standard approach, using theoretical shell-model spectroscopic factors for the structure description coupled with an eikonal model of reaction, has been successful when used in the context of the removal of valence nucleons in nuclei close to stability. However, it has been argued that the reaction theory might need to be revisited in the case of exotic nuclei, more specifically for highly asymmetric nuclei in which the deficient species (neutrons or protons) is being removed. We present here a new formalism for the nucleon-removal and -addition reaction through knockout and transfer reactions, that treats consistently structure and reaction properties using dispersive optical potentials. In particular, our formalism includes the dynamical effects associated with the removal of a nucleon from the projectile, which might explain the long standing puzzle of the quenching of spectroscopic factors in nuclei with extreme neutrons-to-protons ratios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Method and system for storing and generating hydrogen

A method and system for storing and generating hydrogen. The method comprises generating hydrogen and heat from the reaction of a metal or metal compound with water. The heat generated from this reaction may then be converted to other forms of energy such as by passing the heat through a thermal electric device to recover electrical energy for storage in a battery. In an alternative and preferred embodiment, the heat is used to drive additional reactions for generating more hydrogen and is preferably used to drive an endothermic dehydrogenation reaction resulting in increased hydrogen generation and consumption of the heat.

Kindler, Andrew↗
DOE OSTI2021

A high-performance organic cathode customized for sulfide-based all-solid-state batteries

All-solid-state batteries (ASSBs) have become increasingly attractive recently due to their better safety and prospective long-term stability compared with conventional liquid batteries. However, obtaining a sustainable cathode candidate to match the solid electrolyte with regards to operating potential, chemical compatibility, and mechanical property is still an open challenge. In this work, the chemical incompatibility of quinone-based active materials and sulfide-based electrolyte were unveiled for the first time through a heteroconjugate addition reaction mechanism. To develop a quinone cathode customized for sulfide-based ASSBs, poly-(anthraquinonyl sulfide)-graphene (PAQS-G) nanocomposite was reported. The stable polymer framework of PAQS can protect the quinone redox center by preventing nucleophilic attack from sulfide-based solid electrolytes. The graphene additives can ameliorate redox kinetics and improve active material utilization. The PAQS-G cathode exhibited a specific capacity of ~178 mAh g -1 and a high material utilization of ~79%. Excellent cycling stability was achieved with 94 % capacity after 200 cycles in lithium batteries and 95.5 % capacity after 300 cycles in sodium batteries at 0.1C rate, respectively. A promising potential for energy storage applications was demonstrated.

25 ENERGY STORAGE↗
DOE OSTI2020

Plasmon-Mediated Intramolecular Methyl Migration with Nanoscale Spatial Control

Plasmonic materials interact strongly with light to focus and enhance electromagnetic radiation down to nanoscale volumes. Due to this localized confinement, materials that support localized surface plasmon resonances are capable of driving energetically unfavorable chemical reactions. In certain cases, the plasmonic nanostructures are able to preferentially catalyze the formation of specific photoproducts, which offers an opportunity for the development of solar-driven chemical synthesis. Here, using plasmonic environments, we report inducing an intramolecular methyl migration reaction, forming 4-methylpyridine from N-methylpyridinium. Using both experimental and computational methods, we were able to confirm the identity of the N-methylpyridinium by making spectral comparisons against possible photoproducts. This reaction involves breaking a C–N bond and forming a new C–C bond, highlighting the ability of plasmonic materials to drive complex and selective reactions. Additionally, we observe that the product yield depends strongly on optical illumination conditions. This is likely due to steric hindrance in specific regions on the nanostructured plasmonic substrate, providing an optical handle for driving plasmonic catalysis with spatial specificity. Furthermore, this work adds yet another class of reactions accessible by surface plasmon excitation to the ever-growing library of plasmon-mediated chemical reactions.

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

Cool flame chemistry of diesel surrogate compounds: n -Decane, 2-methylnonane, 2,7-dimethyloctane, and n -butylcyclohexane

Elucidating the formation of combustion intermediates is crucial to validate reaction pathways, develop reaction mechanisms and examine kinetic modeling predictions. While high-temperature pyrolysis and oxidation intermediates of alkanes have been thoroughly studied, comprehensive analysis of cool flame intermediates from alkane autoxidation is lacking and challenging due to the complexity of intermediate species produced. Here, jet-stirred reactor autoxidation of four C 10 alkanes: n -decane, 2-methylnonane, 2,7-dimethyloctane, and n -butylcyclohexane, as model compounds of diesel fuel, was investigated from 500 to 630 K using synchrotron vacuum ultraviolet photoionization molecular beam mass spectrometry (SVUV-PIMS). Around 100 intermediates were detected for each fuel. The classes of molecular structures present during the autoxidation of the representative paraffinic functional groups in transport fuels, i.e., n -alkanes, branched alkanes, and cycloalkanes were established and were found to be similar from the oxidation of various alkanes. A theoretical approach was applied to estimate the photoionization cross sections of the intermediates with the same carbon skeleton as the reactants, e.g., alkene, alkenyl keto, cyclic ether, dione, keto-hydroperoxide, diketo-hydroperoxide, and keto-dihydroperoxide. These species are indicators of the first, second, and third O 2 addition reactions for the four C 10 hydrocarbons, as well as bimolecular reactions involving keto-hydroperoxides. Chemical kinetic models for the oxidation of these four fuels were examined by comparison against mole fraction of the reactants and final products obtained in additional experiments using gas chromatography analysis, as well as the detailed species pool and mole fractions of aforementioned seven types of intermediates measured by SVUV-PIMS. This works reveals that the models in the literature need to be improved, not only the prediction of the fuel reactivity and final products, but also the reaction network to predict the formation of many previous undetected intermediates.

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

Single-Ion Conducting Polymer Electrolyte Enabled via Aza-Michael Addition

Single-ion conducting polymer electrolytes present a possible route for achieving high energy density next-generation batteries due to their flexible nature and high cation transport number. The trifluoromethanesulfonimide (TFSI) functional group represents one of the most efficient tethered anions for Li-ion conductivity. However, the covalent attachment of TFSI groups into polymer electrolytes has been challenging and costly due to its synthetic difficulty and limited commercial availability of building blocks. Here, we present a new polymer electrolyte synthesized by the Michael addition reaction between poly(allylamine) (PAA) and a vinyl TFSI anion under mild reaction conditions. The resulting PAATFSI exhibits a lower glass transition temperature and several orders of magnitude higher Li-ion conductivity than similar TFSI-based single-ion conducting homopolymers in the dry state. Moreover, PAATFSI exhibits excellent Li-ion conductivity (2.7 × 10 –4 S/cm at 30 °C) with the addition of a plasticizer (60 wt % of carbonate solvent), which enables stable lithium–metal battery performance. In conclusion, this study demonstrates a facile route for the synthesis of new single-ion conducting polymer electrolytes that opens a library of possibilities to enable the realization of polymer-based batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Activation of Lattice and Adatom Oxygen by Highly Stable Ceria-Supported Cu Single Atoms

Requiring catalysts to be both active yet stable over long periods of time under variable reaction conditions including high and low temperatures is a daunting challenge due to the almost mutual exclusivity of these constraints. Using CO oxidation as a probe reaction, in this work we demonstrate that thermally stable single atom copper catalysts prepared by high temperature synthesis (atom trapping) on ceria can achieve this feat by allowing modulation of the Cu charge state through facile charge transfer between active site and support. This provides the catalysts with an ability to activate either lattice or adatom oxygen atoms, accessing additional reaction channels depending on the reaction conditions. Such adaptability allows dynamic response of catalysts enabling them to remain active under variable reaction conditions. The inherent stability of the catalyst arises from the enhanced strength of the Cu-O interactions from high temperature synthesis that exist even when Cu oxidation state varies, which retards sintering and deactivation. As we show here, one can circumvent the dilemma of designing catalysts that are simultaneously active and stable by matching the redox properties of the active site and support and establishing an environmental adaptability around the active sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

A new synthetic route to a family of non-classical addition-type thermoplastics. I - Concept and demonstration

A synthesis scheme was developed for a host of nonclassical addition-type thermoplastics (ATT) that can be cured by an addition reaction, leading to a linear polymer structure. The synthesis involves the reaction of an acetylene-terminated prepolymer with either a bismaleimide (BMI) or a biscitraconimide. A new polymer, designated LaRC-RP80, synthesized using this scheme, was found to exhibit several significantly improved properties over the commercial BMI, Kerimid 601, including an eight-fold increase in toughness while maintaining a high Tg, a higher (by 167 C) heat stability, and a 50-percent increase in moisture resistance. In addition, LaRC-RP80 has good hot/wet lap shear strength and processes easily at 288 C without voids in the finished product.

Pater, Ruth H.↗
DOE OSTI2020

Quantitative assessment of environmental phenomena on maximum pit size predictions in marine environments

Maximum pit sizes were predicted for dilute and concentrated NaCl and MgCl 2 solutions as well as sea-salt brine solutions corresponding to 40% relative humidity (RH) (MgCl 2 -rich) and 76% RH (NaCl-rich) at 25 °C. A quantitative method was developed to capture the effects of various cathode evolution phenomena including precipitation and dehydration reactions. Additionally, the sensitivity of the model to input parameters was explored. Despite one's intuition, the highest chloride concentration (roughly 10.3 M Cl – ) did not produce the largest predicted pit size as the ohmic drop was more severe in concentrated MgCl 2 solutions. Therefore, the largest predicted pits were calculated for saturated NaCl (roughly 5 M Cl – ). Next, it was determined that pit size predictions are most sensitive to model input parameters for concentrated brines. However, when the effects of cathodic reactions on brine chemistry are considered, the sensitivity to input parameters is decreased. Although there was not one main input parameter that influenced pit size predictions, two main categories were identified. Under similar chloride concentrations (similar RH), the water layer thickness (WL), and pit stability product, (i · x) sf , are the most influential factors. When varying chloride concentrations (RH), changes in WL, the brine specific cathodic kinetics on the external surface (captured in the equivalent current density (i eq )), and conductivity (k o ) are the most influential parameters. Finally, it was noted that dehydration reactions coupled with precipitation in the cathode will have the largest effect on predicted pit size, and cause the most significant inhibition of corrosion damage.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗
DOE OSTI2026

Measurements of the inert Hugoniot and observation of mechanical ignition in Ni(V) + Al reactive multilayers via x-ray diffraction

This study investigates the inert Hugoniot response, mechanical ignition, and reaction dynamics of Ni(V)+Al multilayers during longitudinal, laser-driven shock compression experiments. Ni(V)+Al multilayers, known for their self-propagating exothermic reactions, were subjected to longitudinal stresses exceeding 50 GPa using the laser shock facility within the Dynamic Compression Sector (DCS) at the Advanced Photon Source (APS). In situ x-ray diffraction (XRD) revealed that Ni(V) and Al were not in equilibrium during compression, with stress discrepancies attributed to twinning, grain structure effects, and/or dislocation density. However, the measured inert Hugoniot closely matched prior experimental and computational studies, confirming the utility of XRD for measuring the equation of state of thin, complex materials. Additionally, reaction was observed at significantly higher stresses than reported previously using laser-launched flyers. This discrepancy suggests a strong influence of externally imposed shear stress on reaction thresholds, which likely arose from deviations in flyer planarity during past experiments. Full reaction of the multilayer occurred within 40 ns after shock-wave passage, evidenced by complete melting of the constituents. Eulerian hydrocode simulations replicated experimental conditions, providing insights into equilibrium dynamics and experimental artifacts. The results highlight how even small shear forces facilitate ignition in Ni(V)+Al multilayers at lower stresses.

Crystallography↗
NASA NTRS2001

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.

Kennedy, Christopher A.↗