Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spin lattice models”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Manipulating Growth and Propagation of Correlations in Dipolar Multilayers: From Pair Production to Bosonic Kitaev Models

Here, we study the nonequilibrium dynamics of dipoles confined in multiple stacked two-dimensional layers realizing a long-range interacting quantum spin 1/2 XXX model. We demonstrate that strong in-plane interactions can protect a manifold of collective layer dynamics. This then allows us to map the many-body spin dynamics to bosonic models. In a bilayer configuration we show how to engineer the paradigmatic two-mode squeezing Hamiltonian known from quantum optics, resulting in exponential production of entangled pairs and generation of metrologically useful entanglement from initially prepared product states. In multilayer configurations we engineer a bosonic variant of the Kitaev model displaying chiral propagation along the layer direction. Our study illustrates how the control over interactions, lattice geometry, and state preparation in interacting dipolar systems uniquely afforded by AMO platforms such as Rydberg and magnetic atoms, polar molecules, or trapped ions allows for the control over the temporal and spatial propagation of correlations for applications in quantum sensing and quantum simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗
DOE OSTI2024

Persistence of Spin Coherence in a Crystalline Environment

Here, we analyze quantum interference in the triplet-exciton pair generated by singlet exciton fission in a molecular crystal and introduce transport-induced dephasing (TID) as a key effect that can suppress the expected fluorescence quantum beats when the triplet-exciton wave function can localize on inequivalent sites. TID depends on the triplet-exciton hopping rate between inequivalent sites and on the energy shifts among the stationary states of the entangled triplet pair in different spatial configurations. The theoretical model is confirmed by experiments in rubrene single crystals, where triplet pairs remain entangled for more than 50 ns but quantum beats are suppressed by TID within a few nanoseconds when the magnetic field is misaligned by just a few degrees from specific symmetric directions. Our experiments deliver the zero-field parameters for the rubrene molecule in its orthorhombic lattice and information on triplet-exciton transport, in particular, the triplet-exciton hopping rate between inequivalent sites, which we evaluate to be of the order of 150 ps in rubrene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

fcc lattice, checkerboards, fractons, and quantum field theory

In this work, we consider XY-spin degrees of freedom on an FCC lattice, such that the system respects some subsystem global symmetry. We then gauge this global symmetry and study the corresponding $U(1)$ gauge theory on the FCC lattice. Surprisingly, this $U(1)$ gauge theory is dual to the original spin system. We also analyze a similar $\mathbb{Z}_N$ gauge theory on that lattice. All these systems are fractonic. The $U(1)$ theories are gapless and the $\mathbb{Z}_N$ theories are gapped. We analyze the continuum limits of all these systems and present free continuum Lagrangians for their low-energy physics. Our $\mathbb{Z}_2$ FCC gauge theory is the continuum limit of the well known checkerboard model of fractons. Our continuum analysis leads to a straightforward proof of the known fact that this theory is dual to two copies of the $\mathbb{Z}_2$ X-cube model. We find new models and new relations between known models. The $\mathbb{Z}_N$ FCC gauge theory can be realized by coupling three copies of an anisotropic model of lineons and planons to a certain exotic $\mathbb{Z}_2$ gauge theory. Also, although for $N=2$ this model is dual to two copies of the $\mathbb{Z}_2$ X-cube model, a similar statement is not true for higher $N$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantum spin ice response to a magnetic field in the dipole-octupole pyrochlore Ce 2 Zr 2 O 7

The pyrochlore magnet Ce 2 Zr 2 O 7 has attracted much attention as a quantum spin ice candidate whose novelty derives in part from the dipolar-octupolar nature of the Ce 3+ pseudospin-1/2 degrees of freedom it possesses. We report heat capacity measurements on single crystal samples of Ce 2 Zr 2 O 7 down to T~0.1K in a magnetic field along the [$1,\bar{1},0$] direction. These measurements show that the broad hump in the zero-field heat capacity moves higher in temperature with increasing field strength and is split into two separate humps by the [$1,\bar{1},0$] magnetic field at ~2T. These separate features are due to the decomposition of the pyrochlore lattice into effectively decoupled chains for fields in this direction: One set of chains (α chains) is polarized by the field while the other (β chains) remains free. This situation is similar to that observed in the classical spin ices Ho 2 Ti 2 O 7 and Dy 2 Ti 2 O 7 , but with the twist that here the strong transverse exchange interactions produce substantial quantum effects. Our theoretical modeling suggests that the β chains are close to a critical state, with nearly-gapless excitations. We also report elastic and inelastic neutron scattering measurements on single crystal Ce 2 Zr 2 O 7 in [$1,\bar{1},0$] and [0,0,1] magnetic fields at temperatures down to T=0.03K. The elastic scattering behaves consistently with the formation of independent chains for a [$1,\bar{1},0$] field, while the [0,0,1] field produces a single field-induced elastic magnetic Bragg peak at (0,2,0) and equivalent wavevectors, indicating a polarized spin ice state for fields above ~3T. For both [$1,\bar{1},0$] and [0,0,1] magnetic fields, our inelastic neutron scattering results show an approximately dispersionless continuum of scattering that increases in both energy and intensity with increasing field strength. By modeling the complete set of experimental data using numerical linked cluster and semiclassical molecular dynamics calculations, we demonstrate the dominantly multipolar nature of the exchange interactions in Ce 2 Zr 2 O 7 and the smallness of the parameter θ, which controls the mixing between dipolar and octupolar degrees of freedom. Finally, these results support previous estimates of the microscopic exchange parameters and place strong constraints on the theoretical description of this prominent spin ice candidate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Transport in magnetically doped one-dimensional wires: can the helical protection emerge without the global helicity?

We study the phase diagram and transport properties of arbitrarily doped quantum wires functionalized by magnetic adatoms. The appropriate theoretical model for these systems is a dense one-dimensional Kondo lattice (KL) which consists of itinerant electrons interacting with localized quantum magnetic moments. We discover the novel phase of the locally helical metal where transport is protected from a destructive influence of material imperfections. Paradoxically, such a protection emerges without a need of the global helicity, which is inherent in all previously studied helical systems and requires breaking the spin-rotation symmetry. We explain the physics of this protection of the new type, find conditions, under which it emerges, and discuss possible experimental tests. Our results pave the way to the straightforward realization of the protected ballistic transport in quantum wires made of various materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

How correlations change the magnetic structure factor of the kagome Hubbard model

The kagome Hubbard model (KHM) is a paradigmatic example of a frustrated two-dimensional model. While its strongly correlated regime, described by a Heisenberg model, is of topical interest due to its enigmatic prospective spin-liquid ground state, the weakly and moderately correlated regimes remain largely unexplored. Motivated by the rapidly growing number of metallic kagome materials (e.g., Mn 3 Sn, Fe 3 Sn 2 , FeSn, Co 3 Sn 2 S 2 , Gd 3 Ru 4 Al 12 , and AV 3 Sb 5 with A = K, Rb, Cs), here we study the respective regimes of the KHM by means of three complementary numerical methods: the dynamical mean-field theory, the dynamical vertex approximation, and determinant quantum Monte Carlo. In contrast to the archetypal square lattice, we find no tendencies toward magnetic ordering, as magnetic correlations remain short-range. Nevertheless, the magnetic correlations undergo a remarkable crossover as the system approaches the metal-to-insulator transition. The Mott transition itself does not affect the magnetic correlations. Our equal-time and dynamical structure factors can be used as a reference for inelastic neutron scattering experiments on the growing family of metallic kagome materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Contrasting Trivalent Lanthanide and Actinide Complexation by Polyoxometalates via Solution-State NMR

Deciphering the solution chemistry and speciation of actinides is inherently difficult due to radioactivity, rarity, and cost constraints, especially for transplutonium elements. In this context, the development of new chelating platforms for actinides and associated spectroscopic techniques is particularly important. In this study, we investigate a relatively overlooked class of chelators for actinide binding, namely, polyoxometalates (POMs). We provide the first NMR measurements on americium–POM and curium–POM complexes, using one-dimensional (1D) 31 P NMR, variable-temperature NMR, and spin-lattice relaxation time (T 1 ) experiments. The proposed POM–NMR approach allows for the study of trivalent f-elements even when only microgram amounts are available and in phosphate-containing solutions where f-elements are typically insoluble. The solution-state speciation of trivalent americium, curium, plus multiple lanthanide ions (La 3+ , Nd 3+ , Sm 3+ , Eu 3+ , Yb 3+ , and Lu 3+ ), in the presence of the model POM ligand PW 11 O 39 7– was elucidated and revealed the concurrent formation of two stable complexes, [M III (PW 11 O 39 )(H 2 O) x ] 4– and [M III (PW 11 O 39 ) 2 ] 11– . Interconversion reaction constants, reaction enthalpies, and reaction entropies were derived from the NMR data. The NMR results also provide experimental evidence of the weakly paramagnetic nature of the Am 3+ and Cm 3+ ions in solution. Furthermore, the study reveals a previously unnoticed periodicity break along the f-element series with the reversal of T 1 relaxation times of the 1:1 and 1:2 complexes and the preferential formation of the long T 1 species for the early lanthanides versus the short T 1 species for the late lanthanides, americium, and curium. Furthermore, given the broad variety of POM ligands that exist, with many of them containing NMR-active nuclei, the combined POM–NMR approach reported here opens a new avenue to investigate difficult-to-study elements such as heavy actinides and other radionuclides.

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

Nanoparticle self-assemblies with modern complexity

Thanks to decades of tireless efforts, nanoparticle assemblies have reached at an extremely high level of controllability, sophistication, and complexity, with new insights provided by integration with graph theory, cutting-edge characterization, and machine learning (ML)-based computation and modeling, as well as with ever-diversifying applications in energy, catalysis, biomedicine, optics, electronics, magnetics, organic biosynthesis, and quantum technology. Nanoparticle assemblies can be crystalline, known as superlattices or supracrystals. Their assembly entails a transition from disorder—dispersed nanoparticles—to order, which can be achieved through classical nucleation pathways or nonclassical pathways via prenucleation precursors or particle aggregation. Further, the periodic lattices allow facile manipulations of electrons, phonons, photons, and even spins, leading to advanced device components and metamaterials. Meanwhile, aperiodic assemblies out of nanoparticles, such as gels, networks, and amorphous solids, also start to attract attentions. Despite the loss of periodicity, symmetry-lowering or symmetry-breaking three-dimensional (3D) structures emerge with unique properties, such as chiroptical activity, topological mechanical strength, and quantum entanglement. Real-space imaging such as electron microscopy and X-ray based tomography methods are utilized to characterize these complex structures, while mathematical tools such as graph theories are in need to describe such complex structures. This issue aims to provide a timely review of the efforts in this greatly broadened materials design space including experiment, simulation, theory, and applications. Nine top experts (and their teams) from four countries deliver six review papers, summarizing fundamental mechanistic understandings of nanoparticle assemblies, highlighted with the developments of state-of-the-art in situ characterization tools and ML-assisted reverse engineering, and newly emergent applications of nanoarchitectures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Low-temperature grapho-epitaxial La-substituted BiFeO 3 on metallic perovskite

Bismuth ferrite has garnered considerable attention as a promising candidate for magnetoelectric spin-orbit coupled logic-in-memory. As model systems, epitaxial BiFeO 3 thin films have typically been deposited at relatively high temperatures (650–800 °C), higher than allowed for direct integration with silicon-CMOS platforms. Here, we circumvent this problem by growing lanthanum-substituted BiFeO 3 at 450 °C (which is reasonably compatible with silicon-CMOS integration) on epitaxial BaPb 0.75 Bi 0.25 O 3 electrodes. Notwithstanding the large lattice mismatch between the La-BiFeO 3 , BaPb 0.75 Bi 0.25 O 3 , and SrTiO 3 (001) substrates, all the layers in the heterostructures are well ordered with a [001] texture. Polarization mapping using atomic resolution STEM imaging and vector mapping established the short-range polarization ordering in the low temperature grown La-BiFeO 3 . Current-voltage, pulsed-switching, fatigue, and retention measurements follow the characteristic behavior of high-temperature grown La-BiFeO 3 , where SrRuO 3 typically serves as the metallic electrode. These results provide a possible route for realizing epitaxial multiferroics on complex-oxide buffer layers at low temperatures and opens the door for potential silicon-CMOS integration.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗
DOE OSTI2021

Ligand-Controlled Energetics and Charge Transfer in Pure and Doped Nanocrystals

The research in this award period initially focused on the spectroscopy and dynamics of spherical CdSe nanocrystals and CdSe nanoplatelets. In the later part of the grant period we turned our attention to InP-based nanocrystals. It has long been believed that the transient absorption signal from approximately spherical CdSe nanocrystals is dominated by conduction band state filling. However, this long-held belief has recently been challenged, based on femtosecond absorption measurements. However, the recent studies challenging the conventional wisdom do not account of finite rates of spin-lattice relaxation. Our work showed that this is a crucial error and that the recent results are misinterpreted, that is, the previous prevailing wisdom is correct. Another study focused on CdSe nanocrystal photochemistry used transient absorption (TA) spectroscopy to determine the spatial extents of CdSe nanoplatelet (NPL) excitons. Our work shows that the spatial extents of the excitons in the NPLs are far less than the physical dimension of the NPL. Using a model developed to understand the transient absorption spectroscopy, we obtain an average excitonic area of 21.2 ± 2.5 nm 2 , independent of the nanoplatelet size. Our work on InP/ZnSe and InP/ZnS core/shell nanocrystals shows that when there is a small lattice mismatch (InP-ZnSe, 3.5%) a coherent core-shell interface is obtained. In contrast, the InP-ZnS lattice mismatch is much larger, 8.3%. In this case, the experimental results showed best agreement with calculations in which lattice strain is ignored, indicating that the interfaces in InP/ZnS nanocrystals are largely incoherent.

14 SOLAR ENERGY↗
DOE OSTI2022

Excited quantum anomalous and spin Hall effect: dissociation of flat-bands-enabled excitonic insulator state

Quantum anomalous Hall effect (QAHE) and quantum spin Hall effect (QSHE) are two interesting physical manifestations of 2D materials that have an intrinsic nontrivial band topology. In principle, they are ground-state equilibrium properties characterized by Fermi level lying in a topological gap, below which all the occupied bands are summed to a non-zero topological invariant. In this work, we propose theoretical concepts and models of ‘excited’ QAHE (EQAHE) and EQSHE generated by dissociation of an excitonic insulator (EI) state with complete population inversion (CPI), a unique many-body ground state enabled by two yin-yang flat bands (FBs) of opposite chirality hosted in a diatomic Kagome lattice. The two FBs have a trivial gap in between, i.e. the system is a trivial insulator in the single-particle ground-state, but nontrivial gaps above and below, so that upon photoexcitation the quasi-Fermi levels of both electrons and holes will lie in a nontrivial gap achieved by the CPI-EI state, as demonstrated by exact diagonalization calculations. Then dissociation of singlet and triplet EI state will lead to EQAHE and EQSHE, respectively. Realizations of yin-yang FBs in real materials are also discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Evolution of magnetic field induced ordering in the layered quantum Heisenberg triangular-lattice antiferromagnet Ba 3 CoSb 2 O 9

Quantum fluctuations in the effective spin-12 layered triangular-lattice quantum Heisenberg antiferromagnet Ba 3 CoSb 2 O 9 lift the classical degeneracy of the antiferromagnetic ground state in magnetic field, producing a series of novel spin structures for magnetic fields applied within the crystallographic ab plane, including a celebrated collinear “up-up-down” spin ordering with magnetization equal to 1/3 of the saturation magnetization over an extended field range. Theoretically unresolved, however, are the effects of interlayer antiferromagnetic coupling and transverse magnetic fields on the ground states of this system. Additional magnetic field induced phase transitions are theoretically expected and in some cases have been experimentally observed, but details regarding their number, location, and physical character appear inconsistent with the predictions of existing models. Conversely, an absence of experimental measurements as a function of magnetic-field orientation has left other key predictions of these models untested. To address these issues, in this study we have used specific heat, neutron diffraction, thermal conductivity, and magnetic torque measurements to map out the phase diagram as a function of magnetic field intensity and orientation relative to the crystallographic ab plane. For H||ab, we have discovered an additional magnetic field induced phase transition at low temperature and an unexpected tetracritical point in the high-field phase diagram, which coupled with the apparent second-order nature of the phase transitions eliminates several theoretically proposed spin structures for the high-field phases. Our calorimetric measurements as a function of magnetic field orientation are in general agreement with theory for field-orientation angles close to plane parallel (H||a) but diverge at angles near plane perpendicular; a predicted convergence of two phase boundaries at finite angle and a corresponding change in the order of the field induced phase transition are not observed experimentally. Our results emphasize the role of interlayer coupling in selecting and stabilizing field induced phases, provide guidance on the nature of the magnetic order in each phase, and reveal the need for new physics to account for the nature of magnetic ordering in this archetypal two-dimensional spin-12 triangular-lattice quantum Heisenberg antiferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Disentangling the physics of the attractive Hubbard model as a fully interacting model of fermions via the accessible and symmetry-resolved entanglement entropies

The complicated ways in which electrons interact in many-body systems such as molecules and materials have long been viewed through the lens of local electron correlation and associated correlation functions. However, quantum information science has demonstrated that more global diagnostics of quantum states like the entanglement entropy can provide a complementary and clarifying lens on electronic behavior. One particularly useful measure that can be used to distinguish between quantum and classical sources of entanglement is the accessible entanglement, the entanglement available as a quantum resource for systems subject to conservation laws, such as fixed particle number, due to superselection rules. In this work, we introduce an algorithm and demonstrate how to compute accessible and symmetry-resolved entanglements for interacting fermion systems. This is accomplished by combining an incremental version of the swap algorithm with a recursive auxiliary field quantum Monte Carlo algorithm recently developed by the authors. We apply these tools to study the pairing and charge density waves exhibited in the paradigmatic attractive Hubbard model via entanglement. We find that the particle and spin symmetry-resolved entanglements and their related full probability distribution functions show very clear—and unique—signatures of the underlying electronic behavior even when those features are less pronounced in conventional correlation functions. Altogether, this work provides a systematic means of characterizing the entanglement within quantum systems that can grant a deeper understanding of the complicated electronic behavior that underlies quantum phase transitions and crossovers in many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Magnetic excitations from the hexagonal spin clusters in the 𝑆 = $\frac{1}{2}$ distorted honeycomb lattice antiferromagnet Cu 2 ⁢(pymca)⁢ 3 ⁢(ClO 4 )

Cu 2 ⁢(pymca) ⁢3 (ClO 4 ) (pymca: pyrimidine-2-carboxylate) consists of a slightly distorted honeycomb lattice of Cu 2+ spins, which shows no long-range magnetic order down to 0.6 K. A magnetization study revealed 1/3 and 2/3 plateau phases [A. Okutani et al., J. Phys. Soc. Jpn. 88, 013703 (2019)], which is not expected for regular honeycomb antiferromagnets. Inelastic neutron scattering experiments were performed using a powder sample to investigate the exchange interactions of this material. The spin excitations from the singlet ground state to the first three triplet states, predicted from the antiferromagnetic hexagonal spin cluster interacting with 3.9 meV, were observed. Using the exact diagonalization methods, the intercluster coupling was estimated from the excitation peak width to be about 20% of the intracluster interaction, which is consistent with the previously reported value. Finally, our exchange path model explains the anisotropic exchange interactions in the distorted honeycomb plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Field-tunable BKT and quantum phase transitions in spin-$\frac{1}{2}$ triangular lattice antiferromagnet

Quantum magnetism is one of the most active fields for exploring exotic phases and phase transitions. The recently synthesized Na 2 BaCo(PO 4 ) 2 (NBCP) is an ideal material incarnation of the spin-$\frac{1}{2}$ easy-axis triangular lattice antiferromagnet (TLAF). Experimental evidence shows that NBCP hosts the spin supersolid state with a giant magnetocaloric effect. Theory further predicts that magnetic fields can drive NBCP through Berezinskii-Kosterlitz-Thouless (BKT) and other richer quantum phase transitions. However, detecting these transitions is challenging, as they onset at ultralow temperatures near 60 mK and require high magnetization sensitivity. Using a newly developed gradient force magnetometer in a dilution refrigerator, we mapped the magnetic susceptibility phase diagram down to 30 mK. Our results provide a more comprehensive and accurate understanding of BKT melting of spin supersolidity and several field-tunable quantum phase transitions, which establish NBCP as a model platform for frustrated magnetism and highlight potential applications of its giant magnetocaloric effects.

BKT transition↗