Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Freezing”

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 541 records · Page 30

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗

Global Importance of Secondary Ice Production

Measured ice crystal number concentrations are often orders of magnitude higher than the number concentrations of ice nucleating particles, indicating the existence of secondary ice production (SIP) in clouds. Here, we present the first study to examine the global importance of SIP through the droplet shattering during freezing of rain, ice-ice collision fragmentation, and rime splintering, using a global climate model. Our results show that SIP happens quite uniformly in the two hemispheres and dominates the ice formation in the moderately cold clouds with temperatures warmer than –15°C. SIP decreases the global annual average liquid water path by –14.6 g m –2 (–22%), increases the ice water path by 8.7 g m –2 (23%), improving the model agreement with observations. SIP changes the global annual average shortwave, longwave, and net cloud forcing by 2.1, –1.0, and 1.1 W m –2 , respectively, highlighting the importance of SIP on cloud properties on the global scale.

54 ENVIRONMENTAL SCIENCES↗

Invigoration or Enervation of Convective Clouds by Aerosols?

Are convective clouds strengthened by the addition of aerosol particles? We present new theoretical calculations which starkly contrast previous results. Prior foundational work suggested that aerosols strongly invigorate convective cloud updrafts via changes to cold-phase processes. We show that the peak magnitude of invigoration by this mechanism is substantially reduced for cold-based storms. For warm-based storms, the updrafts are weakened, not strengthened, by aerosol-induced changes to cold-phase processes. Our calculations show that if invigoration occurs, it is driven primarily by changes to warm-phase processes and may be largest well before the storm reaches maturity. The calculations are based on a new formulation of the moist adiabatic lapse rate that accounts for freezing, supersaturation, and condensate loading. The results significantly reshape our understanding of the impact of aerosols on convective updrafts.

58 GEOSCIENCES↗

Mixed‐Phase Clouds Over the Southern Ocean as Observed From Satellite and Surface Based Lidar and Radar

Abstract This study investigates the occurrence of mixed‐phase clouds (MPC, i.e., cloud layers containing both liquid and ice water at sub‐freezing temperatures) over the Southern Ocean (SO) using space‐ and surface‐based lidar and radar observations. The occurrence of supercooled clouds is dominated by geometrically thin (<1 km) layers that rarely contain ice. We diagnose layers that are geometrically thicker than 1 km to contain ice ~65% and ~4% of the time from below by surface remote sensors and from above by orbiting remote sensors, respectively. We examine the discrepancy in MPC occurrence statistics as diagnosed from below and above the cloud layer. From above, we find that MPC occurrence has a gradient associated with the Antarctic Polar Front near 55°S with a rare occurrence of satellite‐derived MPC south of that latitude. In contrast, surface sensors find ice in 33% of supercooled liquid water layers. We infer using observing system simulation experiments and data analysis that space‐based lidar cannot identify the occurrence of MPC except when secondary ice‐forming processes operate in convection that is, sufficiently strong to loft ice crystals to cloud tops. We conclude that the CALIPSO phase statistics of MPC have a severe low bias in MPC occurrence. Based on surface‐based statistics in the SO, we present a parameterization of the frequency of MPC as a function of cloud top temperature that differs substantially from that used in recent climate model simulations.

54 ENVIRONMENTAL SCIENCES↗

Climatology of Melting Layer Heights Estimated From Cloud Radar Observations at Various Locations

A melting layer (ML) detection algorithm for cloud radar with polarimetric capability was developed and applied to the cloud radar data collected from five different sites around the world for several years. The retrieved melting layer top height (MLH) showed a very good correspondence with the ECMWF Reanalysis 5 zero degree level data for all five sites. The ML characteristics were distinctively different for different sites, revealing climatological characteristics of ML forming clouds in different regions. Generally, ML tended to occur more frequently in summer than in winter except for a maritime site, where low stratiform clouds formed frequently in summer, the top of which might be lower than the freezing level. In contrast, at two Arctic sites, ML occurred almost exclusively in summer because it was too cold to have an ML in the other seasons. The MLH also varied significantly from site to site but generally was higher during warmer seasons. Based on MLH, two new indices, bulk temperature lapse rate (BLR) and relative depth (RD) of liquid cloud below MLH (RD) were developed, which were useful to explain the environmental characteristics of the five sites. BLR generally increased with the surface temperature at all sites except at the marine site that showed an opposite trend, where a unique synoptic pattern in winter generated high BLR in this cold season. Finally, these findings confirm that studies on thermodynamic structures using cloud radars can be broadened, taking advantage of BLR and RD information, as these indices can represent environmental thermodynamic characteristics of the clouds that have ML.

54 ENVIRONMENTAL SCIENCES↗

Ice Nucleating Particle Connections to Regional Argentinian Land Surface Emissions and Weather During the Cloud, Aerosol, and Complex Terrain Interactions Experiment

Abstract Here, we present a multi‐season study of ice‐nucleating particles (INPs) active via the immersion freezing mechanism, which took place in north‐central Argentina, a worldwide hotspot for mesoscale convective storms. INPs were measured untreated, after heating to 95°C, and after hydrogen peroxide digestion. No seasonal cycle of INP concentrations was observed. Heat labile INPs, which we define as “biological” herein, dominated the population active at −5 to −20°C, while non‐heat‐labile organic INPs (decomposed by peroxide) dominated at lower temperatures, from −20 to −28°C. Inorganic INPs (remaining after peroxide digestion), were minor contributors to the overall INP activity. Biological INP concentration active around −12°C peaked during rain events and under high relative humidity, reflecting emission mechanisms independent of the background aerosol concentration. The ratio of non‐heat‐labile organic and inorganic INPs was generally constant, suggesting they originated from the same source, presumably from regional arable topsoil based on air mass histories. Single particle mass spectrometry showed that soil particles aerosolized from a regionally common agricultural topsoil contained known mineral INP sources (K‐feldspar and illite) as well as a significant organic component. The INP activity observed in this study correlates well with agricultural soil INP activities from this and other regions of the world, suggesting that the observed INP spectra might be typical of many arable landscapes. These results demonstrate the strong influence of regional continental landscapes, emitting INPs of types that are not yet well represented in global models.

54 ENVIRONMENTAL SCIENCES↗

Influence of Tundra Polygon Type and Climate Variability on CO 2 and CH 4 Fluxes Near Utqiagvik, Alaska

Arctic tundra has the potential to generate significant climate feedbacks, but spatial complexity makes it difficult to quantify the impacts of climate on ecosystem-atmosphere fluxes, particularly in polygonal tundra comprising wetter and drier polygon types on the scale of tens of meters. We measured CO 2 , CH 4 , and energy fluxes using eddy covariance for 7 yr (April to November, 2013–2019) in polygonal tundra near Utqiagvik, Alaska. This period saw the earliest snowmelt, latest snow accumulation, and hottest summer on record. To estimate fluxes by polygon type, we combined a polygon classification with a flux-footprint model. Methane fluxes were highest in the summer months but were also large during freeze-up and increased with the warming trend in August–November temperatures. While CO 2 respiration had a consistent, exponential relationship with temperature, net ecosystem exchange was more variable among years. CO 2 and CH 4 exchange (June–September) ranged between -0.83 (Standard error [SE] = 0.03) and -1.32 (SE = 0.04) μmol m -2 s-1 and 13.92 (SE = 0.26)—23.42 (SE = 0.45) nmol m -2 s -1 , respectively, and varied interannually (p ≤ 0.05). The maximum-influence method effectively attributed fluxes to polygon types. Areas dominated by low-centered polygons had higher CO 2 fluxes except in 2016–2017. Methane fluxes were highest in low-centered polygons 2013–2015 and in flat-centered polygons in subsequent years, possibly due to increasing temperature and precipitation. Sensible and latent heat fluxes also varied significantly among polygon types. Accurate characterization of Arctic fluxes and their climate dependencies requires spatial disaggregation and long term observations.

54 ENVIRONMENTAL SCIENCES↗

Ice‐Nucleating Particles That Impact Clouds and Climate: Observational and Modeling Research Needs

Abstract Atmospheric ice‐nucleating particles (INPs) play a critical role in cloud freezing processes, with important implications for precipitation formation and cloud radiative properties, and thus for weather and climate. Additionally, INP emissions respond to changes in the Earth System and climate, for example, desertification, agricultural practices, and fires, and therefore may introduce climate feedbacks that are still poorly understood. As knowledge of the nature and origins of INPs has advanced, regional and global weather, climate, and Earth system models have increasingly begun to link cloud ice processes to model‐simulated aerosol abundance and types. While these recent advances are exciting, coupling cloud processes to simulated aerosol also makes cloud physics simulations increasingly susceptible to uncertainties in simulation of INPs, which are still poorly constrained by observations. Advancing the predictability of INP abundance with reasonable spatiotemporal resolution will require an increased focus on research that bridges the measurement and modeling communities. This review summarizes the current state of knowledge and identifies critical knowledge gaps from both observational and modeling perspectives. In particular, we emphasize needs in two key areas: (a) observational closure between aerosol and INP quantities and (b) skillful simulation of INPs within existing weather and climate models. We discuss the state of knowledge on various INP particle types and briefly discuss the challenges faced in understanding the cloud impacts of INPs with present‐day models. Finally, we identify priority research directions for both observations and models to improve understanding of INPs and their interactions with the Earth System.

54 ENVIRONMENTAL SCIENCES↗

Ocean Variability Beneath the Filchner-Ronne Ice Shelf Inferred From Basal Melt Rate Time Series

Fourteen phase-sensitive radars (ApRES) were deployed on the Filchner-Ronne Ice Shelf (FRIS) to measure variability in its basal melt rate. Melt rates from sites along the Ronne Depression vary seasonally, consistent with the dynamics of the propagation of seasonal dense water from the western ice front into the cavity. Several sites at the back of the FRIS cavity feature a signal with two seasonal maxima. Sub-ice shelf oceanographic data available from one of the sites indicate that this signal is caused by two different pathways followed by the same source water. Inter-annual variability is strongest along a direct flow pathway between western Ronne Ice Front and western Berkner Island. Highest melting occurred in 1999 and 2018, following anomalously low summer sea-ice concentrations in front of the ice shelf. Inter-annual melt rate variability at the back of the FRIS cavity is limited. If present, it is expressed as a suppression or delay in the arrival of the seasonal melt rate minimum, which can be understood in terms of inter-annual stratification changes and variable inflow pathways toward the sites. Long term mean ApRES melt rates agree with estimates from satellite data over eastern FRIS. However, the satellite estimates overstate the area of active basal freezing in the western part of the ice shelf. Furthermore, the temporal melt rate variability from the satellite estimates exaggerates the range of variability at both seasonal and inter-annual time scales with any correspondence between the in-situ and remotely derived inter-annual variability being limited to a single site.

54 ENVIRONMENTAL SCIENCES↗

Permafrost Carbon: Progress on Understanding Stocks and Fluxes Across Northern Terrestrial Ecosystems

Significant progress in permafrost carbon science made over the past decades include the identification of vast permafrost carbon stocks, the development of new pan-Arctic permafrost maps, an increase in terrestrial measurement sites for CO 2 and methane fluxes, and important factors affecting carbon cycling, including vegetation changes, periods of soil freezing and thawing, wildfire, and other disturbance events. Process-based modeling studies now include key elements of permafrost carbon cycling and advances in statistical modeling and inverse modeling enhance understanding of permafrost region C budgets. By combining existing data syntheses and model outputs, the permafrost region is likely a wetland methane source and small terrestrial ecosystem CO 2 sink with lower net CO 2 uptake toward higher latitudes, excluding wildfire emissions. For 2002–2014, the strongest CO2 sink was located in western Canada (median: -52 g C m -2 y -1 ) and smallest sinks in Alaska, Canadian tundra, and Siberian tundra (medians: -5 to -9 g C m -2 y -1 ). Eurasian regions had the largest median wetland methane fluxes (16–18 g CH4 m -2 y -1 ). Quantifying the regional scale carbon balance remains challenging because of high spatial and temporal variability and relatively low density of observations. More accurate permafrost region carbon fluxes require: (a) the development of better maps characterizing wetlands and dynamics of vegetation and disturbances, including abrupt permafrost thaw; (b) the establishment of new year-round CO 2 and methane flux sites in underrepresented areas; and (c) improved models that better represent important permafrost carbon cycle dynamics, including non-growing season emissions and disturbance effects.

54 ENVIRONMENTAL SCIENCES↗

Experimental Whole-Ecosystem Warming Enables Novel Estimation of Snow Cover and Depth Sensitivities to Temperature, and Quantification of the Snow-Albedo Feedback Effect

Climate change is reducing the amount, duration, and extent of snow across high-latitude ecosystems. But, in landscapes where persistent winter snow cover develops, experimental platforms to specifically investigate interactions between warming and changes in snowpack, and impacts on ecosystem processes, have been lacking. We leveraged a whole-ecosystem warming experiment in a boreal peatland forest to quantify how snow duration, depth, and fractional cover vary with warming of up to +9°C. We found that every snow-related quantity we examined declined precipitously as the amount of warming increased. The importance of deep, continuous snow cover for moderating shallow soil temperature is highlighted by an increase in soil temperature variance and the frequency of short-duration freeze-thaw cycles in the warmer plots. We used a paired-plot approach to estimate the magnitude of the snow-albedo feedback effect, and demonstrate that albedo-driven warming linked to reduced snow cover varies between December (+0.4°C increase in maximum air temperature) and March (+1.2°C increase) because of differences in insolation. Overall, results show that even modest future warming will have profound impacts on northern winters and cold-season ecosystem processes. Plot-level data from this warming experiment, and emergent relationships between warming and quantities related to snow cover and duration, could be of enormous value for testing and improving the representation of snow processes in simulation models, especially under future climate scenarios that are outside of the range of historically observed variability.

54 ENVIRONMENTAL SCIENCES↗

Investigation of the Sensitivity of Tropical Cyclogenesis to Aerosol Intervention

Abstract As risks from tropical cyclones (TCs) are fueled by climate change escalation, there is an urgent need for transformational solutions to complement traditional approaches. Seeding TCs using aerosols can be a promising method to reduce cyclone intensity, supported by theoretical understanding of the microphysical effects of aerosols on TC clouds. The ideal time to intervene effectively in TCs is likely during their initial stage, before TC wind speeds reach their peak. However, studies exploring potential aerosol effects on TC formation remain scarce. This study investigates how a TC embryo responds to the addition of aerosols of varying sizes using the Weather Research & Forecasting (WRF) model coupled with a spectral‐bin microphysics model. We found that aerosols of different sizes and concentrations distinctively affect the pre‐TC vortex's microstructure and dynamics. Fine and ultrafine aerosols enhance the latent heat of condensation, freezing, deposition, and riming, initially intensifying the vortex. However, this results in enhancement of the cold pool, thereby reducing inflow and surface fluxes, subsequently weakening the vortex. Coarse aerosols produce the opposite effect to that of fine and ultrafine aerosols. Coarse aerosols lead to a slower initial acceleration owing to enhanced warm rain. However, the resulting weaker cold pool is insufficient to effectively reduce the strength of the vortex at the later stage. This study provides critical insights into how aerosols of varying sizes and concentrations modulate the energy cascade and impact the evolution of a TC embryo, laying the groundwork for further research on TC risk management through aerosol intervention.

54 ENVIRONMENTAL SCIENCES↗

Global Simulations Suggest Biomass Burning Aerosol Emissions From Grassland Fires Could Be Important Ice Nucleating Particles

Ice nucleating particles (INP) capable of nucleating ice crystals via immersion freezing at temperatures above approximately −35°C may strongly influence cloud glaciation, with implications for global precipitation and climate feedback. In addition to mineral dust, soil dust, and marine organics, laboratory and field measurements suggest biomass burning aerosols (BBA) can act as immersion-mode INP between around −30°C and −15°C. However, the contribution of BBA to the global INP budget remains poorly understood due to poor knowledge of which fuels yield INPs, uncertainties in global coverage of those fuels, and unknown size distributions of the INPs in the BBA. Nonetheless, with some understanding of these uncertainties from sensitivity studies, the relative importance of ice nucleation activity of BBA compared to other INP sources can be quantified. In this work, we investigate the potential global importance of BBA as INP using a global aerosol-climate model, specifically the UK Met Office Unified Model (UM). We evaluate the model using field campaign data sets. We examine potential uncertainties in fuel types and particle sizes on BBA-based INP concentrations. Averaged over June–September between 15°S and 50°S, BBA is a more important INP than dust and marine INP about 30% of the time at altitudes with temperatures between −30°C and −20°C. Our simulations therefore suggest BBA INPs may be at least as important as mineral dust and marine INP over the atmospheric regions and seasons where grassland fires are frequent.

Gohil, Kanishk [Carnegie Mellon University, Pittsb↗

Variability in Ice Nucleating Particles Across Greater Houston Texas

The concentration and cloud-forming potential of a region's ice nucleating particle (INP) population have uncertain impacts on deep convective clouds. Specifically, ice nucleating particles (INPs) may affect various cloud properties related to the formation, lifetime, and precipitation of deep convective clouds. As part of the U.S. Department of Energy's TRacking Aerosol and Convection interaction ExpeRiment (TRACER) campaign, researchers from Texas A&M University deployed three Davis Rotating-drum Universal-size-cut Monitoring (DRUM) samplers throughout Greater Houston, Texas from June through September 2022. Ambient particles, collected at the surface with the DRUM samplers in four aerodynamic diameter size ranges (>3, 3–1.2, 1.2–0.34, and 0.34–0.15 μm), were analyzed in offline cold-stage ice nucleation experiments. The INP population in Greater Houston is complex, varying by site and day, but can be generalized by a weak to moderately efficient mode of INPs at −24°C and an efficient mode at −15°C. Analysis reveals that supermicron particles are largely responsible for ice nucleation warmer than −20°C across the region while submicron particles dominate at temperatures colder than −20°C. Additionally, significant spatial diversity in the INP population was observed, with differences in mean nucleation temperature between sites for nearly every size cut. Although INP concentrations were typically ∼0.08 L −1 at −20°C throughout the campaign, a notable region-wide increase in INP concentration for particles freezing at temperatures warmer than −20°C occurred from mid-August to mid-September. This comprehensive characterization of Greater Houston's INP population, including spatial, temporal, and particle size variations, can help constrain ice microphysics parameterizations in weather and climate models.

Thompson, Seth A. [Texas A & M Univ., College Stat↗

Passive Microwave Radiometry and Active Radar Sounding as Complementary Tools for Geophysical Investigations of Icy Ocean Worlds

Juno Microwave Radiometer (MWR) observations of Europa and Ganymede offer critical insights into the icy shells of these moons ahead of NASA's Europa Clipper and ESA's JUpiter ICy moons Explorer (JUICE) missions. Both missions are equipped with active radar sounders designed to address key unknowns such as ice shell thickness, thermal state, and composition. In this study, we explore how passive microwave radiometry and active radar sounding can constrain ice shell properties, focusing on Europa. Using modeled microwave brightness temperature observations at 0.6 and 1.2 GHz alongside simulated radar attenuation rate observations, we show that each instrument can independently produce robust ice shell thickness constraints under idealized conditions. We then relax these assumptions, quantifying how uncertainties from non-ideal properties—including convective layers, freezing-point depression, and chloride-doped ice—affect thickness estimates. Finally, we demonstrate how combining observations from these complementary techniques breaks degeneracies between ice shell properties, enabling more robust constraints than either method alone. This approach will maximize the science return of Europa Clipper and JUICE, advancing our understanding of the thermophysical structure and habitability of icy ocean worlds.

58 GEOSCIENCES↗

A bioinspired and chemically defined alternative to dimethyl sulfoxide for the cryopreservation of human hematopoietic stem cells

Abstract The cryopreservation of hematopoietic cells using dimethyl sulfoxide (DMSO) and serum is a common procedure used in transplantation. However, DMSO has clinical and biological side effects due to its toxicity, and serum introduces variation and safety risks. Inspired by natural antifreeze proteins, a novel class of ice-interactive cryoprotectants was developed. The corresponding DMSO-, protein-, and serum-free cryopreservation media candidates were screened through a series of biological assays using human cell lines, peripheral blood cells, and bone marrow cells. XT-Thrive-A and XT-Thrive-B were identified as lead candidates to rival cryopreservation with 10% DMSO in serum based on post-thaw cell survival and short-term proliferation assays. The effectiveness of the novel cryopreservation media in freezing hematopoietic stem cells from human whole bone marrow was assessed by extreme limiting dilution analysis in immunodeficient mice. Stem cell frequencies were measured 12 weeks after transplant based on bone marrow engraftment of erythroid, myeloid, B-lymphoid, and CD34 + progenitors measured by flow cytometry. The recovered numbers of cryopreserved stem cells were similar among XT-Thrive A, XT-Thrive B, and DMSO with serum groups. These findings show that cryoprotectants developed through biomimicry of natural antifreeze proteins offers a substitute for DMSO-based media for the cryopreservation of hematopoietic stem cells.

60 APPLIED LIFE SCIENCES↗

Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice

Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale <100 ns and widely separated from ice formation, which begins at times >1 μs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Momentum-independent magnetic excitation continuum in the honeycomb iridate H 3 LiIr 2 O 6

Understanding the interplay between the inherent disorder and the correlated fluctuating-spin ground state is a key element in the search for quantum spin liquids. H 3 LiIr 2 O 6 is considered to be a spin liquid that is proximate to the Kitaev-limit quantum spin liquid. Its ground state shows no magnetic order or spin freezing as expected for the spin liquid state. However, hydrogen zero-point motion and stacking faults are known to be present. The resulting bond disorder has been invoked to explain the existence of unexpected low-energy spin excitations, although data interpretation remains challenging. Here, we use resonant X-ray spectroscopies to map the collective excitations in H 3 LiIr 2 O 6 and characterize its magnetic state. In the low-temperature correlated state, we reveal a broad bandwidth of magnetic excitations. The central energy and the high-energy tail of the continuum are consistent with expectations for dominant ferromagnetic Kitaev interactions between dynamically fluctuating spins. Furthermore, the absence of a momentum dependence to these excitations are consistent with disorder-induced broken translational invariance. Our low-energy data and the energy and width of the crystal field excitations support an interpretation of H 3 LiIr 2 O 6 as a disordered topological spin liquid in close proximity to bond-disordered versions of the Kitaev quantum spin liquid.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗