SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Melting and freezing of a skyrmion lattice

We report comprehensive Monte–Carlo studies of the melting of skyrmion lattices (SkL) in systems of small, medium, and large sizes with the number of skyrmions ranging from 10 3 to over 10 5 . Large systems exhibit hysteresis similar to that observed in real experiments on the melting of SkLs. For sufficiently small systems which achieve thermal equilibrium, a fully reversible sharp solid–liquid transition on temperature with no intermediate hexatic phase is observed. A similar behavior is found on changing the magnetic field that provides the control of pressure in the SkL. Here, we find that on heating the melting transition occurs via a formation of grains with different orientations of hexagonal axes. On cooling, the fluctuating grains coalesce into larger clusters until a uniform orientation of hexagonal axes is slowly established. The observed scenario is caused by collective effects involving defects and is more complex than a simple picture of a transition driven by the unbinding and annihilation of dislocation and disclination pairs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Probing physics beyond the standard model: limits from BBN and the CMB independently and combined

We present new Big Bang Nucleosynthesis (BBN) limits on the cosmic expansion rate or relativistic energy density, quantified via the number N ν of equivalent neutrino species. We use the latest light element observations, neutron mean lifetime, and update our evaluation for the nuclear rates d + d → 3 He + n and d + d → 3 H+p. Combining this result with the independent constraints from the cosmic microwave background (CMB) yields tight limits on new physics that perturbs N ν and η prior to cosmic nucleosynthesis: a joint BBN+CMB analysis gives N ν = 2.898 ± 0.141, resulting in N ν < 3.180 at 2σ. We apply these limits to a wide variety of new physics scenarios including right-handed neutrinos, dark radiation, and a stochastic gravitational wave background. The strength of the independent BBN and CMB constraints now opens a new window: we can search for limits on potential changes in N ν and/or the baryon-to-photon ratio η between the two epochs. The present data place strong constraints on the allowed changes in N ν between BBN and CMB decoupling; for example, we find -0.708 < N ν CMB - N ν BBN < 0.328 in the case where η and the primordial helium mass fraction Yp are unchanged between the two epochs; we also give limits on the allowed variations in η or in (η, N ν ) jointly. We discuss scenarios in which such changes could occur, and show that BBN+CMB results combine to place important constraints on some early dark energy models to explain the H 0 tension. Looking to the future, we forecast the tightened precision for N ν arising from both CMB Stage 4 measurements as well as improvements in astronomical 4 He measurements. Here, we find that CMB-S4 combined with present BBN and light element observation precision can give σ(N ν ) ≃ 0.03. Such future precision would reveal the expected effect of neutrino heating (N eff -3 = 0.044) of the CMB during BBN, and would be near the level to reveal any particle species ever in thermal equilibrium with the standard model. Improved Y p measurements can push this precision even further.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Recycled dark matter

Here, we outline a new production mechanism for dark matter that we dub “recycling”: dark sector particles are kinematically trapped in the false vacuum during a dark phase transition; the false pockets collapse into primordial black holes (PBHs), which ultimately evaporate before Big Bang Nucleosynthesis (BBN) to reproduce the dark sector particles. The requirement that all PBHs evaporate prior to BBN necessitates high scale phase transitions and hence high scale masses for the dark sector particles in the true vacuum. Our mechanism is therefore particularly suited for the production of ultra heavy dark matter (UHDM) with masses above ~10 12 GeV. The correct relic density of UHDM is obtained because of the exponential suppression of the false pocket number density. Recycled UHDM has several novel features: the dark sector today consists of multiple decoupled species that were once in thermal equilibrium and the PBH formation stage has extended mass functions whose shape can be controlled by IR operators coupling the dark and visible sectors.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Hydrodynamic description of D meson production in high-energy heavy-ion collisions *

The large values and constituent-quark-number scaling of the elliptic flow of low-p T D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low-p T charm quarks, we compare the D 0 meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. D 0 mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D 0 mesons at P T < 3 GeV/c as measured in Au+Au collisions at $sqrts_{NN} = 200$ GeV. Though the experimental data on spectra are consistent with the hydrodynamic result at small p T ~ 1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum, GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Finally, our comparative study indicates that charm quarks only approach local thermal equilibrium at small p T , even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate p T .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Evidence for sub-Chandrasekhar Type Ia supernovae from the last major merger

ABSTRACT We investigate the contribution of sub-Chandrasekhar mass Type Ia supernovae to the chemical enrichment of the Gaia Sausage galaxy, the progenitor of a significant merger event in the early life of the Milky Way. Using a combination of data from Nissen & Schuster, the GALactic Archaeology with HERMES (GALAH) Data Release 3 [with 1D non-local thermal equilibrium (NLTE) abundance corrections], and the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we fit analytic chemical evolution models to a nine-dimensional chemical abundance space (Fe, Mg, Si, Ca, Cr, Mn, Ni, Cu, and Zn) in particular focusing on the iron-peak elements, Mn and Ni. We find that low [Mn/Fe] $\sim -0.15\, \mathrm{dex}$ and low [Ni/Fe] $\sim -0.3\, \mathrm{dex}$ Type Ia yields are required to explain the observed trends beyond the [α/Fe] knee of the Gaia Sausage (approximately at [Fe/H] $=-1.4\, \mathrm{dex}$). Comparison to theoretical yield calculations indicates a significant contribution from sub-Chandrasekhar mass Type Ia supernovae in this system (from ${\sim} 60\, \mathrm{per\, cent}$ to $100\, \mathrm{per\, cent}$ depending on the theoretical model with an additional ${\pm} 10\, \mathrm{per\, cent}$ systematic from NLTE corrections). We compare to results from other Local Group environments including dwarf spheroidal galaxies, the Magellanic Clouds, and the Milky Way’s bulge, finding the Type Ia [Mn/Fe] yield must be metallicity dependent. Our results suggest that sub-Chandrasekhar mass channels are a significant, perhaps even dominant, contribution to Type Ia supernovae in metal-poor systems, whilst more metal-rich systems could be explained by metallicity-dependent sub-Chandrasekhar mass yields, possibly with additional progenitor mass variation related to star formation history, or an increased contribution from Chandrasekhar mass channels at higher metallicity.

Sanders, Jason L.↗

DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗

DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗

DOE OSTI2024

Angular correlations of cosmic microwave background spectrum distortions from photon diffusion

ABSTRACT During cosmic recombination, charged particles bind into neutral atoms and the mean free path of photons rapidly increases, resulting in the familiar diffusion damping of primordial radiation temperature variations. An additional effect is a small photon spectrum distortion, because photons arriving from a particular sky direction were originally in thermal equilibrium at various spatial locations with different temperatures; the combination of these different blackbody temperature distributions results in a spectrum with a Compton y-distortion. Using the approximation that photons had zero mean free path prior to their second-to-last scattering, we derive an expression for the resulting y-distortion, and compute the angular correlation function of the diffusion y-distortion and its cross-correlation with the square of the photon temperature fluctuation. Detection of the cross-correlation is within reach of existing arcminute-resolution microwave background experiments such as the Atacama Cosmology Telescope and the South Pole Telescope.

Astronomy & Astrophysics↗

DOE OSTI2025

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, introducing a bulk viscosity into the plasma via the internal degrees of freedom of the molecular composition, directly impacting the decay of compressible modes, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$. This is especially important for small-scale, turbulent dynamo processes in the interstellar medium (ISM), which are known to be sensitive to the effects of compression. To control the viscous properties of $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$, we perform trans-sonic, visco-resistive dynamo simulations with additional bulk viscosity $\nu _{\text{bulk}}$, deriving a new $\nu _{\text{bulk}}$ Reynolds number $\text{Re}_{\text{bulk}}$, and viscous Prandtl number $\text{P}\nu \equiv \text{Re}_{\text{bulk}}/ \text{Re}_{\text{shear}}$, where $\text{Re}_{\text{shear}}$ is the shear viscosity Reynolds number. We derive a framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\boldsymbol {\nabla }\otimes \mathrm{{\boldsymbol {\mathit {v}}}}$, where $\mathrm{{\boldsymbol {\mathit {v}}}}$ is the fluid velocity. We find that $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in compressible dynamos. In the absence of $\nu _{\text{bulk}}$ ($\text{P}\nu \rightarrow \infty$), $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ pile up on small-scales, creating a spectral bottleneck, which disappears for $\text{P}\nu \approx 1$. As $\text{P}\nu$ decreases, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between high-k modes. We emphasize the importance of further understanding the role of $\nu _{\text{bulk}}$ in compressible astrophysical plasmas, which we estimate could be as strong as the shear viscosity in the cold ISM, and highlight that compressible direct numerical simulations without bulk viscosity have unresolved compressible mode dissipation scales.

MHD↗

DOE OSTI2026

Nanoscale Observation and Control of Quasiparticle Induced Magnetic Noise in a Superconducting Resonator

Superconducting circuits are arguably taking a leading role in driving the ongoing quantum technological revolution. A detailed knowledge of the microscopic fluctuating electromagnetic properties plays an important role in advancing the circuitry design, testing, and material integration of cutting-edge superconducting quantum electronics. Here, in this work, we report scanning nitrogen-vacancy (NV) quantum sensing of local magnetic noise environment of an on-chip superconducting resonator. We find that quasiparticle-induced fluctuating magnetic fields can drive NV spin relaxation, which shows a peak value around the superconducting transition point of niobium at the thermal equilibrium state. External microwave driving at the resonator mode frequency significantly increases the quasiparticle density, leading to enhancement of magnetic noise. We further perform optically detected magnetic resonance measurements to demonstrate quasiparticle magnetic noise mediated off-resonant dipole coupling between the NV center and niobium resonator. Our Letter reports experimental observation of the Hebel-Slichter peak signature by an external sensor outside of a superconductor. The presented study also highlights the advantages of quantum sensors in investigating miniaturized superconducting devices, providing insights into their future performance improvements.

Li, Senlei [Georgia Institute of Technology]↗

DOE OSTI2026

Rotating synchrotron radiation: Photon emission from magnetized and rotating quark-gluon plasma

This paper investigates the production of nonprompt photons originating from rotating synchrotron radiation (RoSyRa), specifically the emission of photons by a rigidly rotating quark-gluon plasma in thermal equilibrium, in the presence of an external magnetic field. We compute the nonprompt photon spectrum and its elliptic flow (𝑣 2 ) at midrapidity. In particular, we investigate the finite-volume effects. We find that at low transverse momentum, the magnetic field induces a significant 𝑣 2 , while the plasma rotation boosts the synchrotron radiation of negatively charged quarks. These findings make RoSyRa a viable candidate mechanism to resolve the “direct photon puzzle.”

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Diffusion of acceptor dopants in monoclinic 𝛽−Ga 2⁢ O 3

𝛽−Ga 2 ⁢O 3 is a promising material for next-generation power electronics because of its ultrawide band gap and high critical breakdown voltage. However, realizing its full potential requires precise control over dopant incorporation and stability. In this work, we use first-principles calculations to systematically assess the diffusion behavior of eight potential deep-level substitutional acceptors (Au, Ca, Co, Cu, Fe, Mg, Mn, and Ni) in 𝛽−Ga 2 ⁢O 3 . We consider two key diffusion mechanisms: (i) interstitial diffusion under nonequilibrium conditions relevant to ion implantation, and (ii) trap-limited diffusion (TLD) under near-equilibrium thermal annealing conditions. Our results reveal a strong diffusion anisotropy along the 𝑏 and 𝑐 axes, with dopant behavior governed by competition between diffusion and incorporation (or dissociation) activation energies. Under interstitial diffusion, Ca$^{2+}_{i}$ and Mg$^{2+}_{i}$ show the most favorable combination of low migration and incorporation barriers, making them promising candidates for efficient doping along the 𝑏 and 𝑐 axes, respectively. In contrast, Au$^{+}_{i}$ diffuses readily, but exhibits an incorporation barrier that exceeds 5 eV, rendering it ineffective as a dopant. From a thermal stability perspective, Co$^{2+}_{i}$ shows poor activation but high diffusion barriers, which may suppress undesirable migration at elevated temperatures. Under trap-limited diffusion, the dissociation of dopant-host complexes controls mobility. Mg$^{2+}_{i}$ again emerges as a leading candidate, exhibiting the lowest dissociation barriers along both axes, whereas Co$^{2+}_{i}$ and Fe$^{2+}_{i}$ display the highest barriers, suggesting improved dopant retention under thermal stress. In conclusion, our findings guide dopant selection by balancing activation and thermal stability, essential for robust semi-insulating substrates.

Defects↗

DOE OSTI2025

Hyperstealth dark matter and long-lived particles

A new dark matter candidate is proposed that arises as the lightest baryon from a confining 𝑆⁢𝑈⁡(𝑁) gauge theory which equilibrates with the Standard Model only through electroweak interactions. Surprisingly, this candidate can be as light as a few GeV. The lower bound arises from the intersection of two competing requirements: (i) the equilibration sector of the model must be sufficiently heavy, at least several TeV, to avoid bounds from colliders, and (ii) the lightest dark meson (that may be the dark 𝜂′, 𝜎, or the lightest glueball) has suppressed interactions with the SM and must decay before big bang nucleosynthesis. The low-energy dark sector consists of one flavor that is electrically neutral and an almost electroweak singlet. The dark matter candidate is the lightest baryon consisting of 𝑁 of these light flavors leading to a highly suppressed elastic scattering rate with the Standard Model (SM). The equilibration sector consists of vectorlike dark quarks that transform under the electroweak group, ensuring that the dark sector can reach thermal equilibrium with the SM in the early Universe. The lightest dark meson lifetimes vary between 10 −3 ≲ 𝑐⁢𝜏 ≲ 10 7 m, providing an outstanding target for LHC production and experimental detection. We delineate the interplay between the lifetime of the light mesons, the suppressed direct detection cross section of the lightest baryon, and the scale of equilibration sector that can be probed at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Correlations between electronic order and structural distortions and their ultrafast dynamics in the single-layer manganite P r 0.5 C a 1.5 Mn O 4

Time-resolved x-ray diffraction experiments on the half-doped single-layered manganite Pr 0.5 Ca 1.5 Mn O 4 are used to monitor the ultrafast photoinduced dynamics of the structural distortion associated with the charge and orbital ordering (CO/OO). As in the nonlayered three-dimensional counterpart, the ordered phase melts in less than 100 fs after 800-nm photoexcitation and subsequently partially recovers due to thermal equilibration of electronic and vibrational systems. Photoexciting Pr 0.5 Ca 1.5 Mn O 4 below the transition temperature of a second structural phase transition that occurs around 146 K (deep inside the CO/OO phase) releases this structural transition, but progresses on a much slower timescale. Finally, this additional reduction of crystal symmetry, which we ascribe to a further tilt of the oxygen octahedra, can thus be considered to be only weakly coupled to CO/OO. Furthermore, static hard-x-ray and resonant soft-x-ray diffraction at the Mn L 2,3 edges experiments identify correlations between structural distortions and electronic order in thermal equilibrium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Relaxation of the degenerate one-dimensional Fermi gas

In this work, we study how a system of one-dimensional spin-1/2 fermions at temperatures well below the Fermi energy approaches thermal equilibrium. The interactions between fermions are assumed to be weak and are accounted for within the perturbation theory. In the absence of an external magnetic field, spin degeneracy strongly affects relaxation of the Fermi gas. For sufficiently short-range interactions, the rate of relaxation scales linearly with temperature. Focusing on the case of the system near equilibrium, we linearize the collision integral and find the exact solution of the resulting relaxation problem. Furthermore, we discuss the application of our results to the evaluation of the transport coefficients of the one-dimensional Fermi gas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Landau-Lifshitz-Bloch equation for ferrimagnets with higher-order interaction

In this work, we present a micromagnetic formulation for modeling the magnetization dynamics and thermal equilibrium in ferrimagnetic materials at low and elevated temperatures. The formulation is based on a mean-field approximation (MFA). In this formulation, the ferrimagnet is described micromagnetically by two coupled sublattices with corresponding interactions, including inter- and intrasublattice micromagnetic exchange as well as four-spin interactions described as an intersublattice molecular field with a cubic dependence of the magnetization. The MFA is used to derive a Landau-Lifshitz-Bloch type equation for ferrimagnetic material, including cases with a ferromagnetic-antiferromagnetic phase transitions. For validation, the results obtained via the presented model are compared with recent experimental data for phase transitions in FeRh.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Orbital magnetization of Floquet topological systems

A general expression for the orbital magnetization of a Floquet system is derived. The expression holds for a clean system, and is valid for any driving protocol, and arbitrary occupation of the bands. The orbital magnetization is shown to be large not only for Chern insulators, but also for anomalous phases where the Chern number does not fully account for the topology. In addition, the orbital magnetization is shown to take significant values both for a thermal equilibrium occupation of the Floquet bands, and for occupations determined by a quantum quench from an initial state with zero orbital magnetization. For the latter case, the orbital magnetization is shown to be highly sensitive to van Hove singularities of the Floquet bands.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗