Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational molecular design”

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 361 records · Page 20

Machine learning overcomes human bias in the discovery of self-assembling peptides

Peptide materials have a wide array of functions, from tissue engineering and surface coatings to catalysis and sensing. Tuning the sequence of amino acids that comprise the peptide modulates peptide functionality, but a small increase in sequence length leads to a dramatic increase in the number of peptide candidates. Traditionally, peptide design is guided by human expertise and intuition and typically yields fewer than ten peptides per study, but these approaches are not easily scalable and are susceptible to human bias. Here, in this work, we introduce a machine learning workflow—AI-expert—that combines Monte Carlo tree search and random forest with molecular dynamics simulations to develop a fully autonomous computational search engine to discover peptide sequences with high potential for self-assembly. We demonstrate the efficacy of the AI-expert to efficiently search large spaces of tripeptides and pentapeptides. The predictability of AI-expert performs on par or better than our human experts and suggests several non-intuitive sequences with high self-assembly propensity, outlining its potential to overcome human bias and accelerate peptide discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗

Navigating Transition-Metal Chemical Space: Artificial Intelligence for First-Principles Design

Conspectus The variability of chemical bonding in open-shell transition-metal complexes not only motivates their study as functional materials and catalysts but also challenges conventional computational modeling tools. Here, tailoring ligand chemistry can alter preferred spin or oxidation states as well as electronic structure properties and reactivity, creating vast regions of chemical space to explore when designing new materials atom by atom. Although first-principles density functional theory (DFT) remains the workhorse of computational chemistry in mechanism deduction and property prediction, it is of limited use here. DFT is both far too computationally costly for widespread exploration of transition-metal chemical space and also prone to inaccuracies that limit its predictive performance for localized d electrons in transition-metal complexes. These challenges starkly contrast with the well-trodden regions of small-organic-molecule chemical space, where the analytical forms of molecular mechanics force fields and semiempirical theories have for decades accelerated the discovery of new molecules, accurate DFT functional performance has been demonstrated, and gold-standard methods from correlated wavefunction theory can predict experimental results to chemical accuracy. The combined promise of transition-metal chemical space exploration and lack of established tools has mandated a distinct approach. In this Account, we outline the path we charted in exploration of transition-metal chemical space starting from the first machine learning (ML) models (i.e., artificial neural network and kernel ridge regression) and representations for the prediction of open-shell transition-metal complex properties. The distinct importance of the immediate coordination environment of the metal center as well as the lack of low-level methods to accurately predict structural properties in this coordination environment first motivated and then benefited from these ML models and representations. Once developed, the recipe for prediction of geometric, spin state, and redox potential properties was straightforwardly extended to a diverse range of other properties, including in catalysis, computational “feasibility”, and the gas separation properties of periodic metal–organic frameworks. Interpretation of selected features most important for model prediction revealed new ways to encapsulate design rules and confirmed that models were robustly mapping essential structure–property relationships. Encountering the special challenge of ensuring that good model performance could generalize to new discovery targets motivated investigation of how to best carry out model uncertainty quantification. Distance-based approaches, whether in model latent space or in carefully engineered feature space, provided intuitive measures of the domain of applicability. With all of these pieces together, ML can be harnessed as an engine to tackle the large-scale exploration of transition-metal chemical space needed to satisfy multiple objectives using efficient global optimization methods. In practical terms, bringing these artificial intelligence tools to bear on the problems of transition-metal chemical space exploration has resulted in ML-model assessments of large, multimillion compound spaces in minutes and validated new design leads in weeks instead of decades.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Smoothed particle hydrodynamics modeling and analysis of oxide reduction process for uranium oxides

A common kinetic feature for oxide reduction chemical/electrochemical processes is oxygen transport via a porous metallic layer, which has been considered as a rate-determining step for reducing uranium oxides to metallic uranium. Accounting this kinetic behavior must involve the resolution of the moving reactive interface between shrinking oxide and expanding metal phases. This study presents a numerical model using smoothed particles hydrodynamics (SPH) to effectively deal with the evolution of the shrinking core reaction interface and oxygen transport via mass transfer of lithium oxide (Li 2 O) species in multiple mass transfer domains. We successfully validated the proposed model against a theoretical derivation for the oxide reduction process handling a shrinking oxide core with molten salt and metal ash medium on a simple planar geometry. Armed with successful validation results, the model examined a realistic reactant geometry to extend the arguments beyond the one-dimensional analyses, allowing the proposed model to apply to general application scenarios with multi-dimensional geometry. Here, this study demonstrated that the proposed model could simulate and evaluate an arbitrary reaction basket design without iterative experimental trials, which is prohibitive for a scaled high-temperature molten salt study in an inert environment. The construct potentially provides not only deep insights on multiphysics behaviors governing the process dynamics but also a robust framework for evaluating and screening candidate basket designs in the most cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unification of nonequilibrium molecular dynamics and the mode-resolved phonon Boltzmann equation for thermal transport simulations

Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. Finally, the unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exploring novel HIV ‐1 reverse transcriptase inhibitors with drug‐resistant mutants: A double mutant surprise

Abstract HIV‐1 reverse transcriptase (RT) remains a key target for HIV drug development. As successful management of the disease requires lifelong treatment, the emergence of resistance mutations is inevitable, making development of new RT inhibitors, which remain effective against resistant variants crucial. To this end, previous computationally guided drug design efforts have resulted in catechol diether compounds, which inhibit wildtype RT with picomolar affinities and appear to be promising preclinical candidates. To confirm that these compounds remain potent against Y181C, a widespread mutation conferring resistance to first generation inhibitors, they were screened against the HIV‐1 N119 clinical isolate, reported as a Y181C single mutant. In comparison to a molecular clone with the same mutation, N119 appears less susceptible to inhibition by our preclinical candidate compounds. A more detailed sequencing effort determined that N119 was misidentified and carries V106A in combination with Y181C. While both indolizine and naphthalene substituted catechol diethers are potent against the classical Y181C single mutant, the addition of V106A confers more resistance against the indolizine derivatives than the naphthalene derivatives. Crystal structures presented in this study highlight key features of the naphthyl group, which allow these compounds to remain potent in the double mutant, including stronger interactions with F227 and less reliance on V106 for stabilization of the ethoxy‐uracil ring, which makes critical hydrogen bonds with other residues in the binding pocket.

Biochemistry & Molecular Biology↗

Desalination metamodels and a framework for cross-comparative performance simulations

There is an opportunity to save energy and reduce operational expenses when choosing a suitable desalination method aided by computational modeling. Existing models are not conducive to generalized comparisons between different desalination methods. Therefore, the work in this study developed metamodels for six desalination methods, grouped them into thermal and molecular transport families, and validated their predictive performance within 9% difference from published data. This validated framework allowed comparisons of desalination methods at their prescribed ranges of operational conditions that they were designed for. These conditions specify feed salinity ranges of 1.6 to 2.4 g/kg for Capacitive Deionization and Reverse Osmosis (RO), 2.8 to 4.2 g/kg for Electrodialysis, 28 to 42 g/kg for Thermovapor Compression and Humidification-Dehumidification, and 37 to 55 g/kg for Multi-Effect Distillation (MED). Despite different operational conditions, all models exhibit non-linear, positive correlation between energy consumption and system size in response to feed salinity and production rate. The framework is also employed in a cross-comparative analysis between MED and RO whose results suggest that energy intensity for MED is an order of magnitude greater than RO for the same operational conditions, but actual operational costs are comparable. Overall, the framework is ready for deployment in case studies of actual desalination plants.

42 ENGINEERING↗

Bridge-Mediated Metal-to-Metal Electron and Hole Transfer in a Supermolecular Dinuclear Complex: A Computational Study Using Quantum Electron–Nuclear Dynamics

Bimetallic electron donor–acceptor complexes can facilitate electron and energy transfer with excellent structural control through synthetic design. In this work, we investigate the photochemical dynamics in a Ru–Cu bimetallic complex after photoexcitation of the Ru-centered charge transfer state. The physical underpinnings of the metal-to-metal directional charge transfer process are unraveled via analyses of the quantum electronic dynamics and electron–nuclear trajectories. Here, the effects of molecular vibrations in the photoexcited state on the charge transfer processes are also analyzed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Bottom-Up Approach to Rational Design of Crystalline Materials: Investigation of Vibronic Coherences Underlying Exciton Dynamics in Semiconductors

In this project we uncovered structure-function relationships of donor-acceptor co-crystals used to develop next-generation optoelectronic devices. Unraveling the photodynamics of molecular crystalline materials poses many challenges for spectroscopy due to broad, overlapping features representing numerous underlying dynamical processes. This leads researchers to make many assumptions about the dynamics of a system in choosing an appropriate kinetic fitting model. Computationally, electronic structure methods are either prohibitively expensive or underdeveloped for computing the excited state structure of molecular materials, especially states that exhibit charge transfer. Researchers must therefore perform calculations of excited electronic states using truncated models of molecular materials. Here we present a joint experimental-theoretical approach to bridging the gap between the photodynamics of a molecular material and its constituent molecules. We focus our efforts on quantifying the timescales and mechanisms of photoexcitation in donor-acceptor co-crystals and donor-acceptor dimers where the lowest-lying excited state is characterized by charge transfer from the donor to the acceptor. We employ ultrafast UV pump, UV-Vis probe transient absorption spectroscopy to unravel the time-resolved spectroscopic signatures of the photodynamics in both the crystalline material and donor-acceptor dimers in solution. We perform electronic structure and excited state dynamics calculations of the dimers to inform kinetic fitting models and assign the spectral features. The photodynamics of the crystal vs. dimer systems have many similarities, enabling unprecedented insights into the formation and evolution of charge transfer excitons in the crystalline systems.

36 MATERIALS SCIENCE↗

Effect of Ambient Organic Acids on the Water Structure at ${\rm TiO}_{2}$ Interfaces

A molecular-level understanding of the effects of ambient organic compounds on the wettability of titanium dioxide ( ${\rm TiO}_{2}$ ) surfaces is relevant to many of its energy-related and environmental applications. Herein, we focus on two common atmospheric carboxylic acids, formic and acetic acid, and characterize their adsorption/ desorption at the aqueous interfaces of anatase and rutile ${\rm TiO}_{2}$ using molecular dynamics with an ab initio deep neural network potential. Our simulations show that these acids prefer to be localized in the interfacial water layers close to the ${\rm TiO}_{2}$ surface where they are stabilized by the interaction/exchange of their acid proton with a surface oxygen, rather than chemisorb at the surface Ti sites by displacing the adsorbed water. Notably, these acids make the surface of anatase hydrophobic, whereas the larger fraction of adsorbed water dissociation can offset their effect on rutile. Furthermore, these results provide a picture where carboxylic acids control the wettability of ${\rm TiO}_{2}$ largely through acid-base chemistry at the interface rather than chemisorption on the oxide surface, a finding that can help improve the design of self-cleaning surfaces and photocatalytic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Review of the Effects of Polymer Binder Properties on Microstructure and Irreversible Volume Growth of Plastic Bonded Explosives Formulations

The rational design of effective polymeric binders for the formulation of plastic-bonded explosives (PBX) is challenging due to their inherent compositional complexity. The composites comprise irregularly shaped energetic material (EM) powders coated with low weight fractions of polymer via non-equilibrium processes such as slurry coating. Defects can deleteriously affect PBX stability and performance: nano- to micrometer-scale voids can act as loci for hot spots, lowering deflagration and detonation temperatures in unpredictable ways. Furthermore, some nominally desirable polymer properties are at odds with each other: e. g. good flow characteristics are desirable for coating and adhesion, but mechanical stiffness is needed to prevent deformation and cracking of PBX under mechanical stress. Good binder adhesion is critical, but the best means to predict and measure adhesion in PBX is not obvious. Experimental methods of determining binder adhesion on model surfaces may not capture polymer structural configurations relevant to deposition during coating. Molecular dynamics-based computational models have predicted key observables in PBX formulation, suggesting that they may be powerful tools for binder selection. In this review, primarily recent (~2006 and later) literature on polymeric binders for insensitive HE (IHE) is surveyed. We focus on how binder properties influence observable PBX properties as resistance to irreversible volume growth and void formation in PBX formulations mainly (but not exclusively) featuring 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the EM (e. g. PBX 9502, LX-17 and others). Conclusions from these studies yield useful guidelines for choosing HE binder candidates, as well as for the general design of highly-filled polymer composite materials. Finally, studies describing challenges in PBX formulation with the newer and more energetically dense high explosive, LLM-105, will be discussed.

adhesion↗

A Machine Learning Model for Predicting Composition of Catalytic Coprocessing Products from Molecular Beam Mass Spectra

Demand for the development of an automated and integrated refining process for biofuels has increased in recent years due to the lack of generalized process inspection tools. In bio-oil upgrading processes, all process variables are maintained based on the offline specification of intermediates and products. A lack of real-time product specifications in batch-wise monitoring can cause process failure and wasted resources. Therefore, there is a need for a fast and accurate intermediates/product specification tool that can be used for real-time specification to reduce waste and mitigate the risk of process failure. Here, to address this gap, we developed a machine learning (ML) model for predicting speciated bio-oil composition, including paraffin, iso-paraffins, olefins, naphthene, and aromatics. The model is trained using the mass spectra from upgraded products collected in the vapor phase before condensation and predicts the composition of the condensed product. Training ML models using raw mass spectra is challenging due to numerous overlapped peaks originating from different parent compounds. With this in mind, we propose a protocol that (i) transforms raw mass spectra to chemistry-inspired predefined features and (ii) trains decision tree-based models using these features. Our results show that the random forest model was robust against overfitting and had the highest accuracy compared to other models. Moreover, a stochastic ablation method determined the eight most significant features while maximizing the accuracy. Our protocol facilitates real-time compositional analysis of upgraded bio-oils and thus real-time process monitoring. Additionally, this protocol enables the rational design of efficient catalysts and the determination of optimal process conditions.

09 BIOMASS FUELS↗

Hydrogen-Bonded Organic Frameworks: A Rising Class of Porous Molecular Materials

Hydrogen-bonded organic frameworks (HOFs) are a class of porous molecular materials that rely on the assembly of organic building blocks by means of hydrogen-bonding interactions to form two-dimensional (2D) and three-dimensional (3D) crystalline networks. The reversible nature of the hydrogen-bond formation endows HOFs with the attributes of solution processability and simple regeneration. High-quality single crystals of HOFs can be grown easily for unambiguous superstructure determination by single-crystal X-ray diffraction, which is crucial for the elucidation of superstructure–property relationships. During the past decade, considerable progress has been achieved in realizing stable HOFs with permanent porosities by focusing on the design of molecular building blocks in order to introduce rigidity, auxiliary [π···π] interactions, and interpenetration of their frameworks to sustain the extended networks. The applications of HOFs are far-reaching, spanning catalysis, energy, and biomedical products as well as the storage and separation of fine chemicals. In this paper, we, first of all, provide an overview of the chronological development of HOFs, starting from the seminal work by Marsh and Duchamp in 1969 on the crystal superstructure of the hydrogen-bonded networks of trimesic acid. We identify the development of novel hydrogen-bonding motifs such as diaminotriazine (DTA), the introduction of the concept of molecular tectonics, and the establishment of permanent porosity in HOFs as being some of the milestones, which incentivized the current burgeoning research endeavors on developing HOFs as multifunctional materials. This Account is focused primarily on surveying the strategies for constructing porous 3D HOFs based on organic building blocks with peripheral carboxyl groups. These strategies are presented in the following categories: (1) the polycatenation of 2D networks by trigonal building blocks to form global 3D frameworks, (2) the utilization of building blocks with 3D geometries—tetrahedral and trigonal prismatic—that are predisposed to form 3D networks, and (3) the docking by shape-fitting of geometrically labile building blocks. We emphasize how the molecular geometry of the building blocks plays an important role in modulating the superstructures of extended frameworks so as to address specific applications. Recognizing that the in silico design of HOFs is the ultimate goal of researchers in this field, we also discuss the recent advances in superstructure prediction that lead to the formation of porous supramolecular crystals and assess the complications in implementing computational methods for HOFs with complex superstructures. We hope this Account will inspire the development of new supramolecular designs and creative approaches to crystal engineering that aid and abet the assembly of multifunctional HOFs with customizable properties.

36 MATERIALS SCIENCE↗

Microstructural impacts on ionic conductivity of oxide solid electrolytes from a combined atomistic-mesoscale approach

Abstract Although multiple oxide-based solid electrolyte materials with intrinsically high ionic conductivities have emerged, practical processing and synthesis routes introduce grain boundaries and other interfaces that can perturb primary conduction channels. To directly probe these effects, we demonstrate an efficient and general mesoscopic computational method capable of predicting effective ionic conductivity through a complex polycrystalline oxide-based solid electrolyte microstructure without relying on simplified equivalent circuit description. We parameterize the framework for Li 7- x La 3 Zr 2 O 12 (LLZO) garnet solid electrolyte by combining synthetic microstructures from phase-field simulations with diffusivities from molecular dynamics simulations of ordered and disordered systems. Systematically designed simulations reveal an interdependence between atomistic and mesoscopic microstructural impacts on the effective ionic conductivity of polycrystalline LLZO, quantified by newly defined metrics that characterize the complex ionic transport mechanism. Our results provide fundamental understanding of the physical origins of the reported variability in ionic conductivities based on an extensive analysis of literature data, while simultaneously outlining practical design guidance for achieving desired ionic transport properties based on conditions for which sensitivity to microstructural features is highest. Additional implications of our results are discussed, including a possible connection between ion conduction behavior and dendrite formation.

25 ENERGY STORAGE↗

Potential for exciton condensation in a highly conductive amorphous polymer

An outstanding challenge in synthesis and theory is to develop molecular materials at ambient conditions that exhibit highly efficient energy transfer. Here, in this study, we demonstrate the potential of a recently synthesized, highly conductive amorphous material—a nickel tetrathiafulvalene-tetrathiolate (NiTTFtt) polymer—to become an exciton condensate—a Bose-Einstein condensate of particle-hole pairs, known as excitons, that supports dissipationless flow of excitation energy. While exciton condensates have recently been realized in ordered materials, we show by advanced electronic structure calculations that this highly correlated phenomenon can potentially be realized in molecularly tailored, amorphous materials. In contrast to the Bechgaard salts that support superconductivity at compressed geometries requiring high pressures, we show that the recently synthesized, amorphous NiTTFtt polymer exhibits the computational signature of exciton condensation at experimentally realizable geometries, occurring at ambient pressures. Results suggest that superfluidity in this system and related systems—including van der Waals structures, molecular metals with extended-TTF dithiolate ligands, and Bechgaard salts—may occur via a nontraditional excitonic mechanism tuneable according to system composition, geometry, size, and charge. This study prompts further experimental investigation of the rational design of molecularly scaled exciton condensates with potential applications to efficient transport in technologically relevant materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scalable photonic integrated circuits for high-fidelity light control

Advances in laser technology have driven discoveries in atomic, molecular, and optical (AMO) physics and emerging applications, from quantum computers with cold atoms or ions, to quantum networks with solid-state color centers. This progress is motivating the development of a new generation of optical control systems that can manipulate the light field with high fidelity at wavelengths relevant for AMO applications. These systems are characterized by criteria: (C1) operation at a design wavelength of choice in the visible (VIS) or near-infrared (IR) spectrum, (C2) a scalable platform that can support large channel counts, (C3) high-intensity modulation extinction and (C4) repeatability compatible with low gate errors, and (C5) fast switching times. Here, we provide a pathway to address these challenges by introducing an atom control architecture based on VIS-IR photonic integrated circuit (PIC) technology. Based on a complementary metal–oxide–semiconductor fabrication process, this atom-control PIC (APIC) technology can meet system requirements (C1)–(C5). As a proof of concept, we demonstrate a 16-channel silicon-nitride-based APIC with (5.8±0.4)ns response times and >30dB extinction ratio at a wavelength of 780 nm.

Menssen, Adrian J.↗

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗

teemi: An open-source literate programming approach for iterative design-build-test-learn cycles in bioengineering

Synthetic biology dictates the data-driven engineering of biocatalysis, cellular functions, and organism behavior. Integral to synthetic biology is the aspiration to efficiently find, access, interoperate, and reuse high-quality data on genotype-phenotype relationships of native and engineered biosystems under FAIR principles, and from this facilitate forward-engineering strategies. However, biology is complex at the regulatory level, and noisy at the operational level, thus necessitating systematic and diligent data handling at all levels of the design, build, and test phases in order to maximize learning in the iterative design-build-test-learn engineering cycle. To enable user-friendly simulation, organization, and guidance for the engineering of biosystems, we have developed an open-source python-based computer-aided design and analysis platform operating under a literate programming user-interface hosted on Github. The platform is called teemi and is fully compliant with FAIR principles. In this study we apply teemi for i) designing and simulating bioengineering, ii) integrating and analyzing multivariate datasets, and iii) machine-learning for predictive engineering of metabolic pathway designs for production of a key precursor to medicinal alkaloids in yeast. The teemi platform is publicly available at PyPi and GitHub.

59 BASIC BIOLOGICAL SCIENCES↗