Search NASA⌕ Search

SEARCH · Search NASA

Results for “GROUP DYNAMICS”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

EvoNet: A phylogenomic and systems biology approach to identify genes underlying plant survival in marginal, low‐N soils

The DOE‐BER “EvoNet” project investigates the genetic and molecular basis of plant resilience in extreme environments. We do this by identifying key genes that enable “extreme survivor” species to thrive in the nitrogen-poor soils of Chile’s hyper-arid Atacama Desert. Our collections focus on 32 Atacama extremophile species, including seven grass species with potential biofuel applications. To identify genes-of-importance to survival we compared genomic and transcriptomic profiles of extremophile species that thrive in the Atacama to those of closely related “sister” species from nitrogen-rich arid and mesic regions of California. Deep RNA sequencing and de novo transcriptome assembly across these triplet species sets supported a phylogenomic framework for identifying positively selected genes associated with adaptive divergence. Our integrative analysis combined ecological and environmental data, metagenomics, evolutionary and systems biology, and metabolomics. This enabled us to create an unprecedented framework for systematically understanding how non-model plants have adapted to survive in extreme conditions. Our resulting database of positively selected ortholog groups in the extremophile plants offers promising targets for engineering crop and biofuel species with enhanced resilience to drought and extreme weather. Additionally, our newest dataset explores and exploits a complementary metabolomic approach. This new aspect provides innovative strategies to manipulate plant cell metabolism, further supporting efforts to improve agricultural productivity in the face of extreme climates. Importantly, our combined evolutionary- and metabolomic-based strategies focused on convergent patterns of adaptation, providing a genetic and metabolomic toolkit for improving crop and biofuel resilience across diverse plant species. Finally, our novel exploration of ecological and evolutionary dynamics delivered to the community a phylogenomic computational pipeline called “PhyloGeneious.” Our continued adaptations of this pipeline are publicly available to expedite evolutionary genomic research for future scientific discoveries. In total, our DOE-BER has provided genomic, metabolomic, and computational strategies to understand how extremophile plants provide evolutionary and physiological targets for improving agricultural and biofuel production.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

The Arctic

The Arctic environment in 2024 continued on a trajectory that has put it in a state far different from that of the twentieth century. Ongoing accumulation of greenhouse gases in the atmosphere continues to quickly warm the Arctic, resulting in rapid changes in the cryosphere that are driving cascading impacts to climate, ecological, and societal systems. Many weather- and climate-related impacts in the Arctic are the result of compounding change, such as increased riverbank erosion, which is proximately due to increased river discharge from higher seasonal precipitation, yet is also exacerbated by thawing permafrost. However, even individual storms occur within very different ocean and ice conditions than were typically present in the late twentieth century. As a result, the impacts, including high winds, excessive precipitation, and coastal inundation, may be quite different nowadays, as exemplified by the October 2024 storm in northwest Alaska that produced severe coastal flooding in several communities. To share some of these impacts with a wider audience, select extreme weather impacts around the greater Arctic have been highlighted through the inclusion of sidebars in recent State of the Climate Arctic chapters (e.g., Benestad et al. 2023; Thoman et al. 2024).

Thoman, Richard L. [Univ. of Alaska, Fairbanks, AK↗
DOE OSTI2026

Hidden features in the OH-stretching spectra of amino acid decorated air–water interfaces

Chemical reactivity at the air–water interface is governed by the interfacial solvation of reactive species. For instance, during aqueous amino acid-based CO 2 absorption, water reorganizes around the reactive sites and couples dynamically with reaction pathways, facilitating the reaction. In this context, surface-sensitive vibrational sum-frequency generation (vSFG) spectroscopy can probe the OH stretch vibrations of interfacial water and determine the solvation structures around reactants and products, thereby furthering our understanding of the role of interfacial solvation. However, vSFG spectra of the air–water interface in the presence of charged species can be remarkably complex; key species-bound local water structures with distinct orientations may be hidden beneath prominent vSFG peaks arising from water–water hydrogen bonds and remain difficult to resolve. Here, we measure and compute vSFG spectra of the water OH stretch at air–water interfaces decorated with amino acids in their zwitterionic and anionic forms, as well as equimolar mixtures of these forms with bicarbonate. The latter represents post-CO 2 -absorption conditions. We find that computing depth- and frequency-dependent spectral densities—decomposed into contributions from water molecules hydrogen-bonded exclusively to other water molecules, exclusively to amines, exclusively to carboxylates, or shared between these polar/charged groups—is indispensable for accurate interpretation of the vSFG spectra. Key findings include orientational flip-flop in water sub-layers, strong carboxylate-water H-bonding, and water orientational ordering extending into the bulk aqueous phase induced by anionic amino acids. Here, this study provides a computational spectroscopic platform for improved understanding of interfacial solvation relevant to interfacial reactivity.

Air-water interface↗
DOE OSTI2025

Active Site Dynamics in Molybdenum-Based Silica-Supported Olefin Metathesis Catalysts: Site Renewal and Decay Beyond the Chauvin Cycle

Heterogeneous olefin metathesis catalysts exhibit low active site densities and unpredictable kinetics due to dynamic active site formation and decay processes. Here, in this study, we establish a quantitative framework that captures active site generation, renewal, and decay in olefin metathesis over silica-supported molybdenum oxide catalysts, enabling a mechanistic explanation of catalytic behavior and strategies to achieve high, stable activity. Steady-state active site titrations reveal that 2,3-dimethyl-butene isomers (4MEs) cofeeding increases active site density by up to 4.3-fold, directly correlating with enhanced metathesis rates. Spectroscopic studies demonstrate that 4MEs facilitate Mo(VI) reduction to Mo(IV) and interact strongly with surface Si–OH groups, generating labile protons that promote active site formation via a 1,2-proton shift mechanism. Kinetic modeling indicates that ethylene acts as a decay promoter, shifting kinetic control away from the Chauvin cycle and suppressing metathesis activity. Comparative studies on catalysts with varying Mo loading reveal that promotion is most effective for dispersed molybdate species, with a decline at higher Mo loadings. These findings provide a unified mechanistic framework for heterogeneous olefin metathesis, offering new strategies to enhance active site accessibility, mitigate deactivation, and optimize catalyst design.

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

Anomalies of coset non-invertible symmetries

Anomalies of global symmetries provide important information on the quantum dynamics. We show the dynamical constraints can be organized into three classes: genuine anomalies, fractional topological responses, and integer responses that can be realized in symmetry-protected topological (SPT) phases. Coset symmetry can be present in many physical systems including quantum spin liquids, and the coset symmetry can be a non-invertible symmetry. We introduce twists in coset symmetries, which modify the fusion rules and the generalized Frobenius-Schur indicators. We call such coset symmetries twisted coset symmetries, and they are labeled by the quadruple (G, K, ω D+1 , α D ) in D spacetime dimensions where G is a group and K ⊂ G is a discrete subgroup, ω D+1 is a (D + 1)-cocycle for group G, and α D is a D-cochain for group K. We present several examples with twisted coset symmetries using lattice models and field theory, including both gapped and gapless systems (such as gapless symmetry-protected topological phases). We investigate the anomalies of general twisted coset symmetry, which presents obstructions to realizing the coset symmetry in (gapped) symmetry-protected topological phases. We show that finite coset symmetry G/K becomes anomalous when G cannot be expressed as the bicrossed product G = H ⋈ K , and such anomalous coset symmetry leads to symmetry-enforced gaplessness in generic spacetime dimensions. We illustrate examples of anomalous coset symmetries with A 5 / \Z 2 symmetry, with realizations in lattice models.

Hsin, Po-Shen [King’s College London (United Kingd↗
DOE OSTI2024

Molecular-level design of alternative media for energy-saving pilot-scale fibrillation of nanocellulose

The outstanding mechanical properties, light weight, and biodegradability of cellulose nanofibrils (CNFs) make them promising components of renewable and sustainable next-generation reinforced composite biomaterials and bioplastics. Manufacturing CNFs at a pilot scale requires disc-refining fibrillation of dilute cellulose fibers in aqueous pulp suspensions to shear the fibers apart into their nanodimensional forms, which is, however, an energy-intensive process. Here, in this study, we used atomistic molecular dynamics (MD) simulation to examine media that might facilitate the reduction of interactions between cellulose fibers, thereby reducing energy consumption in fibrillation. The most suitable medium found by the simulations was an aqueous solution with 0.007:0.012 wt.% NaOH:urea, and indeed this was found in pilot-scale experiments to reduce the fibrillation energy by ~21% on average relative to water alone. The NaOH:urea-mediated CNFs have similar crystallinity, morphology, and mechanical strength to those formed in water. The NaOH and urea act synergistically on CNFs to aid fibrillation but at different length scales. NaOH deprotonates hydroxyl groups leading to mesoscale electrostatic repulsion between fibrils, whereas urea forms hydrogen bonds with protonated hydroxyl groups thus disrupting interfibril hydrogen bonds. This suggests a general mechanism in which an aqueous medium that contains a strong base and a small organic molecule acting as a hydrogen-bond acceptor and/or donor may be effectively employed in materials processes where dispersion of deprotonable polymers is required. The study demonstrates how atomic-detail computer simulation can be integrated with pilot-scale experiments in the rational design of materials processes for the circular bioeconomy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Persistent DynAMICS

Persistent DynAMICS is a lightweight software platform that enables distributed sensing systems to operate as a unified, responsive network. It integrates with existing sensors, data systems, and processing tools as an overlay, without requiring hardware replacement or increasing user burden. As data is collected, the system identifies situations that require additional attention and directs selected sensors to gather more information. For example, if a sensor detects an unusual condition along a pipeline, the system can task nearby sensors or different sensing modalities to confirm the event and focus data collection where it is needed. To support this, the platform forms temporary groups of sensors that work together to address specific objectives. Once complete, these groups are dissolved and resources are reassigned as needed. Communication is handled through an integrated messaging framework with client and server components, enabling reliable operation and data continuity even in low-bandwidth or intermittent connectivity environments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Ternary and Quaternary Nickel-Iron-Oxide-Based Thin Films as Oxygen Evolution Reaction Catalysts in Alkaline Media

Alkaline water electrolyzers are a promising technology for energy storage and conversion through the electrochemical production of hydrogen, however the slow kinetics of the oxygen evolution reaction (OER) pose a major challenge. The alkaline environment enables the use of non-Pt-group metals as OER catalysts, in particular NiFe-based materials. In this study, Mn and/or Cr were elementally mixed with NiFe to further improve the OER catalyst properties through a controlled, systematic synthesis process involving thin films prepared by electron beam-induced physical vapor deposition (PVD). By exploring both simultaneous co-deposition and serial deposition of the different metals, the OER performance was investigated as a function of both the location of Mn or Cr in the multi-metal NiFe-based catalysts and the difference between surface versus bulk mixing. The distinct differences in the redox dynamics, surface oxidation, and stability of the seven studied catalysts reveal metal Cr- and/or Mn-overlayers as a promising synthesis approach for fine-tuning of catalyst properties and improving the OER performance.

Alkaline Water Electrolysis↗
DOE OSTI2025

Optomechanical Tuning of Second Harmonic Generation Anisotropy in Janus MoSSe/MoS 2 Heterostructures

Symmetry breaking in van der Waals materials enables the realization of quantum states and advanced device functionalities. Janus transition-metal dichalcogenides (TMDs) exhibit distinctive nonlinear optical properties due to their broken out-of-plane mirror symmetry. However, the dynamic control of second harmonic generation (SHG) anisotropy and resonance behavior via optical excitation remains elusive. Here, in this work, we investigate the SHG response of Janus MoSSe/MoS 2 heterostructures with 2H and 3R stackings. We can tune the SHG response by varying the incident photon wavelength from 800 to 1000 nm, which shows a resonance-dependent enhancement in intensity and a deviation from 6-fold symmetry, indicating wavelength-dependent anisotropy. The ratio between maximum and minimum intensity in the armchair directions, associated with the SHG anisotropy, reaches a value of 1.73 at the excitation wavelength of 1000 nm. Group theory analysis and first-principles calculations reveal that the observed anisotropy arises from optically induced strain. Our findings highlight the role of symmetry breaking and optical resonance contributing to the optomechanical tuning of SHG anisotropy, offering opportunities for developing Janus TMD-based photonic devices for frequency conversion, light generation, and optical switching.

Janus transition metal dichalcogenides↗
DOE OSTI2026

Unravelling chemical pathways of H 2 on Ga 2 O 3 surfaces with spectro-electrochemistry

This work highlights the capability of coupled spectroscopic and electrochemical techniques to probe dynamic surface processes under realistic operating conditions. By simultaneously employing in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and electrochemical impedance spectroscopy (EIS), we elucidate the mechanistic interaction between Ga 2 O 3 and hydrogen under elevated temperatures in a low-oxygen environment. This novel spectro-electrochemical approach allows chemistry to be correlated with the surface charge density of Ga 2 O 3 . Our results reveal a concentration-dependent transition in reaction pathway. At low concentrations, hydrogen reacts with ambient oxygen to form surface hydroxyls. At intermediate concentrations, hydrogen interacts with surface adsorbed oxygen to generate hydroxyl groups along with reducing the surface. Finally, at high H 2 concentrations, hydrogen reduces both hydroxyls and surface oxygen, leading to a highly conductive grain surface. As a result, hydrides form on the reduced Ga 2 O 3 surface. The gained insights are relevant for heterogeneous catalysis and gas sensing.

08 HYDROGEN↗
DOE OSTI2024

Mixing effects on spectroscopy and partonic observables of heavy mesons with logarithmic confining potential in a light-front quark model

Using the variational principle, we systematically investigate the mass spectra and wave functions of both 1⁢𝑆 and 2⁢𝑆 state heavy pseudoscalar (𝑃) and vector (𝑉) mesons within the light-front quark model. This approach incorporates a Coulomb plus logarithmic confinement potential to accurately describe the constituent quark and antiquark dynamics. Additionally, spin hyperfine interactions are introduced perturbatively to compute the masses of pseudoscalar and vector mesons. The present analyses of the 1⁢𝑆 and 2⁢𝑆 states require the consideration of mixing between them to account for empirical constraints. These constraints include the mass gap Δ⁢𝑀 𝑃 >Δ⁢𝑀 𝑉 , where Δ⁢𝑀 𝑃⁡(𝑉) =𝑀$^{2⁢𝑆}_{𝑃⁡(𝑉)}$−𝑀$^{1⁢𝑆}_{𝑃⁡(𝑉)}$ and the hierarchy of the decay constants 𝑓 1⁢𝑆 >𝑓 2⁢𝑆 . We find the optimal value of the mixing angle to be 𝜃 =1⁢8°, significantly enhancing the consistency between our spectroscopic predictions and the experimental data compiled by the Particle Data Group. Furthermore, based on the predicted mass, the newly observed resonance 𝐵 𝐽⁡ (5840) could be assigned as a 2 1⁢ 𝑆 0 state in the 𝐵 meson family. The study also reports various pertinent observables, including twist-two distribution amplitudes, electromagnetic form factors, charge radii, 𝜉 moments, and transition form factors that are found to be consistent with both available lattice simulations and experimental data. In addition, our predicted branching ratios for the channels of 𝐵 + →𝜏 + ⁢𝜈 𝜏 as well as rare decays of 𝐵 0 and 𝐵$^0_𝑠$ appear in accordance with experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Designing Physicochemically‐Ordered Interphases for High‐Performance Composites

To enhance the mechanical properties of carbon fiber‐reinforced polymer composites, a physicochemical scaffold is designed incorporating microscopically architected chemically reactive nanofibers that act as a multiscale bridge between the carbon fibers and the matrix. Thermally activated nanofibers leverage their morphologically driven mechanochemical properties to form covalent bonds with adjacent polymer molecules, creating a co‐continuous network that dramatically enhances fiber‐matrix load transfer. By meticulously controlling the nanofiber architecture through variable surface area, functional group availability, and polymer chain alignment effects, the extent of covalent bonding between nanofibers and the matrix is manipulated ultimately resulting in improved carbon fiber‐matrix adhesion. Further, the concept was validated using polyacrylonitrile nanofibers within an acrylonitrile butadiene styrene matrix in a discontinuous carbon fiber‐reinforced composite system. Nanomechanical studies using atomic force microscopy and low‐field nuclear magnetic resonance spectroscopy confirmed immobilized, chemically transferred, and ordered nanostructures at the interphase. The resulting composites demonstrated ≈56% and ≈175% improvements in tensile strength and toughness, respectively, compared to composites without nanofiber. Comprehensive thermal, rheological, and X‐ray scattering analyses, along side all‐atomic molecular dynamics simulations, revealed the fundamental mechanisms behind these improvements in mechanical behavior. The versatility and efficacy of the approach have the potential to address longstanding interphase challenges in the composite industry.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A BACUS S UMMIT N -body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Predicting the dark matter-baryon abundance ratio

We discuss relaxation solutions to the dark matter-baryon coincidence problem in the context of QCD axion dark matter. In relaxation solutions, a moduli dynamically adjusts the mass of dark matter and baryons until their energy densities are 𝒪⁡(1) the same. Because the QCD axion is heavily connected to QCD, scanning the QCD axion mass inherently also scans the proton mass. In the context of relaxation solutions, this implies that the ratio of dark matter to baryon abundances (Ω DM /Ω B ) is a ratio of beta functions showing that these models can only accommodate discrete values of ΩDM/ΩB thereby “predicting” the ratio of the dark matter to baryon abundances. The original composite axion model has only a single integer degree of freedom 𝑁, the size of the gauge group, and we show that when 𝑁 =8 the observed value of Ω DM /Ω B =5.36 is reproduced to within its percent level error bars. Novel tests of this model include more precise measurements of Ω DM /Ω B , a better lattice determination of the dependence of the proton mass on the high energy QCD gauge coupling, as well as more traditional tests such as fifth force experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Hierarchical semi-Markov models with duration-aware dynamics for activity sequences

Residential electricity demand at granular scales is driven by what people do and for how long. Accurately forecasting this demand for applications like microgrid management and demand response therefore requires generative models for activities that can produce realistic daily activity sequences, capturing both the timing and duration of human behavior. This paper develops a generative model of human activity sequences using nationally representative time-use diaries at a 10-min resolution. We use this model to quantify which demographic factors are most critical for improving predictive performance. We propose a hierarchical semi-Markov framework that addresses two key modeling challenges. First, a time-inhomogeneous Markov router learns the patterns of “which activity comes next.” Second, a semi-Markov hazard component explicitly models activity durations, capturing “how long” activities realistically last. To ensure statistical stability when data are sparse, the model pools information across related demographic groups and time blocks. The entire framework is trained and evaluated using survey design weights to ensure our findings are representative of the U.S. population. On a held-out test set, we demonstrate that explicitly modeling durations with the hazard component provides a substantial and statistically significant improvement over purely Markovian models. Furthermore, our analysis reveals a clear hierarchy of demographic factors: Sex, Day-Type, and Household Size provide the largest predictive gains, while Region and Season, though important for energy calculations, contribute little to predicting the activity sequence itself. The result is an interpretable and robust generator of synthetic activity traces, providing a high-fidelity foundation for downstream energy systems modeling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Thermo-mechanical deterioration and molecular degradation of 3D-printed methacrylate-based polymer in various chemical environments

The advancement of additive manufacturing (AM) has accelerated the development of stereolithographic (SLA) photo-curable resins, particularly methacrylate-based polymers, due to the ability to their high-resolution, robust mechanical properties, and suitability. However, their long-term performance in chemical environments remains poorly understood. This study investigates the extent and mechanisms of degradation on an SLA-printed methacrylate-based polymer subjected to various chemicals, including polar and non-polar solvents, as well as strongly acidic aqueous solutions over a 12-week accelerated aging period. A comprehensive analytical approach incorporating swelling kinetics, surface morphology, tensile and dynamic mechanical analysis (DMA), Fourier-transform infrared (FTIR) spectroscopy, and mass spectrometry (MS) was employed to characterize chemical absorption, structural integrity, and leached products. Results reveal that degradation severity is governed by both the polarity and reactivity of the chemical environment. Notably, exposure to 6 mol L -1 HNO₃ induced the most severe deterioration, with over threefold higher swelling compared to other media, and significant reductions in tensile strength, tensile modulus, and glass transition temperature (T g ). In contrast, specimens aged in non-polar solvents (xylene and dodecane) exhibited negligible chemical interaction and retained mechanical performance. FTIR and MS analyses identified acid-catalyzed hydrolysis of ester groups as prominent degradation pathways in acidic media, while diffusion-controlled plasticization prevailed in polar solvents. Furthermore, this study provides valuable insights into the chemical stability of SLA-printed polymers and develops predictive degradation profiles that are crucial for designing durable polymer systems for advanced industrial use.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Prairie State Generating Company Site: Storage Field Development Plan

This Storage Field Development plan presents the Storage Complex characterization results, construction, monitoring, operational plans, and costs associated with the proposed Prairie State Energy Campus - Carbon Capture and Storage (PSEC-CCS) project in Washington County, Illinois, near Marissa. The proposed storage complex, known as the St. Peter-Everton Storage Complex, comprises the St. Peter Sandstone and Everton sandstone reservoirs, with the Maquoketa Shale serving as the primary confining unit and the Knox Group as lower confining unit. The lowermost Underground Source of Drinking Water (USDW) identified for the site is the Pennsylvanian section.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Best practices in software development for robust and reproducible geoscientific models based on insights from the Global Carbon Budget's dynamic vegetation models

Computational models play an increasingly vital role in scientific research by enabling the numerical simulation of complex processes. Such models are also fundamental in geosciences. For instance, they offer critical insights into the impacts of global change on the Earth system today and in the future. Beyond their value as research tools, models are also software products and should therefore adhere to certain established software engineering standards. However, scientists are rarely trained as software developers, which can lead to potential deficiencies in software quality like unreadable, inefficient, or erroneous code. The complexity of models, coupled with their integration into broader workflows, also often makes it challenging to reproduce results, evaluate processes, and build upon them. In this paper, we review the state and current practices of the development processes of the state-of-the-art land surface models used by the Global Carbon Budget. We combine the experience of modelers from the respective research groups with the expertise of software engineers from tech companies to outline key principles and tools for improving software quality in research. We explore four main areas: (1) model testing and validation, (2) scientific, technical, and user documentation, (3) version control, continuous integration, and code review, and (4) the portability and reproducibility of workflows. Our review reveals that while modeling communities are incorporating many best practices, significant room for improvement remains in areas such as automated testing, automated documentation, and reproducibility. Therefore, we here identify and promote essential software engineering practices, including numerous examples of practices from within the community that can serve as guidelines for other models and could help streamline processes across the entire community. We conclude with an open-source example implementation of these principles, demonstrating portable and reproducible data flows, a continuous integration setup, and web-based visualizations. This example may serve as a practical resource for model developers, users, and all scientists engaged in scientific programming.

Gregor, Konstantin [Technical Univ. of Munich (Ger↗