Search NASASearch

SEARCH · Search NASA

Results for “inference methods”

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.

559 records · Page 6

Evaluation of Low-Energy Hydrogen Separation Method Using Metal-Organic Frameworks (MOFs) for Closed-Loop ECLSS Air Revitalization (CLEAR)

The State-of-the-Art (SOA) air revitalization architecture onboard the International Space Station (ISS) recovered approximately 50% of the oxygen (O 2 ) from metabolic carbon dioxide (CO 2 ) via the Sabatier process from 2011 to 2017. O 2 recovery is currently constrained by the limited availability of reactant hydrogen (H 2 ) preventing complete conversion of CO 2 to H 2 O. Increasing O 2 recovery within Closed-Loop ECLSS is essential to reducing resupply mass for long-duration manned missions; specifically focusing on water (H 2 O) which supplies H 2 for Sabatier via water electrolysis. Past ground test endeavors at Marshall Space Flight Center (MSFC) have attempted to recover H 2 from Sabatier-produced CH 4 using technologies such as carbon vapor deposition (CVD) and plasma pyrolysis. The byproducts of these technologies can act as a catalyst poison or reactor deadload to the Sabatier reactor. Hydrogen separation techniques must be utilized to maintain the Sabatier catalyst during gas recycling and must be scalable, non-energy intensive, and safe to operate in a habitation setting. Research indicated that metal-organic frameworks (MOFs) could meet these criteria and were tested for their capability to capture the various carbon-based gaseous products of CVD and plasma pyrolysis such as acetylene (C 2 H 2 ), ethylene (C 2 H 4 ), ethane (C 2 H 6 ), and carbon monoxide (CO) which would purify the hydrogen gas stream passing through the MOF. A sub-scale adsorption column was developed by Marshall Space Flight Center to test three MOF candidates against a synthetic gas mixture comprised of process-relevant carbonous gases and hydrogen to evaluate the separation capability of the MOFs. The results of the hydrogen separation capability, isothermal desorption capability, and demonstrated cyclic reuse of the MOF are presented in this paper.

Kagen Crawford

Performance Optimization for Lunar Extravehicular Activity Readiness (POLAR) Study: Methods Paper

To better understand which aspects of physical fitness may be most related to performance during Lunar surface operations and thereby help to inform the current NASA fitness standards, much can be learned from fields encompassing the “tactical athlete.” Other physically demanding professions such as law enforcement, military, or rescue professionals often require candidates to meet occupationally-relevant fitness standards. The determination of such standards is a multistep process, including both objective and subjective measures, to determine tasks essential to occupational performance and identify the minimal fitness profile needed to meet physical demands of the job. Notably, fitness is only one component which may contribute to the demands of astronaut selection, flight assignment, and occupational performance. Utilizing a framework to systematically determine which domains of fitness most contribute to relevant job tasks can aid in the refinement of current NASA-STD-3001 fitness standards. Therefore, NASA’s Exercise Physiology & Countermeasures Laboratory conducted the Performance Optimization for Lunar Extravehicular Activity Readiness (POLAR) study to identify and examine a comprehensive list of fitness tests (including NASA-STD-3001 assessments: 1-Repetition Maximum [1-RM] bench press and deadlift) and determine preliminary relationships between identified fitness parameters and novel Artemis-relevant tasks to help inform future investigations for the continued development of aerobic and muscular fitness standards for surface EVAs. This was accomplished through 1) a review of the literature relating fitness assessments to simulated or real EVA performance to identify fitness tests that are most correlated with simulated EVA task performance; 2) a task analysis following a modified framework for physical employment standards development to down select mission critical tasks; and 3) development and pilot testing of a novel, portable Artemis-relevant EVA task circuit to relate to a battery of fitness assessments. The current report describes the methodology used to complete the task analysis, EVA task circuit development, and the pilot study.

Nicole C Strock

Variational Quantum Circuits to Prepare Low Energy Symmetry States

We explore how to build quantum circuits that compute the lowest energy state corresponding to a given Hamiltonian within a symmetry subspace by explicitly encoding it into the circuit. We create an explicit unitary and a variationally trained unitary that maps any vector output by ansatz A(α → ) from a defined subspace to a vector in the symmetry space. The parameters are trained varitionally to minimize the energy, thus keeping the output within the labelled symmetry value. The method was tested for a spin XXZ Hamiltonian using rotation and reflection symmetry and H 2 Hamiltonian within S z = 0 subspace using S 2 symmetry. We have found the variationally trained unitary gives good results with very low depth circuits and can thus be used to prepare symmetry states within near term quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Integration of ultra-low coverage whole-genome sequences for reconstructing the evolutionary history of Galapagos giant tortoises

Genomic data from contemporary and historical samples often need to be coupled for evolutionary reconstructions of multitaxon complexes. However, the genetic data recovered from historical samples may result only in ultra-low coverage whole-genome sequences (ulcWGS; <0.15× depth), leading to inaccurate evolutionary inferences given a preponderance of missing data. Using the Galapagos giant tortoise radiation as a study system (Chelonoidis spp., composed of 13 extant and four extinct lineages), we assembled a novel methodological pipeline that removes potential noise introduced by the missing data and enhances the evolutionary signal from ulcWGS samples. We leveraged existing tools for phylogenomic placement (EPA-ng), population genomic structure (smartsnp) and admixture (Admixfrog, NGSadmix) to demonstrate that the evolutionary history of samples can be uncovered with sequencing depths as low as 0.008–0.139×. Importantly, these approaches do not use genotype imputation of the ulcWGS samples, which would require extensive reference datasets. Our application to two cases of extinct lineages of Galapagos giant tortoises, with and without references from the same lineage, demonstrates the general value of the approach. We confirm where the extinct lineages from San Cristóbal and Santa Fe islands fit into the Galapagos giant tortoise radiation, and that these lineages were evolutionarily distinct entities.

ancient DNA

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven

230Th/234U radiochronometry for uranium materials by alpha spectrometry for nuclear forensics analysis

We have developed an alpha spectrometry method for 230Th/234U radiochronometry to determine the model separation age for uranium materials. This method offers an alternative to the more commonly used, but significantly more resource intensive multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) while providing similar accuracy and precision. The development included radiochemical separation of uranium and thorium, calibration of the 232U tracer and the alpha spectrometry measurement. The method was validated by analyzing a certified reference material, CRM 125-A, which is a uranium oxide pellet assay, isotopic, and radiochronometric standard. The results were compared to those obtained by our routine MC-ICP-MS 230Th/234U radiochronometry analysis method to assess the performance parameters for alpha spectrometry compared to the state-of-the-art technique.

Macsik, Zsuzsanna [Los Alamos National Laboratory

Time evolution of a pumped molecular magnet—A time-resolved inelastic neutron scattering study

Introducing an experimental technique of time-resolved inelastic neutron scattering (TRINS), we explore the time-dependent effects of resonant pulsed microwaves on the molecular magnet Cr 8 F 8 Piv 16 . The octagonal rings of magnetic Cr 3+ atoms with antiferromagnetic interactions form a singlet ground state with a weakly split triplet of excitations at 0.8 meV. A 4.6 tesla field was applied to tune the splitting between two members of the triplet excited level |1$\rangle$ ↔ |2$\rangle$ to resonance with 105 GHz (0.434 meV) microwaves. The time-dependent occupations of the ground state |0$\rangle$, lower lying levels |1$\rangle$ and |2$\rangle$, and higher energy states |λ ≥ 3$\rangle$ were extracted during and after 20 s long microwave pulses incident along the (101) direction of a Cr 8 F 8 Piv 16 crystal held at 1.9 K. At significantly elevated spin temperatures, we found underpopulation relative to thermal equilibrium of |2$\rangle$ and spin-lattice thermalization time scales ranging from 1.6(2) s to 5.7(2) s depending on the power level. This contrasts with the relaxation time τ 1 (T → 0) = 27(5) 𝛍 s inferred for |2$\rangle$ from in situ Electron Spin Resonance measurements. By probing a broad range of excited states during intense microwave pumping, TRINS thus provides a first view of long lived excited states in a molecular antiferromagnet.

Science & Technology - Other Topics

Kondo effect in ferromagnetic quantum critical CeRh 6 ⁢Ge 4

The mechanism of a pressure-induced quantum critical point in the heavy fermion ferromagnet CeRh 6 ⁢Ge 4 has attracted interest, as ferromagnetic quantum criticality in a clean itinerant Ce compound is typically avoided. The localized versus itinerant character of the 4⁢𝑓 electrons is a key aspect for understanding this behavior. We investigated the electronic structure of the 4⁢𝑓 shell in CeRh 6 ⁢Ge 4 using core-level photoelectron and x-ray absorption spectroscopy, demonstrating the hybridization of Ce 4⁢𝑓 with the conduction electrons. Linearly polarized x-ray absorption reveals a temperature-dependent linear dichroism consistent with the crystal-electric-field sequence as inferred from the static susceptibility. This dichroism cannot be described by an ionic full-multiplet model alone, but is reproduced by including the Kondo effect within a single-impurity Anderson model in the noncrossing approximation. The Kondo effect mixes higher-lying crystal-field states into a resulting multiorbital ground state with 4⁢𝑓 occupancy, 𝑛 𝑓 ∼ 0.9. Deviations at low temperatures between the measured linear dichroism and calculated dichroism suggest an orbital-dependent Kondo effect. A scenario in which there is a multiorbital ground state and orbital-dependent Kondo hybridization should be a starting point for a model of pressure-induced criticality in CeRh 6 ⁢Ge 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Electrochemical Deconstruction of Ortho-Phthalate Plasticizers and Recovery of Plasticizing Moieties

Electrochemical methods for the direct reduction of alkyl esters have been understudied but provide a potential advantage for addressing specific waste streams given that such strategies often require only a minimal chemical profile. One such waste stream that could benefit from electrochemical remediation is poly(vinyl chloride) (PVC) plastics. PVC recycling often faces challenges due to the complexity of such waste. Solvent-based recycling methods pose several advantages for addressing post-use PVC; however, such methods are complicated by the high amounts of toxic legacy plasticizers (ortho-phthalates) present in PVC. A potential solution to addressing such phthalates is to address them electrochemically, coupled with recovery of the resulting valuable products. Presented herein is an optimized electrochemical method for ester activation that is leveraged to separate and recover the aromatic and alkyl components of ortho-phthalate plasticizers. The resulting aliphatic alcohols may be reused to prepare other non-toxic plasticizers. This electro-degradation method is demonstrated on six phthalates that have been identified to pose health concerns in addition to plasticizers recovered from commercial samples of PVC.

di-(2-ethylhexyl) phthalate

Up-and-Coming IMCs

While the good oxidation and environmental resistance, high melting points, and comparatively low densities of such ordered intermetallics as Ti 3 Al, NiAl, FeAl, and NbAl 3 render them good candidates for advanced aerospace structures, their poor toughness at low temperatures and low strength at elevated temperatures have prompted the development of fiber-reinforced intermetallic-matrix composites (IMCs) with more balanced characteristics. Fabrication methods for continuous-fiber IMCs under development include the P/M 'powder cloth' method, the foil/fiber method, and thermal spraying. The ultimate success of IMCs depends on fibers truly compatible with the matrix materials.

Randy Bowman

Computing the Critical Temperature of the Affine-Transformed $D=3$ Ising Model Using Masked Autoregressive Flow

The simple Ising model provides a rich environment to build and study lattice field theories. As part of an ongoing project to construct a conformal field theory (CFT) on an arbitrarily curved manifold, in this work we develop methods to measure the critical temperature $β_c$ of the affine-transformed Ising model on the face-centered cubic (FCC) lattice. The main challenge in this endeavor is finding a computationally efficient and accurate method of interpolating and extrapolating Monte Carlo observables with respect to coupling coefficients and temperature. Herein, we compare two such methods. A traditional statistical approach uses the multiple histogram (MH) method, while a newer machine learning approach uses a masked autoregressive flow (MAF) to estimate the underlying probability density function of a set of observables. While the MH method is specifically designed to interpolate and extrapolate Monte Carlo observables, we find that MAF is a viable alternative for measuring $β_c$ with a computational cost that scales more favorably. Furthermore, we comment on additional advantages of MAF relevant to our work, such as extrapolating in system volume.

Svenson, Kai [Texas U.]

Cr plasma-material-interaction in PISCES-RF: D thermal release, retention, and erosion

Pure chromium (Cr) targets were exposed to high-flux deuterium (D) plasmas in the Pisces-RF linear plasma device, and measurements of D retention, release, and erosion were subsequently performed. Post-exposure D retention was quantified using temperature-programmed desorption on Cr targets irradiated by 50 eV ions over a broad range of exposure temperatures (423–873 K) and ion fluences (3 × 10 24 – 3 × 10 26 m −2 ). The retained D inventory was observed to decrease rapidly with increased exposure temperature, from approximately ∼7 × 10 20 m −2 at ∼ 420 K, to then saturate near ∼ 10 20 m −2 for exposure temperatures above ∼550 K. Separately, at fixed exposure temperature (∼450 K), D retention was found to have only a weak dependence on increasing ion fluence. Lastly, the erosion of Cr in D plasma was investigated for ion impact energies in the range 40 ≤ E i ≤ 250 eV. Erosion was inferred using optical emission spectroscopy (OES) from the ratio of emission lines (Cr I (425.4 nm)/D I (656.1 nm)) measured close to the target. Conversion of the OES yield data to net erosion yield was made with singular ion energy target mass-loss measurements. These net erosion yield data were then further corrected to obtain gross erosion yield by accounting for a re-deposition factor, computed using a simple model. The gross erosion yield is found to be 2–4 times lower than predicted by SDTrimSP, consistent with that typically observed for light-ion sputtering under high-flux plasma conditions.

Chromium

Rapid Assessment of Sulfate Resistance in Mortar and Concrete

Extensive research has been conducted on the sulfate attack of concrete structures; however, the need to adopt the use of more sustainable materials is driving a need for a quicker test method to assess sulfate resistance. This work presents accelerated methods that can reduce the time required for assessing the sulfate resistance of mixtures by 70%. Class F fly ash has historically been used in concrete mixtures to improve sulfate resistance. However, environmental considerations and the evolving energy industry have decreased its availability, requiring the identification of economically viable and environmentally friendly alternatives to fly ash. Another challenge in addressing sulfate attack durability issues in concrete is that the standard sulfate attack test (ASTM C1012) is time-consuming and designed for only standard mortars (not concrete mixtures). To expedite the testing process, accelerated testing methods for both mortar and concrete mixtures were adopted from previous work to further the development of the accelerated tests and to assess the feasibility of testing the sulfate resistance of mortar and concrete mixtures rapidly. This study also established criteria for interpreting sulfate resistance for each of the test methods used in this work. A total of 14 mortar mixtures and four concrete mixtures using two types of Portland cement (Type I and Type I/II) and various supplementary cementitious materials (SCMs) were evaluated in this study. The accelerated testing methods significantly reduced the evaluation time from 12 months to 21 days for mortar mixtures and from 6 months to 56 days for concrete mixtures. The proposed interpretation method for mortar accelerated test results showed acceptable consistency with the ACI 318-19 interpretations for ASTM C1012 results. The interpretation methods proposed for the two concrete sulfate attack tests demonstrated excellent consistency with the ASTM C1012 results from mortar mixtures with the same cementitious materials combinations. Metakaolin was shown to improve sulfate resistance for both mortar and concrete mixtures, while silica fume and natural pozzolan had a limited impact. Using 15% metakaolin in mortar or concrete mixtures with Type I/II cement provided the best sulfate resistance.

Chemistry

Elastic-plastic models for multi-site damage

This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.

Actis, Ricardo L.

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING