Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed analytics”

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 1,297 records · Page 72

KiDS-1000 cosmology: Combined second- and third-order shear statistics

Aims.In this work, we perform the first cosmological parameter analysis of the fourth release of Kilo Degree Survey (KiDS-1000) data with second- and third-order shear statistics. This paper builds on a series of studies aimed at describing the roadmap to third-order shear statistics. Methods.We derived and tested a combined model of the second-order shear statistic, namely, the COSEBIs and the third-order aperture mass statistics 〈ℳ ap 3 〉 in a tomographic set-up. We validated our pipeline withN-body mock simulations of the KiDS-1000 data release. To model the second- and third-order statistics, we used the latest version of HMCODE2020 for the power spectrum and BIHALOFITfor the bispectrum. Furthermore, we used an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering only equal smoothing radii for the 〈ℳ ap 3 〉 part of the data vector. This makes it possible to carry out a data analysis of the KiDS-1000 data release using a combined analysis of COSEBIs and third-order shear statistics. Results.We first validated the accuracy of our modelling by analysing a noise-free mock data vector, assuming the KiDS-1000 error budget, finding a shift in the maximum of the posterior distribution of the matter density parameter, ΔΩ m < 0.02 σ Ω m , and of the structure growth parameter, ΔS 8 < 0.05 σ S 8 . Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained Ω m = 0.248 −0.055 +0.062 andS 8 = σ 8 √(Ω m /0.3 )= 0.772 ± 0.022. The geometric average on the errors of Ω m andS 8 of the combined statistics decreases, compared to the second-order statistic, by a factor of 2.2.

Astronomy & Astrophysics↗

A theory of localized excitations in supercooled liquids

Here, a new connection between the structure and dynamics in glass-forming liquids is presented. We show how the origin of spatially localized excitations, as defined by the dynamical facilitation (DF) theory, can be understood from a structure-based framework. This framework is constructed by associating excitation events in the DF theory to hopping events between energy minima in the potential energy landscape (PEL). By reducing the PEL to an equal energy well picture and applying a harmonic approximation, we develop a field theory to describe elastic fluctuations about inherent states, which are energy minimizing configurations of the PEL. We model an excitation as a shear transformation zone (STZ) inducing a localized pure shear deformation onto an inherent state. We connect STZs to T1 transition events that break the elastic bonds holding the local structure of an inherent state. A formula for the excitation energy barrier, denoted as J σ , is obtained as a function of inherent-state elastic moduli and the radial distribution function. The energy barrier from the current theory is compared to the one predicted by the DF theory where good agreement is found in various two-dimensional continuous poly-disperse atomistic models of glass formers. These results strengthen the role of structure and elasticity in driving glassy dynamics through the creation and relaxation of localized excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear elasticity and damping govern ultrafast dynamics in click beetles

Many small animals use springs and latches to overcome the mechanical power output limitations of their muscles. Click beetles use springs and latches to bend their bodies at the thoracic hinge and then unbend extremely quickly, resulting in a clicking motion. When unconstrained, this quick clicking motion results in a jump. While the jumping motion has been studied in depth, the physical mechanisms enabling fast unbending have not. Here, we first identify and quantify the phases of the clicking motion: latching, loading, and energy release. We detail the motion kinematics and investigate the governing dynamics (forces) of the energy release. We use high-speed synchrotron X-ray imaging to observe and analyze the motion of the hinge’s internal structures of four Elater abruptus specimens. We show evidence that soft cuticle in the hinge contributes to the spring mechanism through rapid recoil. Using spectral analysis and nonlinear system identification, we determine the equation of motion and model the beetle as a nonlinear single-degree-of-freedom oscillator. Quadratic damping and snap-through buckling are identified to be the dominant damping and elastic forces, respectively, driving the angular position during the energy release phase. The methods used in this study provide experimental and analytical guidelines for the analysis of extreme motion, starting from motion observation to identifying the forces causing the movement. Here, the tools demonstrated here can be applied to other organisms to enhance our understanding of the energy storage and release strategies small animals use to achieve extreme accelerations repeatedly.

59 BASIC BIOLOGICAL SCIENCES↗

Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born–Oppenheimer potential energy surface. In particular, we focus on ab initio–based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through ( a) parameterization of distributed multipoles, ( b) explicit fitting of the MBE, ( c) inclusion of many-atom features in a neural network, and ( d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear Fuel and Pu Redox Studies from The Glenn T. Seaborg Institute at Idaho National Laboratory

The Glenn T. Seaborg Institute at Idaho National Laboratory (INL-GTSI) focuses on advancing fundamental research in the actinide sciences by providing unique opportunities to early career scientists and engineers to gain experience studying the actinide elements and their associated systems. The INL-GTSI is built from three focus areas that are based on the expertise and supporting infrastructure at INL and include solid state chemistry and physics, solution phase chemistry and physics, and forensic and isotope science. INL is the lead Laboratory for nuclear energy research and development in the U. S. and the research on nuclear fuels performed under the INL-GTSI gives good examples of solid state studies. Uranium-Molybdenum (U-Mo) alloys are leading fuel candidates for conversion of high performance research and test reactors to low-enriched fuels. During irradiation, generated fission gas accumulates into bubbles and self-organizes into a gas bubble superlattice (GBS) that effectively stores fission gases and inhibits fuel swelling. A study on the early self- organizing behavior of the GBS shows that not only grain boundaries but the interfaces between the U-Mo matrix and uranium carbide (UC) impurities are important to GBS formation.[1] In solution, understanding the complex redox behavior of plutonium in aqueous environments is critical for establishing optimized nuclear waste reprocessing solvent systems and storage tank environments. INL-GTSI researchers have produced an experimentally validated multi-scale model of the gamma radiation induced behavior of plutonium ions in concentrated aqueous HNO3 solutions.[2] Here, gamma radiation effected only minimal steady state changes in the redox distribution of the plutonium oxidation states. The redox cycling between Pu(IV) and Pu(III) is demonstrated to be mediated by the •OH/NO3• radical oxidation of Pu(III) and the H2O2/HNO3 driven reduction of Pu(IV). The INL-GTSI offers young researchers the unique chance to work directly with actinide bearing materials in a U. S. National Laboratory environment. Further topical areas of interest to the INL-GTSI include, but are not limited to, fundamental actinide properties, structure/property (electronic, magnetic, thermal) relations, actinide quantum criticality, f- electron interactions, electron correlations, computational studies, new phases, defect effects, interface interactions, isotope production and separation, forensic analytical chemistry, structure and dynamic properties of actinides in non-aqueous media, separations chemistry and kinetics for advanced nuclear fuel cycles, radiation effects, and innovative and advanced ligand design for complexation of the actinides.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Metastable Metallic Phase of a Bilayer Blue Phosphorene Induced by Interlayer Bonding and Intralayer Charge Redistributions

We have carried out diffusion Monte Carlo calculations for an A 1 B -1 -stacked bilayer blue phosphorene to find that it undergoes a semiconductor-metal transition as the interlayer distance decreases. While the most stable bilayer structure is a semiconducting one with two monolayers coupled through weak van der Waals interaction, the metallic bilayer at a shorter interlayer distance is found to be only metastable. This is in contrast to a recent theoretical prediction based on random phase approximation that the metallic phase would be the most stable bilayer configuration of blue phosphorene. Our analysis of charge density distributions reveals that the metastable metallic phase is induced by interlayer chemical bonding and intralayer charge redistributions. This study enriches our understanding of interlayer binding of a blue phosphorene and contributes to establishment of correct energetic order between its different phases, which will be essential in devising an experimental pathway for a metallic phosphorene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis of End-Functionalized Poly(norbornenes) and Poly(ethylidene norbornenes) Using a Pd(II) Catalyst in Combination with Chain Transfer Agents

Synthesis of end-functionalized poly(norbornene)s (PNBs) and poly(ethylidenenorbornene)s (PENBs) were prepared using the palladium catalyst t Bu 3 PPdMeCl/NaBAr f in combination with chain transfer agents (CTAs) 4-penten-1-ol and 10-undecenal. These CTAs produce polymers end-capped with -OH and -CHO groups, respectively, in high yields with M n values in the range of 6K to 17K, depending on the molar ratio of the CTA: monomer. End-functionalized PNBs and PENBs with precisely controlled Mn values and very narrow molecular weight distributions were accessed by using 4-penten-1-ol and 10-undecenal as chain-terminating agents following living polymerizations in which monomer is completely consumed. The utility of the aldehyde-functionalized polymers is illustrated via coupling to a dye molecule and by synthesis of ambiphilic diblock copolymers, PNB-b-PEG and PENB-b-PEG.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Covalently integrated silica nanoparticles in poly(ethylene glycol)-based acrylate resins: thermomechanical, swelling, and morphological behavior

Nanocomposites integrate functional nanofillers into viscoelastic matrices for electronics, lightweight structural materials, and tissue engineering. Herein, the effect of methacrylate-functionalized (MA-SiO 2 ) and vinyl-functionalized (V-SiO 2 ) silica nanoparticles on the thermal, mechanical, physical, and morphological characteristics of poly(ethylene glycol) (PEG) nanocomposites was investigated. The gel fraction of V-SiO 2 composites decreases upon addition of 3.8 wt% but increases with further addition (>7.4 wt%) until it reaches a plateau at 10.7 wt%. The MA-SiO 2 induced no significant changes in gel fraction and both V-SiO 2 and MA-SiO 2 nanoparticles had a negligible impact on the nanocomposite glass transition temperature and water absorption. The Young's modulus and ultimate compressive stress increased with increasing nanoparticle concentration for both nanoparticles. Due to the higher crosslink density, MA-SiO 2 composites reached a maximum mechanical stress at a concentration of 7.4 wt%, while V-SiO 2 composites reached a maximum at a concentration of 10.7 wt%. Scanning electron microscopy, transmission electron microscopy, and small-angle X-ray scattering revealed a bimodal size distribution for V-SiO 2 and a monomodal size distribution for MA-SiO 2 . Although aggregates were observed for both nanoparticle surface treatments, V-SiO 2 dispersion was poor while MA-SiO 2 were generally well-dispersed. Finally, these findings lay the framework for silica nanofillers in PEG-based nanocomposites for advanced manufacturing applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization and Valuation of the Uncertainty of Calibrated Parameters in Microsimulation Decision Models

We evaluated the implications of different approaches to characterize the uncertainty of calibrated parameters of microsimulation decision models (DMs) and quantified the value of such uncertainty in decision making. We calibrated the natural history model of CRC to simulated epidemiological data with different degrees of uncertainty and obtained the joint posterior distribution of the parameters using a Bayesian approach. We conducted a probabilistic sensitivity analysis (PSA) on all the model parameters with different characterizations of the uncertainty of the calibrated parameters. We estimated the value of uncertainty of the various characterizations with a value of information analysis. We conducted all analyses using high-performance computing resources running the Extreme-scale Model Exploration with Swift (EMEWS) framework. The posterior distribution had a high correlation among some parameters. The parameters of the Weibull hazard function for the age of onset of adenomas had the highest posterior correlation of -0.958. When comparing full posterior distributions and the maximum-a-posteriori estimate of the calibrated parameters, there is little difference in the spread of the distribution of the CEA outcomes with a similar expected value of perfect information (EVPI) of $\$$653 and $\$$685, respectively, at a willingness-to-pay (WTP) threshold of $\$$66,000 per quality-adjusted life year (QALY). Ignoring correlation on the calibrated parameters’ posterior distribution produced the broadest distribution of CEA outcomes and the highest EVPI of $\$$809 at the same WTP threshold. Different characterizations of the uncertainty of calibrated parameters affect the expected value of eliminating parametric uncertainty on the CEA. Ignoring inherent correlation among calibrated parameters on a PSA overestimates the value of uncertainty.

97 MATHEMATICS AND COMPUTING↗

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A life cycle analysis (LCA) primer for the agricultural community

Consumer demand for sustainably produced agricultural products is increasing. One approach that can be used to evaluate and compare product sustainability is to conduct a life cycle analysis (LCA). The term LCA is broadly used to describe a suite of analytical resources and standardized methods. Life cycle analyses are a transdisciplinary tool which can and are used by all professions but work best when conducted by teams with diverse skills and backgrounds. The objective of this paper is to provide a primer on LCAs for the agricultural community and for LCA practitioners unfamiliar with agronomy. Agricultural LCAs can differ greatly because agricultural products have multiple end uses, complex socio-economic and environmental trade-offs, and can be generated using a variety of different practices, resources, and production systems. Through worldwide agricultural marketing, intricate agronomic supply chains and support systems have been developed to produce, store, and distribute agricultural products. In many situations, the production practices are not linked to the products consumers purchase. This limitation can be partially overcome by conducting an LCA, which is a systematic summation of the multiple environmental impacts associated with the production, use, and disposal of a product. However, the large number and types of analyses can make it difficult for non-practitioners and decision-makers to use and interpret LCA results. Finally, it is important for members of the agricultural community to be aware of these analyses and how they can be used and how they affect product markets.

54 ENVIRONMENTAL SCIENCES↗

Understanding Relationships between Free Volume and Oxygen Absorption in Ionic Liquids

In this report understanding the factors that govern gas absorption in ionic liquids is critical to the development of high-capacity solvents for catalysis, electrochemistry, and gas separations. Here, we report experimental probes of liquid structure that provide insights into how free volume impacts the O 2 absorption properties of ionic liquids. Specifically, we establish that isothermal compressibility - measured rapidly and accurately through small-angle X-ray scattering - reports on the size distribution of transient voids within a representative series of ionic liquids and is correlated with O 2 absorption capacity. Additionally, O 2 absorption capacities are correlated with thermal expansion coefficients, reflecting the beneficial effect of weak intermolecular interactions in ionic liquids on free volume and gas absorption capacity. Molecular dynamics simulations show that the void size distribution - in particular, the probability of forming larger voids within an ionic liquid - has a greater impact on O 2 absorption than the total free volume. These results establish relationships between the ionic liquid structure and gas absorption properties that offer design strategies for ionic liquids with high gas solubilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formation of Disordered Cocontinuous Phases by Randomly Linked Star Copolymers

Cocontinuous polymeric nanostructures have garnered significant interest due to their ability to combine different properties of two separate polymer domains. Randomly linked copolymer networks have proven to be especially robust for formation of disordered cocontinuous phases across wide composition ranges (≈30 wt % or more). While theoretical treatments of microphase-separated networks have focused primarily on the role of random elastic forces imposed on the self-assembled nanostructures by virtue of the network architecture, experimental studies seeking to disentangle these contributions from other potential effects, such as dispersity in preferred interfacial curvatures, have been scarce. To provide insight into this matter, we here study the self-assembly of randomly linked star copolymers (RSCs), constructed by linking premade polymer arms of polystyrene (PS) and poly(d,l-lactide) (PLA) using 3, 4, and 6-functional connectors. This architecture yields similar distributions of preferred curvature as networks made using corresponding difunctional strands, but lacks the elastic forces imposed by a network architecture. Gravimetry and small-angle X-ray scattering, coupled with scanning electron microscopy, were performed to identify the percolation of PS/PLA RSCs. Remarkably, the 4-arm RSC system exhibited a disordered cocontinuous window of ≈25 wt %, indicating that dispersity in preferred curvature can in some cases be sufficient to robustly drive formation of this morphology. However, the other RSC architectures showed smaller cocontinuous ranges, which we interpret in terms of the influence of homopolymer stars in the 3-arm case and the narrower distribution of preferred interfacial curvatures in the 6-arm case. Finally, thin layers of interconnected porous PS were achieved by solution-processing, suggesting that RSCs have the potential to serve as a robust and easily processable cocontinuous polymeric nanomaterials in both bulk and membrane geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhancing Low-Temperature Syngas Production via Surface Tailoring of Supported Intermetallic Nanocatalysts

An active and coke-resistant silica-encapsulated intermetallic Ni 3 Zn nanoparticle catalyst was developed for low-temperature (450 °C) dry reforming of methane (DRM). The catalyst exhibited a remarkable 4-fold increase in activity (4.5 s –1 ) with over 99% CO selectivity and 3 orders of magnitude less carbonaceous species and demonstrated remarkable stability (70 h) compared to that of a monometallic Ni catalyst. Here, the key is the combined effect of surface ensemble structure and electronic interaction modulation through the surface composition tailoring achieved by off-stoichiometric Ni and Zn loading in controlling surface chemistry for achieving different activities and H 2 /CO ratios. Characterized by ion spectroscopy, X-ray photoelectron spectroscopy, and the neutron pair distribution function, it was revealed that paired Ni δ – –Zn δ+ active sites are crucial for DRM. Transient infrared spectroscopy and isotopic analysis uncovered the synergistic effect of Ni δ – –Zn δ+ sites in activating C–H bonds and dissociating CO 2 to prevent coke formation under low-temperature conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modulating Microphase Separation of Lamellae-Forming Diblock Copolymers via Ionic Junctions

In this work, we present a molecular dynamics simulation study investigating the phase behavior of lamellae-forming diblock copolymers with a single ionic junction on the backbone. Our results show qualitative agreement with experimental findings regarding enhanced microphase separation with the introduction of an ionic junction at the conjunction point, while further revealing nonmonotonic changes in domain spacing and order–disorder transition as a function of the electrostatic interaction strength. This highlights the dominant roles of entropic and binding effects of counterions under weak and strong ionic correlations, respectively. The location of the ionic junction is found to effectively modulate the charge distribution and chain conformation in the ordered domains; its presence in the middle of a block promotes folding of the block, leading to a smaller domain size. These findings demonstrate the interplay of ionic coupling with steric hindrance and chain end effects, which enhances our understanding of the delicate control over the microphase domain features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Butene-Rich Alkene Formation from 2,3-Butanediol through Dioxolane Intermediates

The cost-effective production of sustainable aviation fuels (SAF) remains a major challenge within the energy sector. One approach to address this is the fermentation of biomass feedstocks into oxygenates followed by catalytic conversion to alkenes or other oligomerization precursors. 2,3-Butanediol (BDO) is a promising fermentation product due to its four-carbon nature, its decreased microorganism toxicity and associated higher maximum fermentation titers relative to other alcohols and oxygenates, and its capacity to be readily converted into butene isomers and longer chain alkenes. BDO conversion is currently constrained by separation challenges for BDO isolation due to its high boiling point and hydrophilicity. Here, this work expands upon previous BDO reactive separation via dioxolane formation over a solid acid catalyst by investigating the conversion of dioxolanes into alkene mixtures. Dioxolanes were formed from a range of aldehydes and subsequently converted over a Cu/ZSM-5 catalyst (448–523 K) via an ether cleavage, hydrogenation, and dehydration reaction network to form alkene-rich product mixtures (96% C 3+ alkene yield, 523 K). This selectivity is greater than that of direct BDO conversion to alkenes over an identical catalyst (89%, 523 K). C 3+ alkene selectivity is maximized between 498 and 523 K at complete dioxolane conversion without significant alkene hydrogenation to alkanes. The alkene product distributions can be tailored via both aldehyde selection during dioxolane formation and the dioxolane conversion reaction temperature. Alkene mixtures from dioxolane conversion predominantly reflect the carbon chain length and stereochemistry of BDO and the initial aldehyde at or below 498 K, yet higher reaction temperatures yield alkene mixtures of similar carbon chain distributions, regardless of initial aldehyde selection. Deactivation of the Cu/ZSM-5 catalyst is observed for multiple steps of the overall reaction network but can be minimized by facilitating the complete dioxolane-to-alkene reaction network at temperatures of at least 498 K.

2,3-butanediol↗

Diaminopurine in Nonenzymatic RNA Template Copying

In the RNA World before the emergence of an RNA polymerase, nonenzymatic template copying would have been essential for the transmission of genetic information. However, the products of chemical copying with the canonical nucleotides (A, U, C, and G) are heavily biased toward the incorporation of G and C, which form a more stable base pair than A and U. We therefore asked whether replacing adenine (A) with diaminopurine (D) might lead to more efficient and less biased nonenzymatic template copying by making a stronger version of the A:U pair. As expected, primer extension substrates containing D bound to U in the template more tightly than substrates containing A. However, primer extension with D exhibited elevated reaction rates on a C template, leading to concerns about fidelity. Our crystallographic studies revealed the nature of the D:C mismatch by showing that D can form a wobble-type base pair with C. We then asked whether competition with G would decrease the mismatched primer extension. We performed nonenzymatic primer extension with all four activated nucleotides on randomized RNA templates containing all four letters and used deep sequencing to analyze the products. We found that the DUCG genetic system exhibited a more even product distribution and a lower mismatch frequency than the canonical AUCG system. Furthermore, primer extension is greatly reduced following all mismatches, including the D:C mismatch. Our study suggests that D deserves further attention for its possible role in the RNA World and as a potentially useful component of artificial nonenzymatic RNA replication systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗