Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular”

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 685 records · Page 38

Polarized Resonant X-ray Scattering to Measure Molecular Orientation and Conformation in Organic Nanostructures

The ordering of organic (carbon-based) molecules into nanostructures drives useful materials properties for applications such as organic light emitting diodes, flexible batteries, and polymer nanocarriers. However, there is a current lack of techniques to measure this ordering. X-ray scattering techniques involving molecular resonances have shown promise in identifying molecular types, locations, and orientation. However, resonant soft X-ray scattering (RSoXS) techniques have not been developed even lacking basic optical models. The objectives achieved by this research were the development of methods, models, and analysis techniques that harness these resonant phenomena in X-ray scattering measurements to quantitatively probe ordering of organic molecules within devices and nanostructures. This award resulted in 11 peer reviewed publications with two others soon to be submitted. New physical models for X-ray scattering that accurately describe these complex molecular resonances enable measurements of how molecules assemble, orient and conform at the nanometer scale. These new techniques were applied to both organic electronic devices and aqueous nanoparticles being developed for hydrocarbon sequestration. This new measurement capability will enable characterization of structure-property relationships that lead to disruptive technologies of printable, stretchable, and biocompatible devices based on organic materials.

14 SOLAR ENERGY↗

Determination of Vibrational Motions Driving Photoinduced Electron Transfer Reactions in Molecular Crystals and Organic Thin Films

The objective of this research is to develop an understanding of how nuclear motions can be used to drive and control photo-induced charge transfer and singlet fission reactions in molecular crystals and thin films. By following the structural evolution of reacting molecules and nanostructures on the timescale of their nuclear motion, we will determine what interplay of nuclear coordinates is most efficient in driving electron transfers, a necessary first step in most photovoltaic and photocatalytic processes. We will identify which nuclear coordinates effectively promote desirable electron transfer processes, or ineffectively dissipate photon energy into unwanted thermal energy, thus guiding the design and implementation of molecular systems tailored for high efficiency processes. Specifically, we aim to identify which nuclear coordinates are responsible for driving electron transfer and singlet fission reactions at different points along the reaction coordinate, particularly in regimes which deviate from Marcus theory predictions. We will do so by investigating systems with electron transfer and singlet fission occurring from non-thermalized states due to strong donor-acceptor electronic coupling. First, we will prove that specific vibrational modes are responsible for charge transfer processes in molecular crystals. Secondly, we will investigate how vibrational coherences and specific nuclear motions can be used to drive high yield singlet fission processes. Finally, we will determine how complex film morphology affects long-range charge transport, a key limiting factor in our ability to make gains in efficiency.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bottom-up design of actinide materials from molecular clusters: Demonstration of a general-purpose simulation capability leveraging machine-learned atomic potentials

Actinide thin-film coatings such as uranium dioxide (UO 2 ) play an important role in nuclear reactors and other mission-relevant applications, but realization of their potential requires a deep fundamental understanding of the chemical vapor deposition (CVD) processes used for their growth. The slow experimental progress can be attributed, in part, to the standard safety guidelines associated with handling uranium byproducts, which are often corrosive, toxic, and radioactive. Accurate simulation techniques, when used in concert with experiment, can improve laboratory safety, material durability, and deliverable timeframes. However, state-of-the-art computational methods are either insufficiently accurate or intractably expensive. To remedy this situation, in this project we suggested a machine-learning (ML) accelerated workflow for simulating molecular clustering toward deposition. As a benchmark test case, we considered molecular clustering in steam and assessed independent components of our workflow by comparing with measured thermodynamic properties of water. After analyzing each component individually and finding no fundamental barrier to realization of the workflow, we attempted to integrate the ML component, a Sandia-developed tool called FitSNAP. As this was the first application of FitSNAP to atoms and molecules in the gas phase at Sandia, the method required more fitting data than was originally anticipated. Systematic improvements were made by including in the fit data diatomic potentials, molecular single-bond-breaking curves, and symmetry-constrained intermolecular potentials. We concluded that our strategy provides a feasible pathway toward modeling CVD and related processes, but that extensive training data must be generated before it can be of practical use.

36 MATERIALS SCIENCE↗

DataSet for Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR

Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic-inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. To understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31-calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ~2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31-calcite interactions. These studies provide the first experimental evidence to confirm that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template. However, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31-calcite interface has advanced our fundamental understanding of the atomic interactions at the organic-inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.

Saccuzzo Close, Emily Grace [Pacific Northwest Nat↗

Heterospecific Neighbor Plants Impact Root Microbiome Diversity and Molecular Function of Root Fungi

Within the forest community, competition and facilitation between adjacent-growing conspecific and heterospecific plants are mediated by interactions involving common mycorrhizal networks. The ability of plants to alter their neighbor’s microbiome is well documented, but the molecular biology of plant-fungal interactions during competition and facilitation has not been previously examined. We used a common soil-plant bioassay experiment to study molecular plant-microbial interactions among rhizosphere communities associated with Pinus taeda (native host) and Populus trichocarpa (non-native host). Gene expression of interacting fungal and bacterial rhizosphere communities was compared among three plant-pairs: Populus growing with Populus , Populus with Pinus , and Pinus with Pinus . Our results demonstrate that heterospecific plant partners affect the assembly of root microbiomes, including the changes in the structure of host specific community. Comparative metatranscriptomics reveals that several species of ectomycorrhizal fungi (EMF) and saprotrophic fungi exhibit different patterns of functional and regulatory gene expression with these two plant hosts. Heterospecific plants affect the transcriptional expression pattern of EMF host-specialists (e.g., Pinus -associated Suillus spp.) on both plant species, mainly including the genes involved in the transportation of amino acids, carbohydrates, and inorganic ions. Alteration of root microbiome by neighboring plants may help regulate basic plant physiological processes via modulation of molecular functions in the root microbiome.

59 BASIC BIOLOGICAL SCIENCES↗

Competing Easy-Axis Anisotropies Impacting Magnetic Tunnel Junction-Based Molecular Spintronics Devices (MTJMSDs)

Molecular spintronics devices (MSDs) attempt to harness molecules’ quantum state, size, and configurable attributes for application in computer devices—a quest that began more than 70 years ago. In the vast number of theoretical studies and limited experimental attempts, MSDs have been found to be suitable for application in memory devices and futuristic quantum computers. MSDs have recently also exhibited intriguing spin photovoltaic-like phenomena, signaling their potential application in cost-effective and novel solar cell technologies. The molecular spintronics field’s major challenge is the lack of mass-fabrication methods producing robust magnetic molecule connections with magnetic electrodes of different anisotropies. Another main challenge is the limitations of conventional theoretical methods for understanding experimental results and designing new devices. Magnetic tunnel junction-based molecular spintronics devices (MTJMSDs) are designed by covalently connecting paramagnetic molecules across an insulating tunneling barrier. The insulating tunneling barrier serves as a mechanical spacer between two ferromagnetic (FM) electrodes of tailorable magnetic anisotropies to allow molecules to undergo many intriguing phenomena. Our experimental studies showed that the paramagnetic molecules could produce strong antiferromagnetic coupling between two FM electrodes, leading to a dramatic large-scale impact on the magnetic electrode itself. Recently, we showed that the Monte Carlo Simulation (MCS) was effective in providing plausible insights into the observation of unusual magnetic domains based on the role of single easy-axis magnetic anisotropy. Here, we experimentally show that the response of a paramagnetic molecule is dramatically different when connected to FM electrodes of different easy-axis anisotropies. Motivated by our experimental studies, here, we report on an MCS study investigating the impact of the simultaneous presence of two easy-axis anisotropies on MTJMSD equilibrium properties. In-plane easy-axis anisotropy produced multiple magnetic phases of opposite spins. The multiple magnetic phases vanished at higher thermal energy, but the MTJMSD still maintained a higher magnetic moment because of anisotropy. The out-of-plane easy-axis anisotropy caused a dominant magnetic phase in the FM electrode rather than multiple magnetic phases. The simultaneous application of equal-magnitude in-plane and out-of-plane easy-axis anisotropies on the same electrode negated the anisotropy effect. Our experimental and MCS study provides insights for designing and understanding new spintronics-based devices.

42 ENGINEERING↗

Design, Synthesis, and Utility of Defined Molecular Scaffolds

A growing theme in chemistry is the joining of multiple organic molecular building blocks to create functional molecules. Diverse derivatizable structures—here termed “scaffolds” comprised of “hubs”—provide the foundation for systematic covalent organization of a rich variety of building blocks. This review encompasses 30 tri- or tetra-armed molecular hubs (e.g., triazine, lysine, arenes, dyes) that are used directly or in combination to give linear, cyclic, or branched scaffolds. Each scaffold is categorized by graph theory into one of 31 trees to express the molecular connectivity and overall architecture. Rational chemistry with exacting numbers of derivatizable sites is emphasized. The incorporation of water-solubilization motifs, robust or self-immolative linkers, enzymatically cleavable groups and functional appendages affords immense (and often late-stage) diversification of the scaffolds. Altogether, 107 target molecules are reviewed along with 19 syntheses to illustrate the distinctive chemistries for creating and derivatizing scaffolds. The review covers the history of the field up through 2020, briefly touching on statistically derivatized carriers employed in immunology as counterpoints to the rationally assembled and derivatized scaffolds here, although most citations are from the past two decades. The scaffolds are used widely in fields ranging from pure chemistry to artificial photosynthesis and biomedical sciences.

Sato, Daisuke (ORCID:0000000305472659)↗

Resolved Molecular Gas Observations of MaNGA Post-starbursts Reveal a Tumultuous Past

Post-starburst (PSB) galaxies have recently and rapidly quenched their star formation; thus, they are an important way to understand how galaxies transition from star-forming late types to quiescent early types. The recent discovery of large cold gas reservoirs in PSB galaxies calls into question the theory that galaxies must lose their gas to become quiescent. Optical Integral Field Spectroscopy (IFS) surveys have revealed two classes of PSB galaxies: central PSB (cPSB) galaxies with central quenching regions and ring PSB (rPSB) galaxies with quenching in their outskirts. We analyze a sample of 13 nearby (z < 0.1) PSB galaxies with spatially resolved optical IFS data from the Mapping Nearby Galaxies at Apache Point Observatory survey and matched resolution Atacama Large Millimeter/submillimeter Array observations of 12 CO(1–0). Disturbed stellar kinematics in 7/13 of our PSB galaxies and centrally concentrated molecular gas is consistent with a recent merger for most of our sample. In galaxies without merger evidence, alternate processes may funnel gas inward and suppress star formation, which may include outflows, stellar bars, and minor mergers or interactions. The star formation efficiencies of the PSB regions in nearly half our galaxies are suppressed while the gas fractions are consistent with star-forming galaxies. Active galactic nucleus (AGN) feedback may drive this stabilization, and we observe AGN-consistent emission in the centers of 5/13 galaxies. Finally, our cPSB and rPSB galaxies have similar properties except the ionized and molecular gas in cPSB galaxies is more disturbed. Overall, the molecular gas in our PSB galaxies tends to be compact and highly disturbed, resulting in concentrated gas reservoirs unable to form stars efficiently.

79 ASTRONOMY AND ASTROPHYSICS↗

Rapid Growth of Galactic Supermassive Black Holes through Accreting Giant Molecular Clouds during Major Mergers of Their Host Galaxies

Understanding the formation of the supermassive black holes (SMBHs) present in the centers of galaxies is a crucial topic in modern astrophysics. Observations have detected SMBHs with masses of 10 9 M ⊙ in the high-redshift galaxies with z ~ 7. However, how SMBHs grew to such huge masses within the first billion years after the Big Bang remains elusive. One possible explanation is that SMBHs grow quickly through the frequent mergers of galaxies, which provides sustainable gas to maintain rapid growth. This study presents the hydrodynamics simulations of the SMBHs’ growth with their host galaxies using the GIZMO code. In contrast to previous simulations, we have developed a giant molecular cloud (GMC) model by separating molecular gas particles from atomic gas particles and then evolving them independently. During major mergers, we show that the more massive molecular gas particles’ cloud bears stronger dynamical friction. Consequently, GMCs are substantially accreted onto the galactic centers that grow SMBHs from ~10 7 to ~10 9 M ⊙ within 300 Myr, which explains the rapid growth of SMBHs, and this accretion also triggers a violent starburst at the galactic center. Furthermore, we examine the impact of minor mergers on the bulge of a Milky Way–like galaxy and find that the size and mass of the bulge can increase from 0.92 to 1.9 kpc and from 4.7 × 10 10 to 7 × 10 10 M ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗

Ionic transport in electrostatic Janus Membranes. An explicit solvent molecular dynamic simulation.

Janus –or two-sided, charged membranes offer promise as ionic current rectifiers. In such systems, pores consisting of two regions of opposite charge can be used to generate a current from a gradient in salinity. The efficiency of nanoscale Janus pores increases dramatically as their diameter becomes smaller. However, little is known about the underlying transport processes, particularly under experimentally accessible conditions. In this work, we examine the molecular basis for rectification in Janus nanopores using an applied electric field. Molecular simulations with explicit water and ions are used to examine the structure and dynamics of all molecular species in aqueous electrolyte solutions. For several macroscopic observables, the results of such simulations are consistent with experimental observations on asymmetric membranes. Our analysis reveals a number of previously unknown features, including a pronounced local re-orientation of water molecules in the pores, and a segregation of ionic species that has not been anticipated by previously reported continuum analyses of Janus pores. Using these insights, a model is proposed for ionic current rectification in which electric leakage at pore entrance controls net transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HAMscope: a snapshot Hyperspectral Autofluorescence Miniscope for real-time molecular imaging

We introduce HAMscope, a compact, snapshot hyperspectral autofluorescence miniscope that enables real-time, label-free molecular imaging in a wide range of biological systems. By integrating a thin polymer diffuser into a widefield miniscope, HAMscope spectrally encodes each frame and employs a probabilistic deep learning framework to reconstruct 30-channel hyperspectral stacks (452-703 nm) or directly infer molecular composition maps from single images. A scalable multi-pass U-Net architecture with transformer-based attention and per pixel uncertainty estimation enables high spatio-spectral fidelity (mean absolute error ∼0.0048) at video rates. While initially demonstrated in plant systems, including lignin, chlorophyll, and suberin imaging in intact poplar and cork tissues, the platform is readily adaptable to other applications such as neural activity mapping, metabolic profiling, and histopathology. We show that the system generalizes to out-of-distribution tissue types and supports direct molecular mapping without the need for spectral unmixing. HAMscope establishes a general framework for compact, uncertainty-aware spectral imaging that combines minimal optics with advanced deep learning, offering broad utility for real-time biochemical imaging across neuroscience, environmental monitoring, and biomedicine.

59 BASIC BIOLOGICAL SCIENCES↗

Understanding the Reaction of Nuclear Graphite with Molecular Oxygen

A thorough understanding of oxidation is important when considering the health and integrity of graphite components in graphite reactors. For the next generation of graphite reactors, HTGRs specifically, an unlikely air ingress has been deemed significant enough to have made its way into the licensing applications of many international licensing bodies. While a substantial body of literature exists on nuclear graphite oxidation in the presence of molecular oxygen and significant efforts have been made to characterize oxidation kinetics of various grades, the value of existing information is somewhat limited. Often, multiple competing processes, including reaction kinetics, mass transfer, and microstructural evolution, are lumped together into a single rate expression that limits the ability to translate this information to different conditions. This article reviews the reaction of graphite with molecular oxygen in terms of the reaction kinetics, gas transport, and microstructural evolution of graphite. It also presents the foundations of a model for the graphite-molecular oxygen reaction system that is kinetically independent of graphite grade, and is capable of describing both the bulk and local oxidation rates under a wide range of conditions applicable to air-ingress.

J. Kane, Joshua↗

Molecular understanding of new-particle formation from α-pinene between -50 and +25 °C

Abstract. Highly oxygenated organic molecules (HOMs) contribute substantially to the formation and growth of atmospheric aerosol particles, which affect air quality, human health and Earth's climate. HOMs are formed by rapid, gas-phase autoxidation of volatile organic compounds (VOCs) such as α-pinene, the most abundant monoterpene in the atmosphere. Due to their abundance and low volatility, HOMs can play an important role in new-particle formation (NPF) and the early growth of atmospheric aerosols, even without any further assistance of other low-volatility compounds such as sulfuric acid. Both the autoxidation reaction forming HOMs and their NPF rates are expected to be strongly dependent on temperature. However, experimental data on both effects are limited. Dedicated experiments were performed at the CLOUD (Cosmics Leaving OUtdoor Droplets) chamber at CERN to address this question. In this study, we show that a decrease in temperature (from +25 to -50°C) results in a reduced HOM yield and reduced oxidation state of the products, whereas the NPF rates (J1.7 nm) increase substantially. Measurements with two different chemical ionization mass spectrometers (using nitrate and protonated water as reagent ion, respectively) provide the molecular composition of the gaseous oxidation products, and a two-dimensional volatility basis set (2D VBS) model provides their volatility distribution. The HOM yield decreases with temperature from 6.2 % at 25°C to 0.7 % at -50°C. However, there is a strong reduction of the saturation vapor pressure of each oxidation state as the temperature is reduced. Overall, the reduction in volatility with temperature leads to an increase in the nucleation rates by up to 3 orders of magnitude at -50°C compared with 25°C. In addition, the enhancement of the nucleation rates by ions decreases with decreasing temperature, since the neutral molecular clusters have increased stability against evaporation. The resulting data quantify how the interplay between the temperature-dependent oxidation pathways and the associated vapor pressures affect biogenic NPF at the molecular level. Our measurements, therefore, improve our understanding of pure biogenic NPF for a wide range of tropospheric temperatures and precursor concentrations.

54 ENVIRONMENTAL SCIENCES↗

BETO 2021 Peer Review - Inverse Bioproduct Design Through Machine Learning and Molecular Simulation

This work aims to identify performance advantaged bioproducts (PABPs) through property prediction, which will guide experimental synthesis. The impact of this work will be faster market adoption of bioproducts with greater performance relative to incumbent products. We have identified >106 bioproduct candidates, but only some will have superior performance to create a market pull. High-throughput property prediction, enabled by machine learning, and elucidation of structure-function relationships, enabled by molecular simulation, provide a hypothesis driven approach for down selection of candidate biomolecules to pursue experimentally. To enable machine learning and molecular simulation for bioproduct discovery, automated structure generation and embedding must capture relevant features for prediction, databases must cover domains applicable to biobased products, and best practices for simulation of polymer systems must be developed. To address these challenges, we have established bioproduct relevant datasets, developed high-throughput polymer structure generation, and built end-to-end neural networks that have predicted 8 properties for >1.4 x 106 biopolymers. A molecular simulation pipeline for building, running, and analyzing polymers and polymer additives is being used to predict performance and develop design principles of biobased products. In collaboration with the PABP synthesis project, these computational tools are guiding synthesis and informing design of PABPs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Methods of making porous molecular structures

Disclosed herein is a method of making a porous molecular structure from a solution comprising an insoluble metal containing material and a ligand-providing material. In some embodiments, the porous molecular structure can be a Metal-Organic Framework (MOF). Ionic metal salts are the most common type of metal source for MOF production, but dissolution of metal salts complicates solvent recycling and creates corrosion and oxidation issues through evolved nitrate and chloride anions. Elucidating information that leads toward more efficient production of these versatile nanomaterials, while extending the knowledge base of how MOFs form during reaction, is critical to advancing MOF materials into large-scale use. Disclosed herein are improved methods for controlled synthesis of porous molecular structures such as MOFs comprising a solution-based synthesis with insoluble metallic precursor.

Joshi, Jayraj Nayan↗

An idea to explore: How an interdisciplinary undergraduate course exploring a global health challenge in molecular detail enabled science communication and collaboration in diverse audiences

Abstract Communication and collaboration are key science competencies that support sharing of scientific knowledge with experts and non‐experts alike. On the one hand, they facilitate interdisciplinary conversations between students, educators, and researchers, while on the other they improve public awareness, enable informed choices, and impact policy decisions. Herein, we describe an interdisciplinary undergraduate course focused on using data from various bioinformatics data resources to explore the molecular underpinnings of diabetes mellitus (Types 1 and 2) and introducing students to science communication. Building on course materials and original student‐generated artifacts, a series of collaborative activities engaged students, educators, researchers, healthcare professionals and community members in exploring, learning about, and discussing the molecular bases of diabetes. These collaborations generated novel educational materials and approaches to learning and presenting complex ideas about major global health challenges in formats accessible to diverse audiences.

59 BASIC BIOLOGICAL SCIENCES↗

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗