The research team, which had previously simulated the effects of static magnetic fields on simple inorganic molecules, alkanols, carbohydrates and porphyrins, turned this time to the lipids that matter most for life. They examined seven eighteen-carbon fatty acids: stearic acid, the fully saturated straight chain; oleic, linoleic and linolenic acids, the cis-unsaturated acids common in plant oils; elaidic and vaccenic acids, the trans-unsaturated acids; and trans-linolenic acid, the fully trans polyunsaturated variant. Alongside these they studied five caproyl glycerides, from monoglycerides through diglycerides to the triglyceride that represents the major form of energy storage in animals and plants. These molecules co-constitute biological membranes, cooperate in skin elasticity and dominate the functional properties of foodstuffs, making them a rational focus for any assessment of magnetic risk.

The simulations were performed entirely in silico, in what the authors call the computer vacuum, using advanced computational methods. Molecular structures were drawn with the Fujitsu SCIGRESS 2.0 software, oriented with their principal symmetry axes along the x-axis of a Cartesian system and enclosed in a three-axial ellipsoid. The molecules were then optimized with Gaussian 0.9 software equipped with the 6-31G** basis set, which computed bond lengths, dipole moments, heats of formation, bond energies and total energies. In the decisive step, the influence of static magnetic fields of 0.0, 0.1, 1, 10 and 100 AMFU, an Arbitrary Magnetic Field Unit in which one AMFU exceeds 1000 tesla, was computed with the Amsterdam Modelling Suite software and the non-relativistic orbital momentum L-dot-B method. The field was fixed along the x-axis with its south pole on the left side of the molecule.

The headline finding is that static magnetic fields polarize lipid molecules in a way that depends strongly on the applied flux density, and that at low strengths the fields neither ionize the molecules nor break their valence bonds. At 0.1 and 1 AMFU the molecules survive intact, but their heats of formation rise steadily, evidence of progressive destabilization. Stearic acid, for example, climbs from a heat of formation of minus 560 kilojoules per mole in the field-free state to minus 511 at 1 AMFU, while its dipole moment grows from 3.42 to 3.84 debye. Every lipid acid and glyceride studied showed the same pattern of increasing heat of formation and increasing dipole moment with rising field strength, confirming that the molecules become less stable and more polar as the magnetic field intensifies.

What struck the researchers most was the irregularity of these changes. The orders in which the molecules destabilized bore no relationship to the number of double bonds or the conformation of the carbon chain. Among the fatty acids, trans-linolenic acid destabilized most, followed by linoleic, linolenic, elaidic, vaccenic, stearic and oleic, while the dipole moments declined in a completely different order led by vaccenic acid. Among the glycerides, the triglyceride destabilized most but the dipole moment sequence was again independent of the degree and position of esterification. The explanation, the authors argue, lies in the ability of these flexible molecules to twist out of the plane initially established by the magnetic field and to squeeze their fragments around certain bonds, creating through-space interactions between atoms that involve van der Waals bonding and dispersion forces.

At the highest simulated strengths the situation turns destructive. Flux densities of 10 and 100 AMFU broke some carbon-hydrogen and carbon-carbon bonds in several of the molecules under consideration. In stearic acid, bonds such as C10-C11, C17-H54, C2-H24 and C1-H20 elongated well above 2 angstroms, the calculated equivalent of ceasing to exist as chemical bonds. Oleic acid survived exposure to 10 AMFU but fractured at 100 AMFU, splitting nine separate bonds. Linolenic acid lost the C19-H49 bond at 10 AMFU and the C1-H39 bond at 100 AMFU, while vaccenic and trans-linolenic acids deteriorated only at the highest strength. When a valence bond disintegrates in this way, the binding electron pair splits into two unpaired electrons and a pair of highly chemically active radicals is generated, species that can bind to ambient molecules and change their chemical structure and biological activity.

Only one molecule in the entire study underwent the cis-trans conversion that magnetic fields are known to provoke. Elaidic acid, the trans-monounsaturated acid, transformed from its trans-conformation into the cis-conformer when exposed to fields above 1 AMFU. The authors note that this ability of static magnetic fields to accelerate the transition of a conformer from cis to trans, or in this singular case the reverse, can significantly affect the oxidation of unsaturated fatty acids, because the auxiliary enzyme enoyl-CoA isomerase that participates in the beta-oxidation pathway works precisely by converting cis-conformers into trans. Vaccenic acid, whose double bond sits at a different position along the chain, resisted any such isomerization even at 100 AMFU, demonstrating how sensitive the outcome is to molecular structure.

The metabolic consequences of these distortions depend on which atoms the field affects most. The atoms of the carboxylic group proved the most susceptible to the magnetic field in every fatty acid studied. The first step of fatty acid metabolism, in which coenzyme A bonds to the carbonyl carbon atom, is favoured by a high positive charge density on that atom and a highly polarized carbon-hydroxyl bond. The simulations show that the field shifts these charge densities irregularly: in stearic acid a field of 0.1 AMFU weakly stimulates the first metabolic step while 1 AMFU inhibits it; in oleic acid fields from 0.1 to 10 AMFU inhibit the first step; in trans-linolenic acid 0.1 AMFU strongly stimulates the first step while 10 AMFU inhibits it. For the glycerides, whose metabolism usually involves hydrolysis by lipases, the charge density on the carboxylic carbon atom varied irregularly with field strength, and only in the triglyceride did the field, regardless of its value, always stimulate hydrolysis.

The authors identify two physical factors behind the susceptibility of molecules to static magnetic fields. The first is the Lorentz force, which acts on moving electrons and at high intensity influences the natural geometry of orbitals, an effect too difficult to include fully in the calculations. The second is the ceasing of the coherence of the binding electron pair, which persists despite electrostatic repulsion by magnetic interactions. A very strong external magnetic field competes with the mutual fields of the paired electrons, gradually weakening their pairing until the bond expands and disintegrates on exceeding a critical length. Because the binding electron pair is the fundamental element of the chemical bond, this gradual splitting is what ultimately generates the radicals observed at high field strengths.

The bottom line of the study is a nuanced risk picture rather than a simple warning. Depending on molecular structure and applied flux density, static magnetic fields either stimulate or inhibit the metabolic processes of the lipids under study, and the direction of the effect depends irregularly on field strength. At the low flux densities typical of everyday exposure, the simulated fields distort and polarize lipid molecules without breaking them, shifting the charge densities that govern enzymatic reactions in ways that could speed up or slow down metabolism. At the extreme strengths of 10 AMFU and above, equivalent to many thousands of tesla, the fields begin to snap valence bonds and generate free radicals. As magnetic pollution of the environment grows and the food industry increasingly exploits static magnetic fields to build functional properties into foodstuffs, the authors argue that understanding these molecular interactions is an essential step toward recognizing the potential risk that such fields pose to flora and fauna.

Subject of Research: Biophysics

Article Title: Potential risk resulting from the influence of static magnetic field upon living organisms. Numerically simulated effects of the static magnetic field upon fatty acids and their glycerides

Article References: Ciesielski, W., Kołoczek, H., Oszczęda, Z., Soroka, J. A., & Tomasik, P. (2023). Potential risk resulting from the influence of static magnetic field upon living organisms. Numerically simulated effects of the static magnetic field upon fatty acids and their glycerides. BioRisk, 19, 1-24. https://doi.org/10.3897/biorisk.19.96250

Image Credits: AI Generated

DOI: 10.3897/biorisk.19.96250

Keywords: static magnetic field, fatty acids, glycerides, lipids, molecular simulation, computational chemistry, free radicals, BioRisk, Potential, risk, resulting, influence

The central insight of the new work is deceptively simple. Many types of materials are currently being considered for use in quantum technology, but most can only hold a stable quantum state when chilled to extreme temperatures. Diamond-based materials offer a compelling alternative because they remain functional at room temperature, removing one of the most burdensome engineering requirements in the field. Pure diamond, however, which is made entirely of carbon atoms, will not work on its own. Defects in the diamond are essential for storing and processing information. A defect might mean a carbon atom is missing from its usual spot, a different kind of atom has taken its place, or both. The clusters of atoms closest to a defect form the fundamental information processing unit of a quantum computer or sensor, known as a qubit.

Simulating those qubits has traditionally been a computational nightmare. Because a defect sits inside an extended crystal lattice, faithful simulations have required modeling enormous numbers of atoms, pushing even powerful computers to their limits and making the systematic screening of candidate defects impractical. Mark Martirez, a staff research scientist at PPPL, has solved this problem by carving out a small cluster of just 16 to 30 atoms around the defect. By treating this cluster as an isolated molecule, he can apply highly accurate quantum chemical tools that chemists have been using for decades to study simple compounds, importing a mature computational arsenal into a domain that has traditionally been approached as a problem in solid-state physics.

The payoff is dramatic. Because the qubits are treated as molecules, the cost of calculating these defects becomes equivalent to the cost of modeling a simple molecule, according to Martirez, who described the method as one that really opened the floodgates in terms of how qubits can be assessed. Using the new approach, the researchers were able to simulate how qubits would respond to light and microwave pulses far faster than ever before. That acceleration matters because reading information out of a qubit is a dynamic process: scientists shine laser light on qubits and pulse microwaves at them to read back the information stored within, and predicting those responses accurately is essential to knowing whether a candidate defect will ever be useful in a real device.

The simulations led the team to discover new boron and oxygen-based defects that could make useful qubits, expanding the palette of options beyond the defect that has dominated the field. Most quantum diamond research has focused on the nitrogen-vacancy defect, in which one carbon atom is replaced by a nitrogen atom and a nearby spot, normally holding another carbon atom, is left empty. When scientists shine a laser on this particular kind of qubit, it emits a red light if it is in a particular quantum state. The oxygen-based defects Martirez is studying shine blue, while the boron-based defects emit ultraviolet light. These different color centers could potentially be used in different ways within a single quantum device, opening architectural possibilities that a single defect type cannot provide.

The practical significance of multiple colors is easier to grasp with a concrete picture. One could imagine combining different qubits with different optical responses, so that a researcher can probe one qubit without disturbing the other, Martirez explained. In a multi-qubit device, that selectivity could simplify control and readout, two of the most delicate operations in quantum engineering. Rather than relying on spatial addressing alone, engineers could in principle use wavelength to distinguish qubits, addressing the blue-emitting centers and the ultraviolet-emitting centers independently. Such flexibility could prove valuable as researchers move from isolated qubits toward the densely integrated systems that practical quantum computers and sensors will demand.

Beyond the immediate candidates, fully understanding the chemistry of qubits could transform the very concept of what a qubit structure should be. Scientists have often treated diamond-based qubits as a physics problem, focusing on the quantum mechanical behavior of spin states within a crystal. The new research shows that they can also be understood through chemistry, enabling scientists to design and build qubits the way a chemist designs custom molecules, such as a new medication, for a specific purpose. In that framing, a defect in diamond is not an imperfection to be tolerated but a synthetic target to be engineered, with its atomic composition and local geometry chosen deliberately to produce the desired optical and quantum properties.

Martirez also sees a longer and somewhat recursive horizon for his field. When he thinks about the future, he finds it interesting to note that quantum chemistry may help develop tools to enrich itself in the decades to come. A practical multi-qubit system would vastly improve the level of accuracy of quantum chemistry calculations, he said. With such systems in hand, researchers could describe larger molecules, biomolecules, proteins and everything like that in a more accurate way, and that could improve all the chemical-adjacent industries, such as pharmaceutical science and materials science. Even the fundamental understanding of molecular systems could be enriched, closing a loop in which chemistry builds the quantum devices that in turn sharpen chemistry’s own predictive power.

There is also a disciplinary message in the work. Martirez encouraged other chemists to consider research in solid-state physics, arguing that at this moment in time there is a great emphasis on building the first state-of-the-art qubit computer, but that it is not just physicists who should be participating in this research. Chemists, he noted, have decades of developed intuition about chemical systems that could be easily translated into understanding defects in solids as well. The new method is, in effect, an invitation: by lowering the computational barrier to entry, it allows that chemical intuition, honed on small molecules, to be applied directly to the atomic clusters that function as qubits inside a crystal.

The research was funded by PPPL’s Laboratory Directed Research and Development program under prime contract number DE-AC02-09CH11466, with additional support from the DOE Office of Science, Office of Fusion Energy Sciences. The work reflects the laboratory’s broader mission of using plasma, the fourth state of matter, to address difficult science and technology challenges, including fusion energy, nanoscale fabrication, quantum materials and devices, and sustainability science, with Princeton University managing the laboratory for the Department of Energy’s Office of Science. For a field that has often measured progress in cryogenic kelvins and brute-force computation, the idea that a cluster of a few dozen atoms, treated as a molecule, could unlock faster and more accurate predictions of qubit behavior is a striking reframing. If diamond color centers in red, blue and ultraviolet can be designed like custom molecules, the search for room-temperature quantum devices may become less an expedition into exotic physics and more a discipline of molecular design, with chemists and physicists working side by side on the technology of the future.

Subject of Research: Computational quantum chemistry of diamond color-center defects for solid-state qubits

Article Title: Creating the technology of the future: Quantum chemistry unlocks new possibilities

Article References: Creating the technology of the future: Quantum chemistry unlocks new possibilities. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: quantum computing, qubits, diamond defects, quantum chemistry, color centers, nitrogen-vacancy, room-temperature quantum, quantum sensors, PPPL, solid-state physics, molecular simulation, ACS Central Science

The research, led by Yingying Zhao, Shengchang Zhang, Qibin Xu, Mengjin Jiang, and Pengqing Liu of the College of Polymer Science and Engineering at Sichuan University, together with Krishnan Jayaraman of the University of Auckland, tackles a problem that has stubbornly resisted experimental analysis. At the boundary where a basalt fiber meets a polymer matrix, several factors operate simultaneously: the chemical groups grafted onto the fiber surface, the nanoscale roughness of that surface, and the way resin molecules arrange themselves in response. In a real composite, all of these effects are entangled. Measure the strength of an interface and you cannot easily tell whether the credit belongs to chemical bonding, to mechanical interlocking, or to something subtler still. The Sichuan University team’s solution was to decouple these influences using atomistic simulations, allowing each factor to be varied and evaluated in isolation.

Central to the study is a workhorse of composite chemistry: the silane coupling agent known as gamma-aminopropyl triethoxysilane, or KH550. Coupling agents are molecular intermediaries, molecules designed with one end that bonds to an inorganic surface and another that is compatible with organic polymers. KH550 carries amino groups, which are rich in nitrogen and hydrogen atoms capable of forming hydrogen bonds with surrounding resin chains. Using computational simulations, the researchers showed that KH550 does something rather elegant on the basalt fiber surface: it balances highly positive and highly negative electrostatic regions, creating multiple binding sites where resin molecules can attach. In effect, the coupling agent transforms a chemically patchy surface into a welcoming landscape for polymer chains.

The quantitative results are striking. At the interface between basalt fiber and polyamide 66, a common engineering thermoplastic, KH550 modification increased the number of hydrogen bonds by 10.4 percent. More tellingly, the interfacial interaction energy, a measure of how strongly the two phases attract each other, became more negative, shifting from -384.4 to -487.9 kilocalories per mole. In the language of molecular simulation, a more negative interaction energy means a more stable, more tightly bound interface. That single modification, applied at the molecular scale, produced a measurable and substantial strengthening of the chemical handshake between fiber and resin.

But chemistry alone was only half the story, and this is where the study’s decoupling approach proves its worth. The researchers examined three different composite systems: basalt fiber with polyamide 66, basalt fiber with epoxy resin, and basalt fiber with polypropylene. Across these systems, they separated the contributions of KH550’s amino groups from the contributions of surface morphology, the atomic-scale texture of the fiber surface. The verdict was consistent: the chemical interactions provided by the amino groups contributed more to interfacial mechanical enhancement than surface morphology did. At the basalt fiber and polyamide 66 interface specifically, introducing amino groups enhanced the interfacial mechanical response by a factor of 1.49, while modifying the atomic-scale morphology enhanced it by a factor of 1.29. Both matter, but chemistry leads.

Perhaps the most intriguing finding concerns how surface texture influences the resin itself. The simulations revealed that for both the polyamide 66 and epoxy systems, atomic-level surface morphology regulates the arrangement, or conformation, of resin molecular chains, and that this effect optimizes mechanical properties more effectively than larger surface asperities that simply allow resin to infiltrate surface irregularities. In other words, the classical picture of mechanical interlocking, where rough surfaces grip the polymer like microscopic Velcro, is not the dominant mechanism at these scales. Instead, the fine texture of the surface shapes how polymer chains fold and orient, and those conformations determine how load is transferred across the interface.

KH550 also orchestrates a dramatic reorganization of the resin. According to the simulations, the coupling agent transforms resin molecular chains from conformations that lie aligned along the fiber surface into interpenetrating structures that reach into the matrix. The researchers attribute this to two combined effects: enhanced interaction between the amino groups and the resin, and steric hindrance, the physical crowding that prevents chains from lying flat against the surface. The result is a molecularly interlocked interphase, a region where fiber and polymer are no longer simply touching but genuinely entangled at the chain level.

To connect these molecular pictures to measurable mechanics, the team performed atomic-scale nanoindentation simulations, a computational analogue of the classic hardness test in which a sharp tip presses into a material to probe its resistance to deformation. The results revealed a fundamental shift in how the material fails. In unmodified interfaces, the resin deformed by sliding along the interface, a mode that ultimately leads to adhesive failure, where the bond between fiber and matrix simply gives way. After KH550 modification, deformation was redirected into the bulk of the matrix itself, a transition that shifts failure from adhesive to cohesive. Cohesive failure, in which the polymer itself must tear rather than the bond peeling apart, is generally the more desirable mode because it means the interface is no longer the weakest link. The maximum indentation force rose from approximately 2.9 to 4.0 nanonewtons, a roughly 38 percent increase attributable to the molecular redesign of the interphase.

The significance of this work extends well beyond basalt fiber. Fiber-reinforced composites of every kind, from glass and carbon to natural fibers, depend on the same delicate interfacial chemistry, and designers have traditionally improved interfaces by trial and error, testing coupling agents and surface treatments empirically. By rigorously separating the roles of chemical functionality, atomic-scale morphology, and resin conformation, the study provides something rarer: a mechanistic roadmap for rational interphase design. It tells engineers which lever to pull first. If chemistry contributes more than texture, then investment in tailored coupling agents and surface functionalization should take priority over aggressive roughening treatments, and the choice of coupling agent should be matched to the electrostatic and hydrogen-bonding character of the specific resin.

There are also practical implications for sustainability and materials selection. Basalt fiber is produced from abundant volcanic rock without the energy-intensive precursor chemistry required for carbon fiber, and it is compatible with a range of resins, from commodity polypropylene to high-performance epoxies and polyamides. The finding that the same coupling agent strategy strengthens interfaces across polyamide 66, epoxy, and polypropylene systems suggests a versatile design principle that could accelerate the adoption of basalt composites in automotive parts, wind turbine components, construction materials, and consumer products. As the authors note, the research clarifies the intrinsic interfacial strengthening mechanism and provides guidance for rational interphase design. What was once an opaque tangle of competing effects is now a decoupled, quantified, and ultimately engineerable set of levers, and that clarity may prove to be the most valuable property of all.

Subject of Research: Interfacial strengthening mechanisms in basalt fiber-reinforced polymer composites

Article Title: Decoupling of factors influencing interface characteristics in basalt fiber-reinforced polymer composites

Article References: Zhao, Y., Zhang, S., Xu, Q., Jayaraman, K., Jiang, M., & Liu, P. (2026). Decoupling of factors influencing interface characteristics in basalt fiber-reinforced polymer composites. Advanced Composites and Hybrid Materials. https://doi.org/10.1007/s42114-026-02118-6

Image Credits: AI Generated

DOI: 10.1007/s42114-026-02118-6

Keywords: basalt fiber, polymer composites, silane coupling agent, KH550, interface, hydrogen bonds, molecular simulation, nanoindentation, polyamide 66, epoxy resin, interphase design, surface morphology

The original Perspective had suggested that continuum-scale models, including surface complexation models and reactive transport models, are not designed to incorporate spatially differing reactivities of surfaces, and that introducing probability distributions to describe surface properties would transform how these models are built. Criscenti counters with a detailed historical record. As early as 1978, Davis, James and Leckie reported that the adsorption of metal ions at the oxide-water interface could not be described by assuming a homogeneous solid surface. In 1981, Benjamin and Leckie showed that the adsorption of cadmium, copper, zinc and lead onto amorphous iron oxyhydroxide could only be reproduced by assuming a large dispersion of site affinities across the surface. In other words, the energetic heterogeneity of oxide surfaces, and the statistical machinery needed to describe it, was foundational to the discipline from its earliest days.

The technical literature Criscenti marshals is extensive. Review chapters by van Riemsdijk and Koopal in 1992, Koopal in 1996, and Rudzinski and colleagues in 1996 all document early surface complexation models formulated explicitly with probability distributions describing surface heterogeneity. Koopal and van Riemsdijk went further in 1989, developing electrical double layer theory for three distinct types of spatial heterogeneity: random heterogeneity, isolated patches, and interacting patches. That work directly contradicts the Perspective’s assertion that continuum-scale models cannot handle spatially varying reactivity. Rudzinski and co-workers later fitted electrolyte ion adsorption data at the silica-sodium chloride interface using a single-site triple layer model in which heterogeneity parameters were proportional to the variance of a Gaussian-like function describing the dispersion of adsorption energies, a function originally grounded in studies of argon adsorption on silica. Charmas and colleagues extended this framework in 2002 to cadmium adsorption, fitting both surface charge titrations and electrophoretic mobility measurements across a range of pH values.

Criscenti is careful to acknowledge the parallel track that developed alongside the statistical one. Rather than describing heterogeneity with continuous distributions, some modelers chose to represent it explicitly, assigning distinct reactive groups each with its own affinity. This discrete approach underpins the MUSIC model developed by Hiemstra and colleagues in 1989, in which proton binding constants for different surface site types were estimated a priori from crystallographic considerations. The two strategies, continuous probability distributions and discrete multisite descriptions, have coexisted for decades, and Criscenti’s point is that neither represents the conceptual revolution the Perspective claimed. Both are mature tools with long track records, and choosing between them is a matter of data availability and purpose, not paradigm-shifting novelty.

The rebuttal also highlights a second, often overlooked use of probability distributions in geochemistry: quantifying uncertainty. Normal distributions have long been used to describe error bars in aqueous solution measurements, aqueous speciation equilibrium constants, and adsorption equilibrium constants. Monte Carlo methods have propagated these uncertainties through geochemical and reactive transport calculations, work Criscenti herself contributed to in the 1990s and early 2000s when assessing contaminant migration from uranium mill tailings, hazardous landfills, and potential nuclear waste repository sites. That tradition continues today. A community data-driven database and workflow for raw sorption data, developed by Zavarin and colleagues in 2022, propagates error bars from raw potentiometric titration data through to fitted acidity constants and representative average values for the community database. Even newer hybrid machine learning approaches, such as the random forest method of Chang and colleagues in 2023 for predicting radionuclide sorption, incorporate error propagation based on means and standard deviations of normal probability distributions.

Beyond the question of novelty, Criscenti identifies what she sees as misunderstandings of surface complexation model theory embedded in the Perspective. One claim held that chemical adsorption energies and Coulombic terms cannot be treated as constants because of intrinsic surface heterogeneity. Criscenti responds that these energies are indeed treated as constants for a given set of conditions, but that they properly vary with surface site type, surface coverage, solution composition and pH, and that models account for this through multiple site types and reaction-specific electrostatic terms. Another claim suggested that conventional models use only a single acidity constant and complexation constant. Criscenti points to the widely used diffuse double layer model of Dzombak and Morel, incorporated into the geochemical code PHREEQC, which requires two site types for cation sorption on hydrous ferric oxide, including a low-density, high-affinity site attributed to defects and impurities. The charge distribution MUSIC model, established enough to appear in readily available geochemical codes, allows for many more site types still.

Perhaps the most technically charged disagreement concerns ion pairing at the interface. The Perspective argued that models should account for surface-promoted ion pairing, such as strontium adsorbing as a strontium chloride ion pair subject to reduced lateral Coulomb repulsion, and claimed they currently do not. Criscenti’s rebuttal walks through the evidence that present-day models do exactly this. James and Healy showed as early as 1972, using Born solvation theory, that aqueous metal hydroxide species should adsorb preferentially over bare metal ions. Since then, numerous studies have proposed surface complexes involving the metal cation together with hydroxide or electrolyte anions, a literature comprehensively reviewed by Criscenti and Sverjensky before 2000. In extended triple layer models, the Coulombic term in the equilibrium expression changes according to the reaction, so that complexes carrying different charges and occupying different planes carry different electrostatic corrections. The curvature of metal adsorption isotherms with increasing coverage has long been interpreted as arising from sequential adsorption onto different site types or from changes in surface charge and potential as the dominant complex type shifts.

So where should the field go, if not back to distributions it never really left? Criscenti argues for more comprehensive data collection on bulk systems, since surface complexation theory is constrained by fitting calorimetric, radiometric, titration and adsorption data, and complete datasets remain scarce. Molecular-scale studies were originally recruited precisely to reduce the number of unconstrained variables, for example when several different proposed surface complexes fit the same batch adsorption data and spectroscopy or molecular modeling is needed to narrow the field. Recent work shows how powerful this integration can be: an SCM built using intrinsic acidity constants calculated from first-principles molecular dynamics successfully predicts montmorillonite surface properties as a function of pH and ionic strength, and has been extended to arsenic and uranium adsorption on clay edges. Criscenti suggests this multiscale integration provides more mechanistic information than circling back to a purely statistical representation of heterogeneity, and points to particle-continuum hybrid methods and micro-continuum models, already used for dissolving interfaces, alteration layers and mineral precipitation, as promising but still untested avenues for full multicomponent geochemical systems.

She is equally candid about the weaknesses on every side. A major limitation of large-scale molecular simulations is their still-poor description of aqueous solutions, particularly water dissociation at surfaces and the formation of metal hydroxide and metal-electrolyte complexes. Even the simplest batch adsorption experiment, one solid, one electrolyte and one trace metal, cannot yet be fully described by molecular simulation because aqueous speciation changes with pH in ways the simulations struggle to capture, making direct comparison with macroscopic experiments difficult. Meanwhile, the acknowledged soft spot of surface complexation models is the Stern layer, the near-surface region whose behavior is exquisitely sensitive to surface topography, charge distribution, ion-ion correlations and interfacial hydrogen bonding. Any new molecular-level description of that layer, she stresses, must be coupled back into geochemical equilibrium and reactive transport codes to be useful.

The takeaway from this exchange reaches beyond one contested Perspective. Criscenti’s argument is a reminder that in mature fields, proposals framed as revolutions may be rediscoveries, and that the collective memory of a discipline, preserved in its older literature, is itself a scientific resource. Her prescription is pragmatic: molecular-scale experiments and computations should be designed from the outset to feed continuum-scale constructs, with macroscopic observations guiding the questions asked at the atomic scale. Whether the original authors respond in kind remains to be seen, but the episode has already performed a service, dusting off four decades of heterogeneity modeling and putting the question of how best to represent interfacial reactivity, probabilistically or mechanistically, squarely back on the table.

Subject of Research: Probability distributions and surface complexation modeling of adsorption at oxide-water interfaces

Article Title: Should we return to probability distributions to represent interfacial reactivity?

Article References: Criscenti, L. J. (2026). Should we return to probability distributions to represent interfacial reactivity?. Nature Communications, 17(1), Article 10216. https://doi.org/10.1038/s41467-026-78022-w

Image Credits: AI Generated

DOI: 10.1038/s41467-026-78022-w

Keywords: interfacial reactivity, surface complexation models, probability distributions, surface heterogeneity, oxide-water interface, adsorption, electrical double layer, geochemical modeling, reactive transport, Stern layer, molecular simulation, Monte Carlo uncertainty

The research team, led by Xiaoqing Liu of Shanxi Datong University together with colleagues from Liaoning Technical University, constructed four molecular models of coal spanning nearly the entire coalification spectrum: lignite, long-flame coal, coking coal, and anthracite. These were not arbitrary cartoons of coal. Each model was built from real experimental data, including elemental analysis, X-ray photoelectron spectroscopy, and solid-state carbon-13 nuclear magnetic resonance spectroscopy, yielding molecular formulas such as C124H99O39N for the lignite and C150H65O11NS for the anthracite. The models captured how coalification progressively condenses aromatic ring structures, shortens aliphatic side chains, and strips away oxygen-containing functional groups, producing an increasingly ordered and rigid macromolecular framework as rank increases.

With these validated structures in hand, the researchers applied isotropic compressive stress to the model unit cells at levels ranging from 40 to 200 megapascals, chosen to reflect the in situ stresses encountered in deep coal seams and beyond. The compression produced exactly what one would expect geometrically: pore free volume, pore surface area, and porosity all fell, while the density of the models rose. But the gas behavior that followed was anything but straightforward. As stress increased, the adsorption of both methane and carbon dioxide traced a distinctive decrease–increase–decrease pattern, a nonlinear signature born from two competing effects acting simultaneously at the molecular scale.

The explanation lies in a tug-of-war between space and force. Compression shrinks the effective adsorption space available inside the coal, which should reduce how much gas the coal can hold. At the same time, however, squeezing the pore walls closer together enhances the overlap of the potential energy fields emanating from those walls, deepening the adsorption wells and strengthening the attraction between gas molecules and the coal surface. In certain stress ranges, this second effect briefly wins, producing a recovery in adsorption capacity before the loss of pore space ultimately dominates and drives adsorption down again. The net result across the full stress range was inhibitory: when stress climbed from zero to 0.20 gigapascals, saturated methane adsorption fell by as much as 43.75 percent in the lignite model, and carbon dioxide adsorption fell by similar margins.

One of the study’s most striking discoveries is that methane is more sensitive to stress than carbon dioxide. The reason is deceptively simple: methane molecules are physically larger. Their greater kinetic diameter makes them more vulnerable to being blocked in compressed pores, so when the coal matrix tightens, methane suffers disproportionately. Carbon dioxide, with its smaller molecular footprint and stronger electrostatic interactions with the polar functional groups in coal, consistently out-adsorbed methane under every condition tested. This asymmetry matters enormously for CO2-enhanced coalbed methane recovery, the technique in which injected carbon dioxide displaces methane from the coal while simultaneously locking away a greenhouse gas underground.

The energy analysis reinforced this hierarchy. Adsorbed gases showed unimodal energy distributions whose peaks shifted toward more negative values as stress increased, indicating that individual adsorption sites became more energetically favorable under compression. In the lignite model, for example, the peak adsorption energy of methane deepened from −5.45 to −5.95 kilocalories per mole across the stress range, while carbon dioxide peaked at values between −7.45 and −7.95 kilocalories per mole. Decomposing the interactions revealed that methane adsorption is almost entirely governed by van der Waals forces, accounting for more than 99 percent of the total interaction, whereas carbon dioxide’s interaction includes a substantial electrostatic component alongside van der Waals contributions of roughly 64 to 82 percent.

Diffusion told a parallel but distinct story. The self-diffusion coefficients of both gases fell steadily under stress, with methane dropping by 46 to 65 percent and carbon dioxide by 39 to 60 percent across the four coal models. Visualizing the diffusion trajectories made the mechanism vivid: under stress-free conditions, most gas molecules wandered between pores along irregular, streamline-like paths, colliding frequently with pore walls and covering large distances. Under 0.20 gigapascals of compression, however, the trajectories collapsed into tight clustered patterns, showing that molecules had become trapped oscillating within isolated pores. Stress, in effect, imprisons the gas in its nanoscale cages.

Coal rank added another layer of complexity. Adsorption capacity rose steadily with coalification, mirroring the growth in pore surface area and free volume in the higher-rank models. Diffusion, by contrast, followed a U-shaped trend across the rank spectrum, dipping in the middle-rank coals before recovering in anthracite. Carbon dioxide diffused faster than methane in every model despite its stronger binding to the coal, because its smaller kinetic diameter lets it slip through confined micropore channels more easily. The practical implication is that carbon dioxide can penetrate the coal matrix more rapidly at the molecular scale, giving it an edge in the competitive race to occupy adsorption sites and displace methane during enhanced recovery operations.

The simulations also probed the mechanical integrity of the gas-bearing coal itself. Bulk and shear moduli increased with both compressive stress and coal rank, while compressibility declined, reflecting the closure of compressible pore space and the rigid aromatic framework of mature coals. Anthracite emerged as the most mechanically robust and structurally stable of the four. Cohesive energy density, a molecular-scale gauge of how tightly the organic matter holds together, climbed with both stress and rank, and was consistently higher in carbon dioxide-bearing models than in methane-bearing ones, suggesting that carbon dioxide injection and loading both strengthen intermolecular cohesion within the simulated coal systems.

For a field racing to balance energy production with climate goals, these molecular-level insights arrive at an opportune moment. Coalbed methane is a lower-emission fossil fuel than coal itself, and CO2-enhanced recovery promises a dual dividend: boosting methane output while permanently sequestering carbon dioxide in the same seams. By showing precisely how stress reshapes the adsorption and diffusion landscape—and why methane and carbon dioxide respond so differently to the squeeze—the study gives reservoir engineers a firmer physical foundation for designing extraction and sequestration strategies in the deep, high-stress coal seams where the industry’s future increasingly lies.

Subject of Research: Molecular simulation of methane and carbon dioxide adsorption and diffusion in coals of different ranks under compressive stress loading

Article Title: Molecular simulation of adsorption and diffusion behavior of CH4/CO2 in coals of different ranks under stress loading

Article References: Liu, X., Liang, F., Xu, Q., Zhao, D., Bai, W., & Yang, X. (2026). Molecular simulation of adsorption and diffusion behavior of CH4/CO2 in coals of different ranks under stress loading. Journal of Saudi Chemical Society, 30(5), Article 73. https://doi.org/10.1007/s44442-026-00128-z

Image Credits: AI Generated

DOI: 10.1007/s44442-026-00128-z

Keywords: coalbed methane, molecular simulation, adsorption, diffusion, compressive stress, coal rank, carbon dioxide sequestration, CO2-ECBM, GCMC, molecular dynamics, coal macromolecular model, pore structure