Search NASA⌕ Search

SEARCH · Search NASA

Results for “script”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An Implementation of ICME in Materials Information Exchanging Interfaces

The Integrated Computational Materials Engineering (ICME) approach provides a new paradigm for improving the performance of existing materials or discovering and developing new materials. It focuses on developing effective connections between isolated engineering fields, to bring quantitative processing-structure-property relationships and abundant validated data that populate the knowledge base for accelerating the research of new materials while reducing the cost of development. The data exchanging interfaces play a key role in building such ICME connections among different materials models, simulation tools and individual organizations. With implementations of information exchanging interface between different materials database, and interface between database and applications, this paper concludes that in building effective ICME applications, 1) standards-compliant interface can improve the exchange efficiency; 2) popular web service enables automated online materials data transfer; 3) customized data interoperation scripts can provide flexibility and productivity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Volumetric analysis and mesh generation of real and artificial microstructural geometries

Producing a viable finite element mesh of realistic microstructural structural geometry is a critical step in analyzing the thermo-mechanical behavior of complex multi-material composites. Advancements in imaging technology such as micro computed tomography have allowed modelers to access high resolution mesoscale geometries for direct numerical simulation. However, converting from voxel based 3D images to usable finite element meshes has been challenging. A robust method including algorithms and software scripts for generating finite element meshes from 3D imaged microstructures is presented in this paper. It includes a routine for inserting cohesive elements around material interfaces to enable modeling of interface properties including delamination and damage. The algorithms and procedures presented in this method leverage currently available software packages for processing surface based geometry into volume based meshes. In addition to converting real geometry from physical imaging systems, algorithms for producing numerically generated and statistically equivalent microstructural geometry are also included. These artificial microstructures can be a valuable resource for modelers when physical specimens do not exist or are limited in quantity: Method establishes a workflow from voxel data to viable finite element mesh including interface information, and Includes a method for synthetic geometry generation based on metrics from real microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A practical guide to characterizing irradiated nuclear fuels using FIB tomography

We report focused ion beam (FIB) tomography with combined electron backscatter diffraction (EBSD) and energy dispersive x-ray spectroscopy (EDS) is a technique capable of statistically characterizing the microstructure and spatial compositional variation of nuclear fuel in three-dimensions (3D). The 3D visualization from FIB tomography provides a comprehensive picture of the interconnected microstructural and compositional features that can impact fuel performance. While these features are often characterized with surface examination, the complexity and relationship of fission products and grain boundary networks may not fully be captured by these 2D methods. This work presents a practical guide to FIB tomography that is tailored to nuclear fuel characterization. The steps used to collect and process the data are provided along with the scripts used to process the data. Additionally, suggestions for future characterization efforts utilizing this approach are given.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Microstructure and creep properties of cast near-eutectic Al–Ce–Ni alloys

This study investigates the as-cast and aged microstructures, thermal stability, ambient temperature strengthening, and creep resistance of three ternary Al–Ce–Ni alloys (wt%): near-eutectic Al–10Ce–5Ni (with both eutectic and hypoeutectic regions), hypoeutectic Al-7.5Ce-3.75Ni (with numerous primary Al dendrites), and hypereutectic Al-12.5Ce-6.25Ni (with coarse, blocky primary Al3Ni and Al11Ce3 precipitates and some primary Al dendrites). Depending on the alloy composition and local solidification conditions, the following eutectic morphologies are found: (i) coarse Al–Ce eutectic colonies where Al 11 Ce 3 is in the form of “Chinese script”, (ii) intermingled regions of binary Al–Ce and Al–Ni eutectic colonies, with finer Al 3 Ni and Al 11 Ce 3 fibers, (iii) large ternary eutectic colonies, where the binary Al 3 Ni and Al 11 Ce 3 phases are alternating or intertwining within the individual, fine fibers (diameters of ~60–170 nm, depending on solidification rates), and (iv) ternary eutectic zones (between primary Al dendrites), where fine Al 3 Ni and Al 11 Ce 3 build up a 3D-interconnected network. The high volume fraction of intermetallic phases and extremely fine eutectic spacing/fiber diameter both contribute to high ambient strengthening (higher as-cast microhardness than binary Al–Ce or Al–Ni), and also provide enhanced creep resistance at 300 and 350 °C. Additionally, the alloys are coarsening-resistant up to 425 °C for extended periods, with a gradual decrease in microhardness. The alloys aged at 400 °C to 1050 h show fiber fragmentation and coarsening of the resulting particles, with the faster-diffusing Ni driving more rapid coarsening of the Al 3 Ni particles which engulf finer, more stable Al 11 Ce 3 particles. Severe overaging (performed at 590 °C for 24 h) leads to Al 3 Ni and Al 11 Ce 3 spheroids which remain submicron-sized in eutectic colonies, but micron-sized at colony boundary and at Al dendrite-eutectic interface. Creep resistance at 300 °C of overaged Al–10Ce–5Ni remains substantial, consistent with load-transfer based composite strengthening being an important strengthening mechanism in these alloys, making them excellent candidates for replacement of heavier steel or titanium parts operating under stress up to 300 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Baghdad Atlas: A relational database of inelastic neutron-scattering (n,n ' γ) data

A relational database has been developed based on the original (n,n'γ) work carried out by A. M. Demidov et al., at the Nuclear Research Institute in Baghdad, Iraq (Demidov et al., 1978) for 105 independent measurements comprising 76 elemental samples of natural composition and 29 isotopically-enriched samples. The information from this Atlas includes: γ-ray energies and relative intensities; nuclide and level data corresponding to the residual nucleus and meta data associated with the target sample that allows for the extraction of the flux-weighted (n,n'γ) cross sections for a given transition relative to a defined value. The optimized angular-distribution-corrected fast-neutron flux-weighted partial γ-ray cross section for the production of the 846.8-keV 21+→0gs+γ-ray transition in 56Fe, determined to be $\langle$σγ$\rangle$=143(29) mb, is used for this purpose. However, different values for the adopted cross section can be readily implemented to accommodate user preference based on revised determinations of this quantity. The Atlas (n,n'γ) data has been compiled into a series of CSV-style ASCII data sets and a suite of Python scripts have been developed to build and install the database locally. The database can then be accessed directly through the SQLite engine, or using alternative methods such as the Jupyter Notebook Python-browser interface. Several examples exploiting different interaction methodologies are distributed with the complete software package.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Test and extraction methods for the QC parameters of silicon strip sensors for ATLAS upgrade tracker

The Quality Control (QC) of pre-production strip sensors for the Inner Tracker (ITk) of the ATLAS Inner Detector upgrade has finished, and the collaboration has embarked on the QC test programme for production sensors. This programme will last more than 3 years and comprises the evaluation of approximately 22000 sensors. 8 Types of sensors, 2 barrel and 6 endcap, will be measured at many different collaborating institutes. The sustained throughput requirement of the combined QC processes is around 500 sensors per month in total. Measurement protocols have been established and acceptance criteria have been defined in accordance with the terms agreed with the supplier. For effective monitoring of test results, common data file formats have been agreed upon across the collaboration. To enable evaluation of test results produced by many different test setups at the various collaboration institutes, common algorithms have been developed to collate, evaluate, plot and upload measurement data. This allows for objective application of pass/fail criteria and compilation of corresponding yield data. These scripts have been used to process the data of more than 3000 sensors so far, and have been instrumental for identification of faulty sensors and monitoring of QC testing progress.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

A methodology for customizing implantation profiles of light ions using a single thin foil energy degrader

We report a method was developed to quantify the spatial distribution and implantation depth of energy-degraded light ions with a thin foil rotating energy degrader for use during multiple ion beam irradiation. The methodology covers three physical phenomena: ions passing through the thin foil, ions travelling through the vacuum to the target, and ion implantation into the target, and accounts for the distribution of ions both in depth and in plane. The processes of energy straggling and scattering were calculated using SRIM. The effects of raster-scanning, and the geometry of the system were implemented in scripts handling the SRIM output files. Elastic backscattering (EBS) using 2.38 MeV H+ protons was used to measure the helium depth profiles after implantation with and without thin foil energy degradation. Defect analysis with transmission electron microscopy confirmed the implantation profiles measured with EBS and calculated with SRIM.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

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

FEMA Asteroid Impact Tabletop Exercise Simulations

We describe the computational simulations and damage assessments that we provided in support of a tabletop exercise (TTX) at the request of NASA's Near-Earth Objects Program Office. The overall purpose of the exercise was to assess leadership reactions, information requirements, and emergency management responses to a hypothetical asteroid impact with Earth. The scripted exercise consisted of discovery, tracking, and characterization of a hypothetical asteroid; inclusive of mission planning, mitigation, response, impact to population, infrastructure and GDP, and explicit quantification of uncertainty. Participants at the meeting included representatives of NASA, Department of Defense, Department of State, Department of Homeland Security/Federal Emergency Management Agency (FEMA), and the White House. The exercise took place at FEMA headquarters. Sandia's role was to assist the Jet Propulsion Laboratory (JPL) in developing the impact scenario, to predict the physical effects of the impact, and to forecast the infrastructure and economic losses. We ran simulations using Sandia's CTH hydrocode to estimate physical effects on the ground, and to produce contour maps indicating damage assessments that could be used as input for the infrastructure and economic models. We used the FASTMap tool to provide estimates of infrastructure damage over the affected area, and the REAcct tool to estimate the potential economic severity expressed as changes to GDP (by nation, region, or sector) due to damage and short-term business interruptions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

A new measurement of amine steric hindrance – N exposure

Here, a new measurement of steric hindrance effect of amines, N exposure, is proposed. N exposure outperforms traditional parameters as it is independent of group-specific parame-ters and accounts for the spatial arrangements of pendant groups. Such a measurement can facilitate the screening of sterically hindered amines for carbon capture with improved CO 2 loading capacities. Along with this study, a script is provided to quantify the N exposure for other amine molecules of interest.

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

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A FAIR-compliant parts catalogue for genome engineering and expression control in Saccharomyces cerevisiae

The synthetic biology toolkit for baker's yeast, Saccharomyces cerevisiae, includes extensive genome engineering toolkits and parts repositories. However, with the increasing complexity of engineering tasks and versatile applications of this model eukaryote, there is a continued interest to expand and diversify the rational engineering capabilities in this chassis by FAIR (findable, accessible, interoperable, and reproducible) compliance. In this study, we designed and characterised 41 synthetic guide RNA sequences to expand the CRISPR-based genome engineering capabilities for easy and efficient replacement of genomically encoded elements. Moreover, we characterize in high temporal resolution 20 native promoters and 18 terminators using fluorescein and LUDOX CL-X as references for GFP expression and OD600 measurements, respectively. Additionally, all data and reported analysis is provided in a publicly accessible jupyter notebook providing a tool for researchers with low-coding skills to further explore the generated data as well as a template for researchers to write their own scripts. We expect the data, parts, and databases associated with this study to support a FAIR-compliant resource for further advancing the engineering of yeasts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

UEMtomaton: A Source-Available Platform to Aid in Start-up of Ultrafast Electron Microscopy Labs

The steady rise in the number of ultrafast electron microscopy (UEM) labs, in addition to the opacity and lack of detailed descriptions of current approaches that would enable point-by-point construction, has created an opportunity for sharing common methods and instrumentation for (for example) automating data acquisition to assist in efficient lab start-up and to learn about common and robust protocols. In the spirit of open sharing of methods, we provide here a description of an entry-level method and user interface (UI) for automating UEM experiments, and we provide access to the source code and scripts (source-available) for ease of implementation or as a starting reference point for those entering or seeking to enter the field (https://github.com/CEMSFlannigan/UEMtomaton/releases/tag/v1.0). Core instrumentation and physical connections in the UEM lab at Minnesota are described. Interface communication schemes consisting of duo server-client pairs between critical components – the optical delay stage and the UEM digital camera – are presented, with emphasis placed on describing the logic and communications sequence designed to conduct automated series acquisitions. An application designed and programmed with C++/CLI as Windows Forms in Microsoft Visual Studio – dubbed UEMtomaton – is also presented. Key to the UI layout is centralization of the automation tasks and establishment of communication within the software rather than by interfacing with each individual workstation. It is our hope that this note provides useful insight for current and future UEM researchers, particularly with respect to generalizability and portability of the approach to emerging labs. Here, we note that while this basic, entry-level approach is certainly not the most sophisticated or comprehensive of those currently in use, we feel there is nevertheless value in clearly communicating a proven straightforward method to hopefully lower the barrier to entry into the field.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Systems biology approach to study the role of miRNA in promoter targeting during megakaryopoiesis

Highlights: • miRNAs could regulate Cyclins and CDKs by targeting their promoters. • Perl programming was used to identify putative targets of miRNA in promoters. • miR-1273g-3p and miR-619-5p showed significance with longest seed sequence matches. • Through CDK10, CDK11 and Cyclin, these two miRNAs could regulate cell cycle during megakaryopoiesis. The distinct process of megakaryopoiesis requires occurrence of endomitosis for polyploidization of the megakaryocytes. Although, Cyclins, CDKs and have been described to regulate endomitosis, the exact mechanism still remains an enigma. miRNA which were otherwise known as post transcriptional gene silencers are now emerging with various non-canonical functions including gene regulation at pre-transcriptional level by miRNA binding at promoter region. Out of the many processes they regulate, miRNA have been manifested to play a role in megakaryocyte differentiation. In this study an attempt has been made to identify miRNA that could regulate cell cycle genes (Cyclins and CDKs) by targeting their promoters, during megakaryopoiesis. A new computational algorithm was implemented using Perl programming to identify putative targets of miRNA in CDK and Cyclin promoters. Perl script was also used to check nuclear localizing miRNA based on the presence of a consensus sequence. Real-time PCR was performed to analyze the expression of miRNA and their predicted targets in Dami vs. PMA treated Dami cells. Putative targets of miRNAs with longest, high complementarity matches in CDK/Cyclin promoters were obtained. We identified two significant miRNA, miR-1273g-3p and miR-619-5p with longest seed sequence matches. We further identified three main targets (CDK10, CDK11, Cyclin F) through which these two miRNA could regulate cell cycle during megakaryopoiesis. Our results reinforce the role of promoting targeting miRNA in regulation of cell cycle through certain CDK/Cyclins to support the process of endomitosis during megakaryopoiesis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Comparing Sensor Fusion and Multimodal Chemometric Models for Monitoring U(VI) in Complex Environments Representative of Irradiated Nuclear Fuel

Optical sensors and chemometric models were leveraged for the quantification of uranium(VI) (0–100 μg mL –1 ), europium (0–150 μg mL –1 ), samarium (0–250 μg mL –1 ), praseodymium (0–350 μg mL –1 ), neodymium (0–1000 μg mL –1 ), and HNO 3 (2–4 M) with varying corrosion product (iron, nickel, and chromium) levels using laser fluorescence, Raman scattering, and ultraviolet–visible–near-infrared absorption spectra. In this paper, an efficient approach to developing and evaluating tens of thousands of partial least-squares regression (PLSR) models, built from fused optical spectra or multimodal acquisitions, is discussed. Each PLSR model was optimized with unique preprocessing combinations, and features were selected using genetic algorithm filters. The 7-factor D-optimal design training set contained just 55 samples to minimize the number of samples. The performance of PLSR models was evaluated by using an automated latent variable selection script. PLS1 regression models tailored to each species outperformed a global PLS2 model. PLS1 models built using fused spectra data and a multimodal (i.e., analyzed separately) approach yielded similar information, resulting in percent root-mean-square error of prediction values of 0.9–5.7% for the seven factors. Further, the optical techniques and data processing strategies established in this study allow for the direct analysis of numerous species without measuring luminescence lifetimes or relying on a standard addition approach, making it optimal for near-real-time, in situ measurements. Nuclear reactor modeling helped bound training set conditions and identified elemental ratios of lanthanide fission products to characterize the burnup of irradiated nuclear fuel. Leveraging fluorescence, spectrophotometry, experimental design, and chemometrics can enable the remote quantification and characterization of complex systems with numerous species, monitor system performance, help identify the source of materials, and enable rapid high-throughput experiments in a variety of industrial processes and fundamental studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Practical Guide to Chemometric Analysis of Optical Spectroscopic Data

The methodology and mathematical treatment of several classic multivariate methods for the analysis of spectroscopic data is demonstrated in a straightforward way that can be used as a basis for teaching an undergraduate introductory course on chemometric analysis. The multivariate techniques of classical least squares (CLS), principal component regression (PCR), and partial least squares (PLS), as well as the univariate Beer’s law method have been described and compared, building students’ understanding by starting with the univariate method and progressing step by step into the multivariate methods. Equations for the production of regression vectors from training set spectral data is described and their use demonstrated for the prediction of constituent concentrations on a separate validation set of spectra. Extreme care is taken to ensure consistency in variable formatting of data matrices. This provides a key foundation to understanding how spectral data are manipulated using these different mathematical approaches for building quantitative regression models. Each method is applied to a real-world data set, and the results are discussed to show students the types of information that can be gleaned from each method. A training set comprised of 20 infrared absorbance spectra containing 3 constituents (benzene, polystyrene, and gasoline) of known composition are used to demonstrate the matrix operations for each regression method. A separate set of 12 real-world napalm samples (containing benzene, polystyrene and gasoline) are used as a validation set to demonstrate the ability to utilize the regression models on an unknown dataset. A toolbox (PNNL Chemometric Toolbox) written in MATLAB language is supplied in the Supplemental Information file and can be used as a companion for understanding the development and deployment of the chemometric algorithms described in this paper. The datasets of the infrared spectra are also supplied, allowing users to build and inspect the chemometric models on their own. Finally, the Toolbox includes scripts to assist users in loading their own datasets into MATLAB and performing CLS, PCR, and PLS on their data.

Upper-Division Undergraduate, Analytical Chemistry↗
DOE OSTI2025

Catalight─An Open-Source Automated Photocatalytic Reactor Package Illustrated through Plasmonic Acetylene Hydrogenation

An open-source and modular Python package, Catalight, is developed and demonstrated to automate (photo)catalysis measurements. (Photo)catalysis experiments require studying several parameters to evaluate performance, including the temperature, gas flow rate and composition, illumination power, and spectral profile. Catalight orchestrates measurements over this complicated parameter space and systematically stores, analyzes, and visualizes the results. To showcase the capabilities of Catalight, we perform an automated apparent activation barrier measurement of acetylene hydrogenation over a plasmonic AuPd catalyst on an Al 2 O 3 support, simultaneously varying laser power, wavelength, and temperature in a multiday experiment controlled by a simple Python script. Our chemical results unexpectedly show an increased activation barrier upon light excitation, contrary to previous findings for other plasmonic reactions and catalysts. We show that the reaction rate order with respect to both acetylene and hydrogen remains unchanged upon illumination, suggesting that molecular surface coverage is not changed by light. By analyzing the inhomogeneity of the laser-induced heating, we attribute these results to a partial photothermal effect combined with a photochemical/hot electron-driven mechanism. In conclusion, our findings highlight the capabilities of a new experiment automation tool; explore the photocatalytic mechanism for an industrially relevant reaction; and identify systematic sources of error in canonical photocatalysis experimental procedures.

Catalysts↗
DOE OSTI2018

Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis

We report that replica-exchange molecular dynamics (REMD) has been proven to efficiently improve the convergence of free energy perturbation (FEP) calculations involving considerable reorganization of their surrounding. We previously introduced the FEP/(λ,H)-REMD algorithm for ligand binding, in which replicas along the alchemical thermodynamic coupling axis λ were expanded as a series of Hamiltonian boosted replicas along a second axis to form a two-dimensional (2D) replica-exchange exchange map [Jiang, W.; Roux, B., J. Chem. Theory Comput. 2010, 6 (9), 2559-2565]. Aiming to achieve a similar performance at a lower computational cost, we propose here a modified version of this algorithm in which only the end-states along the alchemical axis are augmented by boosted replicas. The reduced FEP/(λ,H)-REMD method with one-dimensional (1D) unbiased alchemical thermodynamic coupling axis λ is implemented on the basis of generic multiple copy algorithm (MCA) module of the biomolecular simulation program NAMD. The flexible MCA framework of NAMD enables a user to design customized replica-exchange patterns through Tcl scripting in the context of a highly parallelized simulation program without touching the source code. Two Hamiltonian tempering boosting scheme were examined with the new algorithm: a first one based on potential energy rescaling of a pre-identified “solute”, and a second one via the introduction of flattening torsional free energy barriers. As two illustrative examples with reliable experiment data, the absolute binding free energies of pxylene and n-butylbenzene to the nonpolar cavity of the L99A mutant of T4 lysozyme were calculated. Lastly, the tests demonstrate that the new protocol efficiently enhances the sampling of torsional motions for backbone and side chains around the binding pocket and accelerates the convergence of the free energy computations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗