Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GROUND STATES”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2022

Circuit complexity near critical points

Here, we consider the Bose–Hubbard model in two and three spatial dimensions and numerically compute the quantum circuit complexity of the ground state in the Mott insulator and superfluid phases using a mean field approximation with additional quadratic fluctuations. After mapping to a qubit system, the result is given by the complexity associated with a Bogoliubov transformation applied to the reference state taken to be the mean field ground state. In particular, the complexity has peaks at the O(2) critical points where the system can be described by a relativistic quantum field theory. Given that we use a Gaussian approximation, near criticality the numerical results agree with a free field theory calculation. To go beyond the Gaussian approximation we use general scaling arguments that imply that, as we approach the critical point t → t c , there is a non-analytic behavior in the complexity c 2 (t) of the form |c 2 (t) – c 2 (t c )| ~ |t – t c | νd , up to possible logarithmic corrections. Here d is the number of spatial dimensions and ν is the usual critical exponent for the correlation length ξ ~ |t – t c | –ν . As a check, for d = 2 this agrees with the numerical computation if we use the Gaussian critical exponent $v=\frac{1}{2}$. Finally, using AdS/CFT methods, we study higher dimensional examples and confirm this scaling argument with non-Gaussian exponent ν for strongly interacting theories that have a gravity dual.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Structural, electronic, and magnetic properties of nearly ideal J eff = 1/2 iridium halides

Heavy transition metal magnets with J e f f = 1 2 electronic ground states have attracted recent interest due to their penchant for hosting new classes of quantum spin liquids and superconductors. Unfortunately, model systems with ideal J e f f = 1 2 states are scarce due to the importance of noncubic local distortions in most candidate materials. In this work, we identify a family of iridium halide systems [i.e., K 2 Ir Cl 6 , K 2 Ir Br 6 , ( NH 4 ) 2 Ir Cl 6 , and Na 2 Ir Cl 6 · 6 ( H 2 O ) ] with Ir 4 + electronic ground states exhibiting extremely small deviations from the ideal J e f f = 1 2 limit. We also find ordered magnetic ground states for the three anhydrous systems, with single-crystal neutron diffraction on K 2 Ir Br 6 revealing type-I antiferromagnetism. This spin configuration is consistent with expectations for significant Kitaev exchange in a face-centered-cubic magnet. This work establishes that incorporating isolated Ir X 6 octahedra in materials, where X is a halogen ion with a low electronegativity, is an effective design principle for realizing unprecedented proximity to the pure J e f f = 1 2 state. However, at the same time, we highlight undeniable deviations from this ideal state, even in clean materials with ideal Ir X 6 octahedra as inferred from the global cubic crystal structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Resource-Optimized Fermionic Local-Hamiltonian Simulation on a Quantum Computer for Quantum Chemistry

The ability to simulate a fermionic system on a quantum computer is expected to revolutionize chemical engineering, materials design, nuclear physics, to name a few. Thus, optimizing the simulation circuits is of significance in harnessing the power of quantum computers. Here, we address this problem in two aspects. In the fault-tolerant regime, we optimize the R z and T gate counts along with the ancilla qubit counts required, assuming the use of a product-formula algorithm for implementation. We obtain a savings ratio of two in the gate counts and a savings ratio of eleven in the number of ancilla qubits required over the state of the art. In the pre-fault tolerant regime, we optimize the two-qubit gate counts, assuming the use of the variational quantum eigensolver (VQE) approach. Specific to the latter, we present a framework that enables bootstrapping the VQE progression towards the convergence of the ground-state energy of the fermionic system. This framework, based on perturbation theory, is capable of improving the energy estimate at each cycle of the VQE progression, by about a factor of three closer to the known ground-state energy compared to the standard VQE approach in the test-bed, classically-accessible system of the water molecule. The improved energy estimate in turn results in a commensurate level of savings of quantum resources, such as the number of qubits and quantum gates, required to be within a pre-specified tolerance from the known ground-state energy. We also explore a suite of generalized transformations of fermion to qubit operators and show that resource-requirement savings of up to more than 20 % , in small instances, is possible.

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

Theory of topological defects and textures in two-dimensional quantum orders with spontaneous symmetry breaking

In this article, we consider two-dimensional (2d) quantum many-body systems with long-range orders, where the only gapless excitations in the spectrum are Goldstone modes of spontaneously broken continuous symmetries. To understand the interplay between classical long-range order of local order parameters and quantum order of long-range entanglement in the ground states, we study the topological point defects and textures of order parameters in such systems. We show that the universal properties of point defects and textures are determined by the remnant symmetry enriched topological order in the symmetry-breaking ground states with a nonfluctuating order parameter, and provide a classification for their properties based on the inflation-restriction exact sequence. We highlight a few phenomena revealed by our theory framework. First, in the absence of intrinsic topological orders, we show a connection between the symmetry properties of point defects and textures to deconfined quantum criticality. Second, when the symmetry-breaking ground state has intrinsic topological orders, we show that the point defects can permute different anyons when braided around. They can also obey projective fusion rules in the sense that multiple vortices can fuse into an Abelian anyon, a phenomenon for which we coin “defect fractionalization.” Finally, we provide a formula to compute the fractional statistics and fractional quantum numbers carried by textures (skyrmions) in Abelian topological orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Magnetoresistance effects in a spin-fermion model for multilayers

Here, we consider a spin-fermion model consisting of free electrons coupled to classical spins where the latter are embedded in a quasi-one-dimensional superlattice structure consisting of spin blocks separated by spinless buffers. Using a spiral ansatz for the spins, we study the effect of the electron mediated Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction on the $\textit{T}$ = 0 ground state of the system. We find that the RKKY interaction can lead to ferromagnetic, antiferromagnetic, or intermediate spiral phases for different system parameters. When the width is much larger than the length of the individual blocks, the spiral phases are suppressed, and the ground state oscillates between ferromagnetic and antiferromagnetic orders as the size of the buffer regions is varied. This is accompanied by a corresponding oscillation in the Drude weight reflecting an increased conductivity in the ferromagnetic state compared to the antiferromagnetic one. These results are reminiscent of classic giant magnetoresistance phenomena observed in a similar geometry of thin sandwiched magnetic and nonmagnetic layers. Our analysis provides a robust framework for understanding the role of the RKKY interaction on the ground-state order and corresponding transport properties of such systems, extending beyond the conventional perturbative regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Iron( IV ) alkyl complexes: electronic structure contributions to Fe–C bond homolysis and migration reactions that form N–C bonds from N 2

High-valent iron alkyl complexes are rare, as they are prone to Fe–C bond homolysis. Here, we describe an unusual way to access formally iron(IV) alkyl complexes through double silylation of iron(I) alkyl dinitrogen complexes to form an NNSi 2 group. Spectroscopically validated computations show that the disilylehydrazido(2−) ligand stabilizes the formal iron(IV) oxidation state through a strongly covalent Fe–N π-interaction, in which one π-bond fits an “inverted field” description. This means that the two bonding electrons are localized more on the metal than the ligand, and thus an iron(II) resonance structure is a significant contributor, similar to the previously-reported phenyl analogue. However, in contrast to the phenyl complex which has an S = 1 ground state, the ground state of the alkyl complex is S = 2, which places one electron in the π* orbital, leading to longer and weaker Fe–N bonds. The reactivity of these hydrazido(2−) complexes is dependent on the steric and electronic properties of the specific alkyl group. When the alkyl group is the bulky trimethylsilylmethyl, the formally iron(IV) species is stable at room temperature and no migration of the alkyl ligand is observed. However, the analogous complex with the smaller methyl ligand does indeed undergo migration of the carbon-based ligand to the NNSi 2 group to form a new N–C bond. This migration is followed by isomerization of the hydrazido ligand, and the product exists as two isomers that have distinct η 1 and η 2 binding of the hydrazido group. Lastly, when the alkyl group is benzyl, the Fe–C bond homolyzes to give a three-coordinate hydrazido(2−) complex which is likely due to the greater stability of a benzyl radical compared to that for methyl or trimethylsilylmethyl. These studies demonstrate the availability of a hydrocarbyl migration pathway at formally iron(IV) centers to form new N–C bonds directly to N 2 , though product selectivity is highly dependent on the identity of the migrating group.

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

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2005

Atomic Clock Based on Opto-Electronic Oscillator

A proposed highly accurate clock or oscillator would be based on the concept of an opto-electronic oscillator (OEO) stabilized to an atomic transition. Opto-electronic oscillators, which have been described in a number of prior NASA Tech Briefs articles, generate signals at frequencies in the gigahertz range characterized by high spectral purity but not by longterm stability or accuracy. On the other hand, the signals generated by previously developed atomic clocks are characterized by long-term stability and accuracy but not by spectral purity. The proposed atomic clock would provide high spectral purity plus long-term stability and accuracy a combination of characteristics needed to realize advanced developments in communications and navigation. In addition, it should be possible to miniaturize the proposed atomic clock. When a laser beam is modulated by a microwave signal and applied to a photodetector, the electrical output of the photodetector includes a component at the microwave frequency. In atomic clocks of a type known as Raman clocks or coherent-population-trapping (CPT) clocks, microwave outputs are obtained from laser beams modulated, in each case, to create two sidebands that differ in frequency by the amount of a hyperfine transition in the ground state of atoms of an element in vapor form in a cell. The combination of these sidebands produces a transparency in the population of a higher electronic level that can be reached from either of the two ground-state hyperfine levels by absorption of a photon. The beam is transmitted through the vapor to a photodetector. The components of light scattered or transmitted by the atoms in the two hyperfine levels mix in the photodetector and thereby give rise to a signal at the hyperfine- transition frequency. The proposed atomic clock would include an OEO and a rubidium- or cesium- vapor cell operating in the CPT/Raman regime (see figure). In the OEO portion of this atomic clock, as in a typical prior OEO, a laser beam would pass through an electro-optical modulator, the modulated beam would be fed into a fiber-optic delay line, and the delayed beam would be fed to a photodetector. The electrical output of the photodetector would be detected, amplified, filtered, and fed back to the microwave input port of the modulator. The laser would be chosen to have the same wavelength as that of the pertinent ground-state/higher-state transition of the atoms in the vapor. The modulator/ filter combination would be designed to operate at the microwave frequency of the hyperfine transition. Part of the laser beam would be tapped from the fiberoptic loop of the OEO and introduced into the vapor cell. After passing through the cell, this portion of the beam would be detected differentially with a tapped portion of the fiber-optically-delayed beam. The electrical output of the photodetector would be amplified and filtered in a loop that would control a DC bias applied to the modulator. In this manner, the long-term stability and accuracy of the atomic transition would be transferred to the OEO.

Maleki, Lute↗
DOE OSTI2024

GdWN 3 is a nitride perovskite

Nitride perovskites ABN3 are an emerging and highly underexplored class of materials that are of interest due to their intriguing calculated ferroelectric, optoelectronic, and other functional properties. Incorporating novel A-site cations is one strategy to tune and expand such properties; for example, Gd 3+ is compelling due to its large magnetic moment, potentially leading to multiferroic behavior. However, the theoretically predicted ground state of GdWN 3 was a non-perovskite monoclinic structure. Here, we experimentally show that GdWN 3-y crystallizes in a perovskite structure. High-throughput combinatorial sputtering with activated nitrogen is employed to synthesize thin films of Gd 2-x W x N 3-y O y with oxygen content y < 0.05. Ex situ annealing crystallizes a polycrystalline perovskite phase in a narrow composition window near x = 1. LeBail fits of synchrotron grazing incidence wide angle x-ray scattering data are consistent with a perovskite ground-state structure. Refined density functional theory calculations that included antiferromagnetic configurations confirm that the ground-state structure of GdWN 3 is a distorted Pnma perovskite with antiferromagnetic ordering, in contrast to prior predictions. Initial property measurements find that GdWN 3-y is paramagnetic down to T = 2 K with antiferromagnetic correlations and that the absorption onset depends on cation stoichiometry. This work provides an important path toward both the rapid expansion of the emerging family of nitride perovskites and understanding their potential multiferroic properties.

14 SOLAR ENERGY↗
NASA NTRSDate not supplied

Self-consistent Quantum Iteratively Sparsified Hamiltonian Algorithm (SQuISH)

Due to coherence time limitations, reducing the resources required to run quantum algorithms and simulate physical systems on a quantum computer is crucial. With regards to Hamiltonian simulation, a significant effort has focused on building efficient algorithms using various factorizations and truncations, typically derived from the Hamiltonian alone. We introduce a new paradigm for improving Hamiltonian simulation and reducing the cost of ground state problems based on ideas recently developed for classical chemistry simulations. The key idea is that one can find efficient ways to reduce resources needed by quantum algorithms by making use of two key pieces of information: the Hamiltonian operator and an approximate ground state wavefunction. We refer to our algorithm as the self-consistent quantum iteratively sparsified Hamiltonian (SQuISH). By performing our scheme iteratively, one can drive SQuISH to create an accurate wavefunction using a truncated, resource-efficient Hamiltonian. By utilizing this more compact Hamiltonian, our algorithm provides an approach to reduce the gate complexity of ground state calculations on quantum hardware. As proof of principle, we implement SQuISH using configuration interaction for small molecules and coupled cluster for larger systems. Through our combination of approaches, we demonstrate how it performs on a range of systems, the largest of which would require more than 200 qubits to run on quantum hardware.

Diana Chamaki↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2021

Wigner SU(4) symmetry, clustering, and the spectrum of $$^{12}$$C

Abstract We present lattice calculations of the low-lying spectrum of $$^{12}$$ 12 C using a simple nucleon–nucleon interaction that is independent of spin and isospin and therefore invariant under Wigner’s SU(4) symmetry. We find strong signals for all excited states up to $$\sim 15$$ ∼ 15 MeV above the ground state, and explore the structure of each state using a large variety of $$\alpha $$ α cluster and harmonic oscillator trial states, projected onto given irreducible representations of the cubic group. We are able to verify earlier findings for the $$\alpha $$ α clustering in the Hoyle state and the second $$2^+$$ 2 + state of $$^{12}$$ 12 C. The success of these calculations to describe the full low-lying energy spectrum using spin-independent interactions suggest that either the spin-orbit interactions are somewhat weak in the $$^{12}$$ 12 C system, or the effects of $$\alpha $$ α clustering are diminishing their influence. This is in agreement with previous findings from ab initio shell model calculations.

Shen, Shihang (ORCID:0000000280516466)↗
DOE OSTI2024

Ultrafast electronic relaxation pathways of the molecular photoswitch quadricyclane

Abstract The light-induced ultrafast switching between molecular isomers norbornadiene and quadricyclane can reversibly store and release a substantial amount of chemical energy. Prior work observed signatures of ultrafast molecular dynamics in both isomers upon ultraviolet excitation but could not follow the electronic relaxation all the way back to the ground state experimentally. Here we study the electronic relaxation of quadricyclane after exciting in the ultraviolet (201 nanometres) using time-resolved gas-phase extreme ultraviolet photoelectron spectroscopy combined with non-adiabatic molecular dynamics simulations. We identify two competing pathways by which electronically excited quadricyclane molecules relax to the electronic ground state. The fast pathway (<100 femtoseconds) is distinguished by effective coupling to valence electronic states, while the slow pathway involves initial motions across Rydberg states and takes several hundred femtoseconds. Both pathways facilitate interconversion between the two isomers, albeit on different timescales, and we predict that the branching ratio of norbornadiene/quadricyclane products immediately after returning to the electronic ground state is approximately 3:2.

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

Filling constraints on translation invariant dipole conserving systems

Systems with conserved dipole moment have drawn considerable interest in light of their realization in recent experiments on tilted optical lattices. Here, an important question for such systems is delineating the conditions under which they admit a unique gapped ground state that is consistent with all symmetries. Here, we study one-dimensional translation-invariant lattices that conserve U(1) charge and Z L dipole moment, where discreteness of the dipole symmetry is enforced by periodic boundary conditions, with L the system size. We show that in these systems, a symmetric, gapped, and non-degenerate ground state requires not only integer charge filling, but also a fixed value of the dipole filling, while other fractional dipole fillings enforce either a gapless or symmetry-breaking ground state. In contrast with prior results in the literature, we find that the dipole filling constraint depends both on the charge filling as well as the system size, emphasizing the subtle interplay of dipole symmetry with boundary conditions. We support our results with numerical simulations and exact results.

1-dimensional systems↗
NASA NTRS1987

Detection of interstellar CH in the far-infrared

The first astronomical detection of CH in the far-infrared has been made. A ground state of rotational transition was observed in absorption against the far-infrared continuum peak of Sgr B2. The lines are resolved at a velocity resolution of 62 km/s, have a line width of roughly 250 km/s, and a line center optical depth of about 0.29. The inferred total column density of CH in the ground state along the line of sight is roughly 1.6 x 10 to the 15th/sq cm. Comparison of the far-infrared profiles to the 3 GHz emission lines confirms that the ground-state Lambda-doublet levels are inverted and gives an accurate estimate of the excitation temperature. The excitation temperature of the 3264 MHz line varies from cloud to cloud along the line of sight, the levels being most inverted in the Sgr B2 molecular cloud. The large intensity of the 3264 MHz line in this cloud relative to other clouds along the line of sight may thus be primarily an excitation effect.

Stacey, Gordon J.↗
DOE OSTI2025

Topological and magnetic properties of a noncollinear spin state on a honeycomb lattice in a magnetic field

Here, this paper studies the topological and magnetic properties of a noncollinear spin state on a honeycomb lattice that evolves from coplanar to ferromagnetic with a magnetic field applied along the z axis. The coplanar state is stabilized by nearest-neighbor ferromagnetic interactions, single-ion anisotropy along z, and DzyaloshinskiiMoriya interactions between next-nearest-neighbor sites. Below the critical field H$_c$ that aligns the spins, the magnetic unit cell contains six sites and the spin dynamics contains six magnon modes. Although the classical energy is degenerate with respect to the twist angle φ between nearest-neighbor spins, the dependence of the free energy on φ at low temperatures is dominated by the magnon zero-point energy, which contains extremum at φ = πl/3 for integer l. The only unique ground states GS(φ) have l = 0 or 1. For H < H$'_c$, the zero-point energy has minima at even l and the ground state is GS(0); for H$'_c$ < H < H$'_c$, the zero-point energy has minima at odd l and the ground state is GS(π/3). In GS(0), the magnon density of states exhibits five distinct topological phases with increasing field associated with the opening and closing of energy gaps between two or three magnonic bands. While the Berry curvature vanishes for the coplanar φ = 0 phase in zero field, the Berry curvature and Chern numbers exhibit signatures of the five topological phases below H$'_c$. Whereas the Berry curvature and Chern number are sensitive to changes in the magnon density of states within GS(π/3), the inelastic spectrum S(k,ω) is sensitive to changes in the intensity of the magnon modes in the different magnetic phases GS(0) and GS(π/3) rather than the five topological phases within GS(π/3).

Fishman, Randy S. [Oak Ridge National Laboratory (↗
DOE OSTI2025

Magnetic dilution in the triangular lattice antiferromagnet NaYb 1−𝑥 ⁢Lu 𝑥 ⁢O 2

The delafossitelike compound NaYbO 2 hosts a triangular lattice of Yb 3+ moments and is a promising candidate for the realization of a quantum spin liquid ground state—an exotic, quantum-disordered magnetic phase featuring long-range entanglement of spins. Tuning this system away from this quantum-disordered regime toward classical order or spin freezing is a powerful approach to shed light on the nature of the parent ground state. Here we leverage the substitution of nonmagnetic Lu 3+ onto the Yb 3+ sites to study the effects of magnetic disorder in NaYbO 2 using low-temperature ac susceptibility, heat capacity, and muon spin relaxation (𝜇⁢SR) measurements. Here, our 𝜇⁢SR measurements reveal resilient, correlated magnetic fluctuations that persist to at least 15% dilution, precluding conventional spin freezing and magnetic inhomogeneity. Heat capacity and magnetic susceptibility resolve a rapid suppression of the field-induced “up-up-down” magnetic order upon dilution and a crossover in the power-law behavior of the low-temperature magnetic excitations associated with the zero-field quantum disordered ground state. Taken together, these results support the notion of a robust network of entangled moments in NaYbO 2 and provide experimental validation of several models of a Heisenberg triangular lattice antiferromagnet in the presence of disorder.

Muon spin relaxation & rotation↗