Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mass diffusivity”

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 127 records · Page 7

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Soil porous microstructure control over soil organic matter mobility: A multimethod workflow for understanding chemistry-dependent organic matter binding in soil

Soil organic matter (SOM) has attracted a great deal of interest; particularly for its potential to mitigate human derived CO 2 emissions. Studies have demonstrated that SOM plays a critical role in carbon storage and CO 2 sequestration. However, the sorption properties of SOM, which influence its transport in pore water and stabilization within the soil, remain poorly understood. This study develops a workflow to: (1) examine compound-specific advective and diffusive transport and desorption behaviors, (2) quantify desorption rates through stop-flow and continuous-flow column experiments, and (3) evaluate the impact of soil microporosity on SOM mobility using high-resolution imaging and extractions. Intact core column experiments were conducted on Uncultivated (Natural) and Cultivated soil samples, both were arid soils, collected in Washington State. X-ray computed tomography was employed to measure porosity and pore connectivity, while Fourier-transform ion cyclotron resonance mass spectrometry was used to analyze SOM composition. The findings revealed that cultivation increased total carbon and nitrogen levels due to irrigation and fertilization, enhancing carbon capture potential in arid soils. In contrast, the Natural soil, characterized by higher porosity and connectivity, contained more oxidized carbon. Pore network analysis indicated that soil compaction in the Cultivated soil may lead to longer diffusion pathways, significantly influencing SOM transport and stability.

Hydraulic Properties↗

Porous transport electrodes for oxygen evolution reaction in proton exchange membrane water electrolysis -cells: Materials, designs, and diagnoses

H 2 production using proton exchange membrane (PEM) water electrolysis (PEMWE) cells has received considerable attention because of the high efficiencies of these cells and no harmful emissions from the related process. In PEMWE cells, porous transport electrodes (PTEs) composed of a catalyst layer (CL) comprising O 2 evolution reaction (OER) catalysts, porous transport layer (PTL), and PEM play key roles in the stack performance and lifetime. Herein, Ir-based and non-precious-metal OER catalysts that are highly active and stable at low pH values and high anodic potentials are reviewed to understand their OER mechanisms. Various strategies are proposed for engineering CLs and PTLs to improve the interfacial properties and mass transfers of reactants and products to and from the active sites. Additionally, diagnoses of PTEs is significantly crucial for interpreting electrochemical processes and addressing their current challenges. Therefore, half-cell analyses, including diffusion electrode (DE), floating electrode (FE), and modified rotating disk electrode (MRDE) techniques, are explored, and membrane electrode assembly (MEA)-based analyses, such as the polarization technique, electrochemical impedance spectroscopy, and magnetic field analysis, are established. In conclusion, this study aims to provide an overview of recent technologies used for the engineering and diagnostic tools of PEMWE cells and insights into the advanced components and systems to be developed in this field.

Diagnosis of PEMWE Cells↗

Tailoring the Physicochemical Properties of Nb Thin Films via Surface Engineering Methods

The modification of surface oxide layers formed on niobium (Nb) thin films via chemical mechanical planarization (CMP) and accelerated neutral atom beam (ANAB) processing provides a promising route toward tailoring their emergent properties and performance when used as superconducting qubits. Here, in this study, we show that CMP- and ANAB-formed Nb oxides are significantly thinner and smoother than the native oxide, as revealed by transmission electron microscopy (TEM) and atomic force microscopy. Scanning TEM and energy-dispersive X-ray spectroscopy along with X-ray photoelectron spectroscopy identified an oxidation gradient within the native and surface-engineered oxides. The topside layer is dominated by Nb 5+ (Nb 2 O 5 ), with various Nb suboxides present closer to the oxide/metal interface. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) depth profiling confirmed the presence of an oxygen content gradient and demonstrated the enhanced resistance of the CMP- and ANAB-formed oxides to oxygen surface exchange and subsequent diffusion via 18 O 2 isotopic labeling experiments. ToF-SIMS also identified an interfacial layer containing trapped hydrogen (H)-containing species at the Nb oxide/metal interface. In situ ToF-SIMS and TEM revealed migration of the H/OH interfacial layer coinciding with decomposition of the surface oxide. Furthermore, our density functional theory calculations indicated that both H from moisture present in ambient air and bulk H in Nb films tend to segregate at the interface. These findings underscore the importance of understanding surface oxidation mechanisms, hydrogen incorporation, and their impact on the designed functionalities of Nb-based devices.

36 MATERIALS SCIENCE↗

The generation of a multiphase medium in ‘Splash’ bridge systems: towards an understanding of star formation suppression in turbulent galaxy systems

ABSTRACT Cloud–cloud collisions in splash bridges produced in gas-rich disc galaxy collisions offer a brief but interesting environment to study the effects of shocks and turbulence on star formation rates in the diffuse intergalactic medium, far from the significant feedback effects of massive star formation and active galactic nucleus. Expanding on our earlier work, we describe simulated collisions between counter-rotating disc galaxies of relatively similar mass, focusing on the thermal and kinematic effects of relative inclination and disc offset at the closest approach. This includes essential heating and cooling signatures, which go some way towards explaining the luminous power in H$_2$ and [C ii] emission in the Taffy bridge, as well as providing a partial explanation of the turbulent nature of the recently observed compact CO-emitting clouds observed in Taffy by the Atacama Large Millimeter Array (ALMA). The models show counter-rotating disc collisions result in swirling, shearing kinematics for the gas in much of the post-collision bridge. Gas with little specific angular momentum due to collisions between counter-rotating streams accumulates near the centre of mass. The disturbances and mixing in the bridge drive continuing cloud collisions, differential shock heating, and cooling throughout. A wide range of relative gas phases and line-of-sight velocity distributions are found in the bridges, depending sensitively on initial disc orientations, and the resulting variety of cloud collision histories. Most cloud collisions can occur promptly or persist for quite a long duration. Cold and hot phases can largely overlap throughout the bridge or can be separated into different parts of the bridge.

79 ASTRONOMY AND ASTROPHYSICS↗

Active Sites in the Dealuminated Beta Zeolite-Supported Cobalt Catalyst for Non-Oxidative Ethane Dehydrogenation

Dispersed metal species in siliceous zeolites have been actively studied for non-oxidative dehydrogenation of ethane (NDE). Fundamental insights into the dynamics of metal species in zeolites under reaction conditions have rarely been explored. Herein, we report an atomic level understanding of the dynamics and activity of cobalt (Co) sites in dealuminated Beta zeolite (DeAl-BEA) for NDE during induction and reaction conditions with extensive characterization techniques such as diffuse reflectance UV–vis, solid state nuclear magnetic resonance and X-ray photoelectron, X-ray diffraction along with in situ Fourier transform infrared and X-ray absorption spectroscopy. For a catalyst with 0.5 mass % Co loading, tetrahedral Co 2+ mononuclear sites, di-coordinated to the zeolite framework and with two silanol groups in vicinity (i.e., (≡SiO) 2 Co(HO–Si≡) 2 ), form upon exposure to hydrogen during induction and persist through the NDE reaction. Increasing the Co loading to 3.0 mass % yielded Co sites with similar electronic and coordination structures but slightly elongated Co–O bonds. Upon cooling to room temperature, the Co sites persisted in the same coordination environment, though the disappearance of a feature in the Co K-edge near-edge region revealed changes in the active site’s electronic structure coinciding with modest shifts in bond lengths. The electronic structure and activity of (≡SiO) 2 Co(HO–Si≡) 2 sites were studied comparatively to a few other hypothetical Co 2+ coordination structures, using electronic structure calculations and microkinetic simulations. The simulations showed that NDE is controlled by β-hydride elimination following C–H bond activation and that Co-sites possessing flexibility because of neighboring silanol defects are more active. Interestingly, dinuclear Co–O–Co sites (i.e., (≡SiO)Co(HO–Si≡) 2 –O–(HO–Si≡) 2 Co(≡SiO)) were more active than the mononuclear (≡SiO) 2 Co(HO–Si≡) 2 sites because of favorable hydrogen bonding with the vicinal silanol groups. In conclusion, the present study bridges the gap between the knowledge acquired by ex-situ characterizations and the active sites under the reaction conditions in alkane dehydrogenation chemistry.

36 MATERIALS SCIENCE↗

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗

Interpenetrating 3D Electrodes for High-Rate Alkaline Water Splitting

Mass transfer is critical for the reaction kinetics and efficiency of alkaline water splitting (AWS). For AWS to operate at high current densities (hundreds of mA/cm 2 ), the device architecture must ensure a large catalytic surface area, rapid ion diffusion, and minimal solution and charge transfer resistances. Effective electrodes should also facilitate gas bubble detachment and release. 3D-printed electrodes have shown promise, but stacking them increases the ion diffusion length and solution resistance. Here we demonstrate a new device architecture with interpenetrating gyroid electrodes, providing a large ion-accessible surface area and gas diffusion channels. This design significantly reduces the interelectrode distance, lowering ion diffusion length and solution resistance. Simulations show faster ion diffusion and higher current density in the interpenetrating configuration compared with separate electrodes. This improved performance, especially at low temperatures and high current densities, highlights a promising strategy for enhancing AWS and other electrochemical systems limited by slow ion diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Translating Fundamental Insights into Ag-Based Bimetallic Electrocatalysts to Anion-Exchange Membrane Fuel Cells

To address challenges of high costs and scalability for fuel cells, it is essential to develop af-fordable and earth-abundant materials that are Pt-group-metal (PGM)-free. Recent advance-ments in PGM-free electrocatalysis for the oxygen reduction reaction (ORR) in alkaline media show that high catalyst loadings are needed to achieve high reaction rates; however, there are associated mass transport limitations. To address this issue, we develop low-loading ionomer-less Ag-bimetallic thin films by alloying with 3d block elements or Sn. We bridge knowledge from fundamental rotating disk electrode (RDE) studies to gas-diffusion cathodes in high-temperature anion-exchange membrane fuel cells (HT-AEMFCs), resulting in high-activity ORR catalysis and peak power densities up to 1.2 W cm-² geo for the Ag-Sn and Ag-Co alloys. Here, experimental post-characterization and theoretical calculations reveal small Co or Sn oxide nano-islands on Ag as likely active sites enhancing ORR activity, ultimately offering these low-loading PGM-free ORR materials as a viable path towards sustainable energy conversion.

Catalysts↗

Decoupling of Nitrogen and Oxygen Impurities in Nitrogen Doped SRF Cavities

The performance of superconducting radiofrequency (SRF) cavities is critical to enabling the next generation of efficient high-energy particle accelerators. Recent developments have focused on altering the surface impurity profile through \textit{in-situ} baking, furnace baking, and doping to introduce and diffuse beneficial impurities such as nitrogen, oxygen, and carbon. However, the precise role and properties of each impurity are not well understood. In this work, we attempt to disentangle the role of nitrogen and oxygen impurities through time-of-flight secondary ion mass spectrometry of niobium samples baked at temperatures varying from 75-800 $^\circ$C with and without nitrogen injection. From these results, we developed treatments recipe that decouple the effects of oxygen and nitrogen in doping treatments. Understanding how these impurities and their underlying mechanisms drive further optimization in the tailoring of impurity profiles for high-performance SRF cavities.

43 PARTICLE ACCELERATORS↗

Decoupling of Nitrogen and Oxygen Impurities in Nitrogen Doped SRF Cavities

The performance of superconducting radiofrequency (SRF) cavities is critical for enabling the next generation of efficient high-energy particle accelerators. Recent developments have focused on altering the surface impurity profile through \textit{in-situ} baking, furnace baking, and doping to introduce and diffuse beneficial impurities such as nitrogen, oxygen, and carbon. However, the precise role and properties of each impurity are not well understood. In this work, we attempt to disentangle the role of nitrogen and oxygen impurities through time-of-flight secondary ion mass spectrometry using niobium cavity cutouts baked at temperatures varying from 120-800$^\circ$C with and without nitrogen injection. Understanding the underlying mechanisms behind these impurities will drive further optimization in the tailoring of nitrogen and oxygen profiles for high-performance SRF cavities.

43 PARTICLE ACCELERATORS↗

Impact of Site Identity, Location, and Accessibility on Polyethylene Conversion Rates and Product Selectivities over Metal-Free MFI Zeolites

Brønsted (BAS), Lewis (LAS), and surface Brønsted (SBAS) acid sites have been investigated for polyethylene (PE) upcycling by zeolite catalysts, but there is no clear consensus regarding their catalytic roles, partly due to the complexity of the catalysts used and varying reaction conditions across studies. This work systematically determined how these sites impact PE conversion rates and product distributions by utilizing a suite of microporous MFI catalysts with varying Si/Al ratios, acid site densities, and inherent mesoporosities. PE conversion rates did not trend with total BAS or LAS densities due to a combination of internal mass transfer limitation and the apparent inability of LAS alone to cleave C–C bonds, but a strong, statistically significant correlation was present with respect to SBAS density and mesopore surface area, jointly, owing to accelerated polymer activation on external surfaces to smaller diffusion-limited chains. However, ingress of these SBAS-derived fragments ultimately remained rate limiting, as demonstrated by solid conversion rates that increased with mesopore surface area at similar SBAS density and likewise increased with SBAS density at similar mesopore surface area. In batch PE cracking reactions, light gaseous product selectivities were most sensitive to total BAS, with higher densities generally exhibiting higher selectivity to C 3 and linear C 4 –C 7 products and higher alkane/alkene product ratios, consistent with increased β-scission turnovers. Insights from this work help systematically clarify the roles of BAS, LAS, SBAS, and mesopores in PE cracking reactions and inform the development of tailored zeolite catalysts for efficient polyolefin upcycling.

accessibility↗

In Vitro Antimicrobial Activity of Volatile Compounds from the Lichen Pseudevernia furfuracea (L.) Zopf. Against Multidrug-Resistant Bacteria and Fish Pathogens

Lichens are symbiotic organisms with unique secondary metabolism. Various metabolites from lichens have shown antimicrobial activity. Nevertheless, very few studies have investigated the antimicrobial potential of the volatile compounds they produce. This study investigates the chemical composition and antimicrobial properties of volatile compounds from Pseudevernia furfuracea collected in two regions of Morocco. Hydrodistillation was used to obtain volatile compounds from samples collected in the High Atlas and Middle Atlas. Gas chromatography–mass spectrometry (GC-MS) analysis identified phenolic cyclic compounds as the primary constituents, with atraric acid and chloroatranol being the most abundant. Additionally, eight compounds were detected in lichens for the first time. The antimicrobial activity of these compounds was assessed using disc diffusion and broth microdilution methods. Both samples demonstrated significant antimicrobial effects against multidrug-resistant human bacteria, reference microorganisms, fish pathogens, and Candida albicans, with minimum inhibitory concentrations (MICs) ranging from 1000 µg/mL to 31.25 µg/mL. This study provides the first report on the volatile compounds from Pseudevernia furfuracea and their antimicrobial effects, particularly against fish pathogens, suggesting their potential as novel antimicrobial agents for human and veterinary use. Further research is warranted to explore these findings in more detail.

Essadki, Yasser (ORCID:0009000648460075)↗

Estimation of Fission Product Transport Parameters for Cesium in the AGR-3/4 TRISO Fuel Experiment

A one-dimensional (1D) finite-element model of fission product transport in the AGR-3/4 experiment has been developed using the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and implemented in the fuel performance code, BISON. The model resolves capsule-specific geometries, materials, and temperature histories and simulates radial migration of fission products from the fuel compact through the inner ring, outer ring, and into the sink ring. Model parameters governing diffusion and sorption were estimated for key fission products – cesium (Cs), and europium (Eu) – by simultaneously fitting modeled isotopic concentration profiles and total ring inventories to a post-irradiation experimental measurement. These data include gamma scanning, liquid scintillation for Sr-90, radial deconsolidation leach-burn-leach analysis, tomographic reconstructions, and destructive physical sampling. A mortar-based interfacial sorption framework was implemented to enforce physically consistent mass transfer and flux conservation across gas gaps. Two classes of parameter sets were derived: a least-squares best-fit, and a safety-oriented conservative-fit, what applies strong penalties for underprediction of sink inventories. Across all twelve capsules, the model successfully reproduces the dominant radial transport trends for Cs, Sr, with decreasing concentrations from the compact outward through successive rings. Cs behavior is captured most consistently, while strontium predictions reveal systematic trade-offs between compact accuracy and conservative sink-ring bounding. The results demonstrate that sink ring weighted calibration provides conservative, safety relevant bounds on low temperature fission product transport, but at the cost of underpredicting compact inventories for Sr isotopes. These discrepancies highlight the need for additional physics, including fast-slow diffusion model, incorporating trapping mechanism in the transport behavior. Overall, this work establishes a robust, capsule-specific modeling framework for AGR-3/4 fission product transport and provides a defensible basis for parameter selection in source-term and fuel performance analyses for high temperature gas-cooled reactors.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗

Oxygen vacancies in niobium pentoxide as a source of two-level system losses in superconducting niobium

Two-level systems (TLS) have long been a catch-all explanation for RF loss and quantum decoherence in superconducting devices. In our study, the first to directly link TLS losses to a specific physical mechanism, we demonstrate that oxygen vacancies in the naturally formed Nb₂O₅ on oxidized niobium are a major driver of such dissipation. We performed sequential in situ vacuum-baking treatments on niobium superconducting radio-frequency (SRF) cavities and used time-of-flight secondary ion mass spectrometry (TOF-SIMS) to reveal a nonmonotonic evolution in cavity quality factor (Q₀). This behavior correlates with the interplay of Nb₂O₅ vacancy generation and oxide-thickness dissolution. We localize this effect to the oxide itself and present the insignificant role of diffused interstitial oxygen in the underlying Nb by regrowing the oxide via wet oxidation, revealing a mitigation of aggravated TLS losses. We hypothesize that such vacancies in the pentoxide serve as magnetic impurities and are a source of TLS-driven rf loss. Although our measurements center on 3-D SRF cavities, the insights gained here have significant implications for mitigating decoherence in 2-D superconducting qubits.

Bafia, Daniel P. [Fermilab]↗

The impact of environment on size: Galaxies are 50% smaller in the Fornax Cluster compared to the field

Size is a fundamental parameter for measuring the growth of galaxies and the role of the environment on their evolution. However, the conventional size definitions used for this purpose are often biased and miss the diffuse, outermost signatures of galaxy growth, including star formation and gas accretion. We address this issue by examining low surface brightness truncations or galaxy ‘edges’ as a physically motivated tracer of size based on star formation thresholds. Our total sample consists of ∼900 galaxies with stellar masses ranging from 10 5 M ⊙ < M ⋆ < 10 11 M ⊙ . This sample of nearby cluster, group satellite, and nearly isolated field galaxies was compiled using multi-band imaging from the Fornax Deep Survey, deep IAC Stripe 82, and Dark Energy Camera Legacy Surveys. We find that the edge radii scale as R edge ∝ M ⋆ 0.42 , with a very small intrinsic scatter (∼0.07 dex). The scatter is driven by the morphology and environment of galaxies. In both the cluster and field, early-type dwarfs are systematically smaller by approximately 20% compared to late-type dwarfs. However, galaxies in the Fornax cluster are the most impacted. At a fixed stellar mass, edges in the cluster can be found at about 50% smaller radii, and the average stellar surface density at the edges is a factor of two higher, ∼1 M ⊙ /pc 2 . Our findings support the rapid removal of loosely bound neutral hydrogen (H I ) in hot, crowded environments, which truncates galaxies outside-in earlier, preventing the formation of more extended sizes and lower density edges. Our results highlight the importance of deep imaging surveys to the study of low surface brightness imprints of the large-scale structure and environment on galaxy evolution.

Astronomy & Astrophysics↗

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗