SEARCH · Search NASA

Results for “Order fractionalization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Entropy of Strongly Correlated Electrons in a Partially Filled Landau Level

We use high-resolution chemical potential measurements to extract the entropy of monolayer and bilayer graphene in the quantum Hall regime via the Maxwell relation [(𝑑⁢𝜇)/(𝑑⁢𝑇)]| 𝑁 =−[(𝑑⁢𝑆)/(𝑑⁢𝑁)]| 𝑇 . Measuring the entropy from 𝑇 = 300 K down to 𝑇 = 200 mK, we identify the sequential emergence of quantum Hall ferromagnetism, fractional quantum Hall states, and various charge orders by comparing the measured entropy in different temperature regimes with theoretical models. At the lowest temperature of 𝑇 ≈ 200 mK we perform a detailed study of the entropy near even-denominator fractional quantum Hall states in bilayer graphene, and comment on the possible topological origin of the observed excess entropy.

Assouline, Alexandre [University of California, Sa↗

DOE OSTI2025

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗

DOE OSTI2021

Dynamic fingerprint of fractionalized excitations in single-crystalline Cu3Zn(OH)6FBr

Abstract Beyond the absence of long-range magnetic orders, the most prominent feature of the elusive quantum spin liquid (QSL) state is the existence of fractionalized spin excitations, i.e., spinons. When the system orders, the spin-wave excitation appears as the bound state of the spinon-antispinon pair. Although scarcely reported, a direct comparison between similar compounds illustrates the evolution from spinon to magnon. Here, we perform the Raman scattering on single crystals of two quantum kagome antiferromagnets, of which one is the kagome QSL candidate Cu 3 Zn(OH) 6 FBr, and another is an antiferromagnetically ordered compound EuCu 3 (OH) 6 Cl 3 . In Cu 3 Zn(OH) 6 FBr, we identify a unique one spinon-antispinon pair component in the E 2g magnetic Raman continuum, providing strong evidence for deconfined spinon excitations. In contrast, a sharp magnon peak emerges from the one-pair spinon continuum in the E g magnetic Raman response once EuCu 3 (OH) 6 Cl 3 undergoes the antiferromagnetic order transition. From the comparative Raman studies, we can regard the magnon mode as the spinon-antispinon bound state, and the spinon confinement drives the magnetic ordering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Unavoidable multilevel biaxial symmetry breaking in chiral hybrid liquid crystals

Chiral nematic or cholesteric liquid crystals (LCs) are fluid mesophases with long-ranged orientational order featuring a quasilayered periodicity imparted by a helical director configuration but lacking long-range positional order. Doping molecular cholesteric LCs with strongly anisotropic uniaxial colloidal particles adds another level of complexity because of the interplay between weak surface-anchoring boundary conditions and bulk-based elastic distortions near the particle-LC interface. Using cylindrical colloidal disks and rods with different geometric shapes and surface conditions, we demonstrate that these colloidal inclusions generically exhibit biaxial orientational probability distributions which may impart anomalously strong local biaxiality onto the hybrid cholesteric LC structure. Unlike nonchiral hybrid molecular-colloidal LCs, where biaxial order emerges only at critical colloid volume fractions exceeding some uniaxial-biaxial transition value, the orientational probability of the colloidal inclusions immersed in chiral nematic hosts is unambiguously biaxial even at infinite dilution. We demonstrate that the colloids induce local biaxial perturbations within the molecular orientational order of the LC host medium which strongly enhances the weak but native biaxial order of chiral nematic LC induced by the chiral symmetry breaking of the director field. With the help of analytical modeling and computer simulations based on the Landau–de Gennes free energy of the host LC around the colloids, we rationalize the observed multilevel biaxial order and conclude that it is not only unavoidable but also strongly enhanced compared to both achiral hybrid LCs and purely molecular cholesteric LCs.

biaxial liquid crystals↗

DOE OSTI2022

Phase evolution of U-Zr system in a thermal cycling neutron diffraction experiment: as-cast U-35Zr and U-50Zr

Here, phase evolution of as-cast U-35 wt% Zr and U-50 wt% Zr alloys during thermal cycling (303–1073 K) was investigated using in-situ neutron diffraction. Analysis was performed using Rietveld crystal and microstructure refinements from time-of-flight neutron diffraction data, with a focus on evolution of lattice parameter, atom ordering, unit cell volume, and weight fractions during the thermal cycling. Disordered δ-UZr 2 and residual γ-(U,Zr) phase are retained in the field of δ of the U-Zr equilibrium phase diagram for the U-50Zr sample, and in the field between δ and α-U for the U-35Zr sample. The α phase that is present in the reported phase diagrams of U-Zr was not observed in the diffraction patterns. The evolutions of lattice parameter and unit cell volume of δ and γ are affected by both thermal expansion and chemistry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

A seamless multiscale operator neural network for inferring bubble dynamics

Modelling multiscale systems from nanoscale to macroscale requires the use of atomistic and continuum methods and, correspondingly, different computer codes. Here, we develop a seamless method based on DeepONet, which is a composite deep neural network (a branch and a trunk network) for regressing operators. In particular, we consider bubble growth dynamics, and we model tiny bubbles of initial size from 100 nm to 10 $\mathrm {\mu }\textrm {m}$ , modelled by the Rayleigh–Plesset equation in the continuum regime above 1 $\mathrm {\mu }\textrm {m}$ and the dissipative particle dynamics method for bubbles below 1 $\mathrm {\mu }\textrm {m}$ in the atomistic regime. After an offline training based on data from both regimes, DeepONet can make accurate predictions of bubble growth on-the-fly (within a fraction of a second) across four orders of magnitude difference in spatial scales and two orders of magnitude in temporal scales. The framework of DeepONet is general and can be used for unifying physical models of different scales in diverse multiscale applications.

Mechanics↗

DOE OSTI2025

Simultaneous Absorption and Desorption Isotherms of Various Hydrogen and Deuterium Mixtures between 20 and 120 °C Are Used to Determine the Activity Coefficients for Palladium-Hydride Solutions

A test bed was constructed to measure the absorption and desorption isotherms for palladium hydride using H 2 , D 2 , and various H 2 /D 2 mixtures for temperatures in the range 20 °C ≤ T ≤ 120 °C. The pressure–composition–temperature isotherms were measured. The pressures obtained with mixtures between each pure isotope were monotonic, yet nonlinear. This nonlinear dependence of total pressure with feed gas protium concentration reveals the mixed isotope hydride system behaves nonideally. A thermodynamic model was adapted from the literature, which accounts for the nonideal nature of the mixed-isotope system. The measured data were used as constraints in this model in order to calculate the protium mole fractions in the hydride phase and the protium activity coefficients for palladium hydride at various temperatures and protium concentrations in the system. Knowing the protium mole fraction and activity coefficients allows for a priori calculation of the isotopologue distribution in the gas phase and the isotope distribution in the hydride phase, given a palladium temperature and equilibrium pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Tunable spin and valley excitations of correlated insulators in Γ-valley moiré bands

Moiré superlattices formed from transition metal dichalcogenides (TMDs) have been shown to support a variety of quantum electronic phases that are highly tunable using applied electromagnetic fields. While the valley character of the low-energy states dramatically affects optoelectronic properties in the constituent TMDs, this degree of freedom has yet to be fully explored in moiré systems. Here, we establish twisted double bilayer WSe 2 as an experimental platform to study electronic correlations within Γ-valley moiré bands. Through a combination of local and global electronic compressibility measurements, we identify charge-ordered phases at multiple integer and fractional moiré band fillings ν. By measuring the magnetic field dependence of their energy gaps and the chemical potential upon doping, we reveal spin-polarized ground states with novel spin polaron quasiparticle excitations. In addition, an applied displacement field allows us to realize a new mechanism of metal-insulator transition at ν=–1 driven by tuning between Γ- and K-valley moiré bands. Together, our results demonstrate control over both the spin and valley character of the correlated ground and excited states in this system.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Laser cooling of antihydrogen atoms

The photon—the quantum excitation of the electromagnetic field—is massless but carries momentum. A photon can therefore exert a force on an object upon collision 1 . Slowing the translational motion of atoms and ions by application of such a force 2,3 , known as laser cooling, was first demonstrated 40 years ago 4,5 . It revolutionized atomic physics over the following decades 6–8 , and it is now a workhorse in many fields, including studies on quantum degenerate gases, quantum information, atomic clocks and tests of fundamental physics. However, this technique has not yet been applied to antimatter. Here we demonstrate laser cooling of antihydrogen 9 , the antimatter atom consisting of an antiproton and a positron. By exciting the 1S–2P transition in antihydrogen with pulsed, narrow-linewidth, Lyman-α laser radiation 10,11 , we Doppler-cool a sample of magnetically trapped antihydrogen. Although we apply laser cooling in only one dimension, the trap couples the longitudinal and transverse motions of the anti-atoms, leading to cooling in all three dimensions. We observe a reduction in the median transverse energy by more than an order of magnitude—with a substantial fraction of the anti-atoms attaining submicroelectronvolt transverse kinetic energies. We also report the observation of the laser-driven 1S–2S transition in samples of laser-cooled antihydrogen atoms. The observed spectral line is approximately four times narrower than that obtained without laser cooling. The demonstration of laser cooling and its immediate application has far-reaching implications for antimatter studies. A more localized, denser and colder sample of antihydrogen will drastically improve spectroscopic 11–13 and gravitational 14 studies of antihydrogen in ongoing experiments. Furthermore, the demonstrated ability to manipulate the motion of antimatter atoms by laser light will potentially provide ground-breaking opportunities for future experiments, such as anti-atomic fountains, anti-atom interferometry and the creation of antimatter molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Estimating the local dark matter density in a non-axisymmetric wobbling disc

The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Fate of multiparticle resonances: From Q-balls to 3 He droplets

We study a system of N nonrelativistic particles which form a near-threshold resonance. Assuming no subset of these particles can form a bound state, the resonance can only decay through an “explosion” into N particles. We find that the decay width of the resonance scales as E Δ –5/2 in the limit when the energy E of the resonance goes to zero, where Δ is the ground-state energy of a system of N particles in a spherical harmonic trap with unit frequency. Here, the formula remains valid when some pairs of final particles have zero-energy s-wave resonance, but the Efimov effect is not present. In the limit of large N, we show that the final particles follow a Maxwell-Boltzmann distribution if they are bosons and a semicirclelike law if they are fermions. We expect our general result to be applicable to various systems that exist in nature. In particular, we argue that metastable 3 He droplets exist with the lifetime varying over many orders of magnitude ranging from a fraction of a nanosecond to values greatly exceeding the age of the Universe.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Observation of the Rare Decay of the 𝜂 Meson to Four Muons

A search for the rare η → μ + μ − μ + μ − double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb −1 . A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the η → μ + μ − decay as normalization, the branching fraction B(η → μ + μ − μ + μ − ) = [5.0 ± 0.8 (stat) ± 0.7 (syst) ± 0.7 (B 2μ )] × 10 −9 is measured, where the last term is the uncertainty in the normalization channel branching fraction. This work achieves an improved precision of over five orders of magnitude compared to previous results, leading to the first measurement of this branching fraction, which is found to agree with theoretical predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Above-ground net primary productivity in regenerating seasonally dry tropical forest: Contributions of rainfall, forest age and soil

Identifying factors controlling forest productivity is critical to understanding forest-climate change feedbacks, modelling vegetation dynamics and carbon finance schemes. However, little research has focused on productivity in regenerating tropical forests which are expanding in their fraction of global area have an order of magnitude larger carbon uptake rates relative to older forest. Here We examined above-ground net primary productivity (ANPP) and its components (wood production and litterfall) over 10 years in forest plots that vary in successional age, soil characteristics and species composition using band dendrometers and litterfall traps in regenerating seasonally dry tropical forests in northwestern Costa Rica. We show that the components of ANPP are differentially driven by age and annual rainfall and that local soil variation is important. Total ANPP was explained by a combination of age, annual rainfall and soil variation. Wood production comprised 35% of ANPP on average across sites and years, and was explained by annual rainfall but not forest age. Conversely, litterfall increased with forest age and soil fertility yet was not affected by annual rainfall. In this region, edaphic variability is highly correlated with plant community composition. Thus, variation in ecosystem processes explained by soil may also be partially explained by species composition. These results suggest that future changes in annual rainfall can alter the secondary forest carbon sink, but this effect will be buffered by the litterfall flux which varies little among years. In determining the long-term strength of the secondary forest carbon sink, both rainfall and forest age will be critical variables to track. We also conclude that detailed understanding of local site variation in soils and plant community may be required to accurately predict the impact of changing rainfall on forest carbon uptake. Synthesis. We show that in seasonally dry tropical forest, annual rainfall has a positive relationship with the growth of above-ground woody tissues of trees and that droughts lead to significant reductions in above-ground productivity. These results provide evidence for climate change—carbon cycle feedbacks in the seasonal tropics and highlight the value of longitudinal data on forest regeneration.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Improved Value of the Gasoline and Fuel Oil Co-Product Fractions Generated by the PNNL/LanzaTech Alcohol-to-Jet Process (Final Report)

The BETO-funded PNNL collaboration with LanzaTech, Inc. utilizes patented technology to address the need for low aromatic jet fuel blendstocks. The PNNL/LanzaTech alcohol-to-jet (ATJ) process converts ethanol from LanzaTech’s syngas fermentation process to jet-range isoparaffins. The aim of this Cooperative Research Development Agreement (CRADA) project was to improve the commercial viability of the LanzaTech/ PNNL ATJ process by developing two new co-product options. Two Focus Areas (FA) were explored: FA1) Increasing the RON of the lighter-than-jet gasoline fraction above 98; and FA2) Creating a synthetic lubricant base oil from the heavier-than-jet fraction. For FA 1, multiple processing approaches were evaluated, and RONs ranging from 96-100 were obtained, each comprising different processing costs. The most cost effective approach identified was a single-step process and yielded a liquid product of 97. Thus, additional development is required to obtain a RON > 98. We will continue work on this focus area with separate Direct Funding Opportunity (DFO) funding within the DOE-BETO Co-Optima Consortium (expected to begin Q2-FY21). Experimental work will focus on i) increasing single pass conversion to > 50%, and ii) obtaining a > 98 RON product. Technoeconomic analysis performed in this project suggests cost competitiveness provided these technical targets can be met. For FA2, multiple approaches were evaluated for increasing the viscosity index (VI) of the heavier-than-jet fraction to at least 120 in order to meet ASTM specifications for a Group III Base Oil (ASTM D6074). All other ASTM specifications are already met. A viscosity index of 115 was obtained, and with minimal undesirable cracking products. We believe a VI of 115 is about the highest possible from a catalytic approach, per review of the patent literature. In order to further increase the VI to > 120 we believe additives and/or liquid-liquid extraction is required. TEA performed earlier in the project suggests cost competitiveness with market prices assuming cracking byproducts are kept < 10 wt.%. If TEA updated with this additional processing still projects economic feasibility next steps could be taken in a subsequent effort.

02 PETROLEUM↗

DOE OSTI2024

Predictions for a Low-mass Cutoff for the Primordial Black Hole Mass Spectrum

In this note, we outline how a modest violation in the conservation of mass during the merger of two PBHs affects the PBH mass spectrum that we previously obtained using a Boltzmann equation model for the evolution of the mass spectrum with no mass loss. We find that if the initial cosmological redshift is of the order of 10 12 , then the fraction of primordial holes with masses greater than 10 3 solar masses appears to be close to what is required to provide the seeds for galaxies. In addition, we note that as a result of rapid collisions and strong coupling to electromagnetic radiation for temperatures > GeV, there will be an effective low-mass cutoff in the mass spectrum for PBH masses less than a certain PBH mass less than 0.1M ⊙ . We also point out that this cutoff in the mass spectrum below ~ 0.1M ⊙ can be confirmed by combining future microlensing observations from the Roman Space Telescope and the Vera C. Rubin Observatory with astrometric observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Assessing the Behavior of Moderately Volatile Elements on the Moon in Order to Constrain Processes of Magmatic Evolution and Planetary Accretion

Heterogeneous isotopic compositions of Zn and Ga in lunar basalts and highlands rocks indicate that secondary processing affected the stable isotopic composition of major lunar reservoirs. It follows that these fractionating processes must be identified and characterized in order to make robust estimates for the isotopic composition of the mare basalt source regions or the bulk silicate Moon (BSM). Contrary to previous assertions (e.g., Kato & Moynier, 2017), MVEs such as Zn and Ga do not behave similarly on the Moon, in line with differences in their chemical behavior and volatility. There is clear evidence that impact reworking fractionates zinc isotopes but, in addition, zinc isotopes are also highly sensitive to the formation of Ti-bearing oxides such as ilmenite. In contrast, gallium is far less susceptible to volatile reworking, but its isotopic composition is fractionated during incorporation into plagioclase (Wimpenny et al., 2022; Render et al., 2023). Based on these findings we conclude that current isotopic compositions of Zn and Ga in lunar rocks, such as those measured for the mare basalt suite, are unlikely to be representative of the bulk Moon when it accreted. As such, estimates for the MVE isotopic composition of the BSM derived from mare basalts may not be accurate, which may have implications for the incorporation of these data into dynamical models of the Moon’s formation. In the following, we discuss the main findings from this project in more detail, focusing on the behaviour of zinc and gallium isotopic systematics in lunar rocks and the processes that control their isotopic fractionation on the Moon.

58 GEOSCIENCES↗