Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Neutron diffraction discriminates between models for the nanoarchitecture of graphene sheets in glassy carbon

Glassy carbon is a chemically inert, biocompatible, disordered material with the graphene sheet as its basic building block. Structural characterisation techniques have so far been unable to provide definitive distinction between the many proposed models for its structure. Computer based simulation methods have made a step forward by predicting structure from interaction potentials, but the results are sensitive to the choice of the potential. Here we use the white neutron beam of the Spallation Neutron Source at Oak Ridge National Laboratory coupled with the SNAP time-of-flight diffractometer to sample the reciprocal space of glassy carbon and calculate accurate radial distribution functions. From the radial distribution information, we determine graphene sheet dimensions, registration, curvature, and assess defect content to rank the proposed structures in terms of their agreement with experimental data. In this work, we find that the recent models derived from empirical potentials give the best agreement.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗
DOE OSTI2019

Generation and characterization of carbon fiber microstructure in atomistic simulations

A clear understanding of carbon fiber (CF) microstructure is necessary for the development of high strength CFs. Here, we present an atomistic approach for generating and characterizing realistic mi- crostructures of CFs. Large-scale reactive molecular dynamics simulations are used to generate a set of distinct CF microstructures. Comprehensive characterization of the simulated microstructures is enabled by the development of a suite of computational structural analysis tools capable of evaluation of hy- bridization states of carbon atoms, populations and orientations of individual carbon rings, degree of graphitization, and pore size distribution. The calculation of X-ray diffraction profiles provides a direct link between the structural features of simulated samples and experimental data available for CFs. The CF generation algorithm is shown to produce microstructures with experimental densities and with structural characteristics matching those of PAN-based CFs. The key structural features affecting the properties of CFs, such as the relative fractions of graphitic, turbostratic and amorphous micro- constituents, degree of alignment, pore size distributions, and chemical cross-linking can be effectively controlled in simulations, thus enabling efficient exploration of structureeproperties relationships in CFs. The capabilities of the developed approach are illustrated by performing computational analysis of the mechanical deformation and fracture of CFs under axial tensile loading.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural basis of transcription: RNA polymerase II substrate binding and metal coordination using a free-electron laser

Catalysis and translocation of multisubunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near-atomic resolution and precise arrangement of key active site components have been elusive. Here, we present the free-electron laser (FEL) structures of a matched ATP-bound Pol II and the hyperactive Rpb1 T834P bridge helix (BH) mutant at the highest resolution to date. The radiation-damage-free FEL structures reveal the full active site interaction network, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg 2+ , and, more importantly, a putative third (site C) Mg 2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structures indicate that the third Mg 2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/BH interactions induce conformational changes that could allow translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the T834P mutant reveal rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Science & Technology - Other Topics↗
DOE OSTI2022

Role of Molecular Simulations in the Design of Metal–Organic Frameworks for Gas-Phase Thermocatalysis: A Perspective

Here, metal organic frameworks (MOFs) are highly tunable porous crystalline solids with spatially and electronically isolated catalytically active sites that have been demonstrated for a variety of thermo-, redox-, and photocatalytic reactions. Their tunable natures and relatively well-defined active sites make them advantageous for catalyst design, and molecular simulations have proven highly valuable in this endeavor. However, complexities in the MOF structure require advanced simulation strategies that accurately capture quantum chemistry at the strongly correlated transition metal cation active sites, compositional and structural changes caused by finite temperature and pressure reaction conditions, and transport effects in variously sized pore environments. Luckily, several groups have started using such simulation strategies, paving the way for rich opportunities to design MOF catalysts. Herein, we highlight such examples and provide a perspective on the needs for molecular simulations in the design of MOF catalysts moving forward.

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

Quantifying Membrane Structure and Dynamics during Bioproduct Production in Zymomonas mobilis by Molecular Simulation

The conversion of lignocellulosic biomass into biofuels and bioproducts by microbial biorefineries is central to a sustainable chemical industry. Zymomonas mobilis is one such biorefinery chassis and is resistant to ethanol stress, leading to its use in biomass conversion to biofuels and bioproducts. However, Z. mobilis growth is often inhibited by organic acids, aldehydes, alcohols, ketones, and amides found in biomass hydrolysate. The resulting slow growth inhibits production and as a result drives up the price for the resulting products. One hypothesis is that these molecules interact with or disrupt the bacterial membrane, triggering stress responses and hindering growth. To test this hypothesis at the molecular level, we employ all-atom molecular dynamics (MD) simulations to investigate lignocellulose-derived small molecules and their impact on a biologically relevant Z. mobilis membrane model. Simulations were conducted across a range of inhibitor concentrations from 0 to 2.5 mol %, analyzing key membrane properties such as area per lipid (APL), membrane thickness, lipid-order parameter (−S CH ), lateral diffusion coefficient (D xy ), and permeability coefficient (Pm). From simulation, we observed altered membrane structure and dynamics at these modest small molecule concentrations commonly found in hydrolysates. Generally, the membranes become thinner, with a higher area per lipid and lower-order parameter as the small molecule concentration increases. These trends are stronger for more hydrophobic molecules with greater hydrophobic bulk, as isobutanol, propanol, and propanoic acid showed greater membrane perturbations as the concentration increased compared to other small molecules. Tracking small molecule distributions directly in our equilibrium simulations allows us to determine concentration-dependent free energy profiles for these molecules. While the trends are noisy, generally the barriers to crossing the membrane decrease as the concentration increases, indicating that the membranes become leakier as small molecule concentrations rise. Comparing between native Z. mobilis membranes with hopanoids and membranes sharing the same phospholipid composition but without hopanoids, hopanoids stabilize and order the membrane for smaller molecules to maintain membrane structure but appear insufficient for larger hydrophobic molecules like isobutanol. These findings provide a mechanistic understanding of how small molecules found in biomass degradation streams interact with the Z. mobilis membrane, offering valuable insights for future strain engineering efforts to optimize biofuel and bioproduct synthesis from biomass feedstocks by highlighting limits to small molecule tolerance. This knowledge can guide the modification of membrane composition to develop more robust microbes, thereby improving microbial survival and yields in industrial contexts.

Singh, Nitin Kumar [Michigan State Univ., East Lan↗
DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Real-time dynamics and structures of supported subnanometer catalysts via multiscale simulations

Understanding the performance of subnanometer catalysts and how catalyst treatment and exposure to spectroscopic probe molecules change the structure requires accurate structure determination under working conditions. Experiments lack simultaneous temporal and spatial resolution and could alter the structure, and similar challenges hinder first-principles calculations from answering these questions. Here, we introduce a multiscale modeling framework to follow the evolution of subnanometer clusters at experimentally relevant time scales. We demonstrate its feasibility on Pd adsorbed on CeO 2 (111) at various catalyst loadings, temperatures, and exposures to CO. We show that sintering occurs in seconds even at room temperature and is mainly driven by free energy reduction. It leads to a kinetically (far from equilibrium) frozen ensemble of quasi-two-dimensional structures that CO chemisorption and infrared experiments probe. CO adsorption makes structures flatter and smaller. High temperatures drive very rapid sintering toward larger, stable/metastable equilibrium structures, where CO induces secondary structure changes only.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Future Electron-Ion Collider

Generalized Parton Distributions include rich information and became a powerful tool for studying hadron structure. Deeply Virtual Compton Scattering is the golden channel to access GPDs. The development of high luminosity and high-acceptance detectors (Electron-Ion Collider) allows physicists to overcome the difficulty of DVCS measurements. The outstanding performance of the accelerator performs the eA collisions at the center of mass energy from 20 to 140 GeV, while the luminosity at the scale of 1034 cm?2s?1. The new type of high-density crystal, PWO-II, produced by CRYTUR in 2×2×20 cm3 will be utilized in the calorimeter. The transparency of PWO-II crystals is > 70% at 620 nm, > 60% at 420 nm, and >35% at 360 nm. The crystals also show good radiation hardness under 30 Gy radiation exposure. Ten PWO-II crystals have a uniform light yield at 30 p.e./MeV, providing sufficient light yield to reduce the fluctuation of the energy measurements. Near 3000 crystals were constructed inside the 12-sided polygon supporting structure in the simulation. The island clustering algorithm was used for reconstructed the energy. The energy resolution study by the particle gun shows the stochastic term and constant term are 1.8% and 1.2%, respectively. The spatial resolution of NEEMC varies from 5% to 15% of the crystal?s width depending on the particle?s incident angle. The pion rejection of NEEMC can reach 103 with electron efficiency > 85% when the particle?s energy is larger than 1GeV. The Pi0- identified efficiency study can be interpreted as finding the local maxima in the single cluster caused by two high energy close photons. The study results show that efficiency is nearly 100% for pi0 energy is smaller than 10 GeV, and efficiency drops to 30% with 20 GeV pi0. The new type of 3x3 pixelated AC-LGAD were wire bonded to ALTIROC for performance test. The cross-talk between the adjacent channels is about 20% for VPA and 10% for TZ. A series of the TDC characteristic measurements show that the jitter for both preamplifiers is about 20 ps, and the time-walk effect is mild for inject charge > 12 pF. Furthermore, the beta source radiation results quantify the sharing scale of the 3x3 pixels AC-LGAD (? 20%). The electronics simulation study results show no significant changes in the spatial resolution, whether ADC resolution is 8, 10, or 12 bits. The 8-bit ADC is decided to use in the EICROC as it has a smaller size and power consumption than the 10-bit and 12-bit ADC. The ECCE is one of the full detector proposals for EIC. The new type of DVCS generator, called TOPEG, is used to generate the high acceptance beam configurations of 18×110 GeV2 electron and 4He. The acceptance study shows a -1.8 < ? < -1.4 gap between the BEMC and EEMC. The 10?x,y geometry cut is applied on the Roman Pots, so the acceptance of Roman Pots quickly drops to 0% when the polar angle of the recoiled 4He < 2 mrad. This primary ECCE simulation study suggests extending the acceptance of BEMC longitudinally as the structure limits the radial size of EEMC. The acceptance of the Roman Pots is still challenging. In 2023, The overall design of ePIC was finalized by merging two full detector proposals (ECCE and ATHENA) after a series of intensive simulation studies.

Wang, Pu-Kai↗
DOE OSTI2025

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Convection–Kelvin Wave Coupling in a Global Convection-Permitting Model

A convectively coupled equatorial Kelvin wave (CCKW) was observed over the equatorial Indian Ocean in early November 2011 during the DYNAMO field campaign. This study examines the structure of the CCKW event using two simulations made using the MPAS model: one with 3-km grid spacing without convective parameterization and another with a 15-km grid and parameterized convection. Both simulations qualitatively capture the observed structure of the CCKW, including its vertical tilt and progression of cloud/precipitation structures. The two simulations, however, differ substantially in the amplitude of the CCKW-associated precipitation. While the 3-km run realistically captures the observed modulation of precipitation by the CCKW, the 15-km simulation severely underestimates its magnitude. To understand the difference between the two MPAS simulations regarding wave–convection coupling within the CCKW, the relationship of precipitation with convective inhibition, saturation fraction, and surface turbulent fluxes is investigated. Results show that the 15-km simulation underestimates the magnitude of the CCKW precipitation peak in association with its unrealistically linear relationship between moisture and precipitation. Precipitation, both in observations and the 3-km run, is predominantly controlled by saturation fraction and this relationship is exponential. In contrast, the parameterized convection in the 15-km run is overly sensitive to convective inhibition and not sensitive enough to environmental moisture. The implications of these results on CCKW theories are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Structure–property relations of sodium iron phosphate nuclear waste glasses: Effects of iron redox ratio and glass composition

Iron phosphate glasses, known for their exceptional chemical durability and potential applicability in nuclear waste management, have gained significant attention over the years. The structures of these glasses are complicated by the coexistence of Fe 3+ and Fe 2+ , which plays a crucial role in determining their structures and properties. Here, this work uses molecular dynamics simulations to study the structural changes in Na 2 O–Fe 2 O 3 –P 2 O 5 glasses with varying glass composition and Fe 2+ /Fe 3+ redox ratio. It was found that the redox ratio and modifier contents significantly affected the short-range and medium-range orders in the glasses. Significant changes in the local environments around P 5+ and Fe 3+ were observed, as reflected by the bond distances and coordination numbers. Na + cations are found to preferentially associate with Fe 3+ (rather than Fe 2+ ), whereas Fe 2+ has stronger association with P 5+ than Na + , confirming the structural role of Fe 2+ as a glass modifier. The disruptions in P–O–P linkages upon increasing FeO suggest that FeO causes glass depolymerization. These glasses achieved higher connectivity with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, conerting phosphorous Q 2 to Q 3 units and iron Q 5 units to Q 4 units. The decrease of nonbridging oxygen fractions with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, through creating P–O–Fe linkages, is the main reason of enhanced network connectivity. Quantitative structure–property relationship analyses with different structural descriptors were used to correlate with measured properties. The analyses provided valuable insights into structure–property relationships, emphasizing the importance of choosing relevant energy parameters and defining glass network connectivity, particularly in F net descriptors. It was found the Fe–O–P linkage density exhibits strong correlations to measured dissolution rates, supporting the importance of these linkages in improving the chemical durability in iron phosphate glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

On the Structural Transformation of Ni/BaH 2 During a N 2 -H 2 Chemical Looping Process for Ammonia Synthesis: A Joint In Situ Inelastic Neutron Scattering and First-Principles Simulation Study

The demand for decarbonizing the ammonia industry by using renewable energy has invoked increasing research interests into catalyst development for effective N 2 reduction under mild conditions. Hydride-based materials are among some of the emerging catalysts for ammonia synthesis at ambient pressure and low temperatures (< 673 K). A recent chemical looping process based on Ni/BaH 2 showed the most promise as it can realize ammonia production at a temperature as low as 373 K and under ambient pressure. However, the chemical transformation of the hydride catalyst at the molecular level remains unclear in this process. Here, we report detailed in situ neutron spectroscopy and diffraction investigations along with first-principles simulations on the structural transformation of Ni/BaH 2 during the nitridation and hydrogenation steps in the chemical looping process for ammonia synthesis. It was shown that a ball-milling process of the starting Ni/BaH 2 could significantly decrease the size of BaH 2 and increase the density of defects, thus potentially enhancing the reactivity of the hydride. The evolution from BaH 2 to barium imide (BaNH) was evidenced in the inelastic neutron scattering (INS) and neutron diffraction results during the N 2 reaction step. During the hydrogenation study, in addition to the recovery of BaH 2 , a possible intermediate species, N-deficient barium imide, was also detected. In comparing the N 2 and H 2 reaction steps, the neutron results indicate that the hydrogenation step appears more difficult than the nitridation step, confirming the facile N 2 fixation property of Ni/BaH 2 catalyst in ammonia synthesis.

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

Insights into microbial community structure and function from a shallow, simulated CO 2 -leakage aquifer demonstrate microbial selection and adaptation

Geological carbon storage is likely to be a part of a comprehensive strategy to minimize the atmospheric release of carbon dioxide (CO 2 ), raising concerns that injected CO 2 will leak into overlying freshwater aquifers. CO 2(aq) leakage may impact the dominant microbial community responsible for important ecosystem functions such as nutrient cycling, metal cycling and carbon conversion. Here, we examined the impact of an experimental in situ CO 2 leakage on a freshwater aquifer microbial community. High-throughput 16S rRNA gene sequencing demonstrated lower microbial diversity in freshwater wells with CO 2 concentrations above 1.15 g/L. Metagenomic sequencing and population genome binning were used to evaluate the metabolic potential of microbial populations across four CO 2 exposed samples and one control sample. Population genome binning resulted in the recovery and annotation of three metagenome assembled genomes (MAGs). Two of the MAGs, most closely related to Curvibacter and Sulfuricurvum, had the functional capacity for CO 2 utilization via carbon fixation coupled to sulfur and iron oxidation. The third draft genome was an Archaea, most closely related to Methanoregula, characterized by the metabolic potential for methanogenesis. Together, these findings show that CO 2 leakage in a freshwater aquifer poses a strong selection, driving both microbial community structure and metabolic function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Azimuthally averaged wind and thermodynamic structures of tropical cyclones in global climate models and their sensitivity to horizontal resolution

Characteristics of tropical cyclones (TCs) in global climate models (GCMs) are known to be influenced by details of the model configurations, including horizontal resolution and parameterization schemes. Understanding model-to-model differences in TC characteristics is a prerequisite for reducing uncertainty in future TC activity projections by GCMs. This study performs a process-level examination of TC structures in eight GCM simulations that span a range of horizontal resolutions from 1° to 0.25°. A recently developed set of process-oriented diagnostics is used to examine the azimuthally averaged wind and thermodynamic structures of the GCM-simulated TCs. Results indicate that the inner-core wind structures of simulated TCs are more strongly constrained by the horizontal resolutions of the models than are the thermodynamic structures of those TCs. As expected, the structures of TC circulations become more realistic with smaller horizontal grid spacing, such that the radii of maximum wind (RMW) become smaller, and the maximum vertical velocities occur off the center. However, the RMWs are still too large, especially at higher intensities, and there are rising motions occurring at the storm centers, inconsistently with observations. The distributions of precipitation, moisture, radiative and surface turbulent heat fluxes around TCs are diverse, even across models with similar horizontal resolutions. At the same horizontal resolution, models that produce greater rainfall in the inner-core regions tend to simulate stronger TCs. Further, when TCs are weak, the radial gradient of net column radiative flux convergence is comparable to that of surface turbulent heat fluxes, emphasizing the importance of cloud-radiative feedbacks during the early developmental phases of TCs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Model Structure and Climate Data Uncertainty in Historical Simulations of the Terrestrial Carbon Cycle (1850–2014)

Abstract The divergence among Earth system models in the terrestrial carbon cycle has prompted interest in how to reduce uncertainty. Previous studies have identified model structural uncertainty arising from process parameterizations and parameter values. The current study highlights the importance of climate forcing in generating carbon cycle uncertainty. We use simulations in which three models (Community Land Model version 4 (CLM4), CLM4.5, CLM5) with substantially different carbon cycles are forced with two climate reconstructions (CRUNCEPv7, Global Soil Wetness Project 3 version 1 (GSWP3v1)) to examine the contributions of model structure and climate to uncertainty in the carbon cycle over the period 1850–2014. Climate uncertainty for global annual net biome production exceeds one third of total uncertainty (defined as the sum of climate and model structure uncertainty) in the first half of the twentieth century, but declines after the 1950s. Global annual gross primary productivity, net primary productivity, heterotrophic respiration, and vegetation and soil carbon stocks have substantial climate uncertainty (relative to total uncertainty) throughout the simulation period. Climate forcing contributes more than one half of total uncertainty for these carbon cycle fluxes and stocks throughout boreal North America and Eurasia, some midlatitude regions, and in eastern Amazonia and western equatorial Africa during the decade 2000–2009. Comparison with observationally based data sets of the carbon cycle using model benchmarking methods provides insight into strengths and deficiencies among models and climate forcings, but we caution against overreliance on benchmarking to discriminate among models. The conceptualization of uncertainty arising from this study implies embracing multiple feasible model simulations rather than focusing on which model or simulation is best.

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

TeraChem protocol buffers ( TCPB ): Accelerating QM and QM/MM simulations with a client–server model

The routine use of electronic structures in many chemical simulation applications calls for efficient and easy ways to access electronic structure programs. Here, we describe how the graphics processing unit (GPU) accelerated electronic structure program TeraChem can be set up as an electronic structure server, to be easily accessed by third-party client programs. We exploit Google’s protocol buffer framework for data serialization and communication. The client interface, called TeraChem protocol buffers (TCPB), has been designed for ease of use and compatibility with multiple programming languages, such as C++, Fortran, and Python. To demonstrate the ease of coupling third-party programs with electronic structures using TCPB, we have incorporated the TCPB client into Amber for quantum mechanics/molecular mechanics (QM/MM) simulations. The TCPB interface saves time with GPU initialization and I/O operations, achieving a speedup of more than 2× compared to a prior file-based implementation for a QM region with ~250 basis functions. We demonstrate the practical application of TCPB by computing the free energy profile of p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI - )—a model chromophore in green fluorescent proteins—on the first excited singlet state using Hamiltonian replica exchange for enhanced sampling. All calculations in this work have been performed with the non-commercial freely-available version of TeraChem, which is sufficient for many QM region sizes in common use.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗