Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stabilization method”

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 1,189 records · Page 66

Laser Speckle Pattern Stability in Hypersonic Regimes for Experimental Mechanics and Metrological Studies

Here, an empirical study assessing the suitability of laser speckle patterns for use in experimental mechanics research conducted at hypersonic conditions is presented. Speckle patterns generated via coherent optical scattering from a gold-coated graphite surface inserted into a Mach 10 air stream were collected and analyzed using Computer-Aided Speckle Interferometry using two methodologies: a typical spanwise approach involving increasing temporal span between image pairs, and a sequential approach that utilizes a combination of temporal re-referencing and spatial resampling. Results were characterized for spatial-temporal stability and variance. Our observations demonstrate that optical methods relying on stable and properly formed speckle patterns can be used in a hypersonic regime at various Reynolds numbers and common camera frame rates, but the technique is sensitive to the conditions of the airflow and the processing methodology chosen.

42 ENGINEERING↗

Subcell limiting strategies for discontinuous Galerkin spectral element methods

Here, we present a general family of subcell limiting strategies to construct robust high-order accurate nodal discontinuous Galerkin (DG) schemes. The main strategy is to construct compatible low order finite volume (FV) type discretizations that allow for convex blending with the high-order variant with the goal of guaranteeing additional properties, such as bounds on physical quantities and/or guaranteed entropy dissipation. For an implementation of this main strategy, four main ingredients are identified that may be combined in a flexible manner: (i) a nodal high-order DG method on Legendre–Gauss–Lobatto nodes, (ii) a compatible robust subcell FV scheme, (iii) a convex combination strategy for the two schemes, which can be element-wise or subcell-wise, and (iv) a strategy to compute the convex blending factors, which can be either based on heuristic troubled-cell indicators, or using ideas from flux-corrected transport methods. By carefully designing the metric terms of the subcell FV method, the resulting methods can be used on unstructured curvilinear meshes, are locally conservative, can handle strong shocks efficiently while directly guaranteeing physical bounds on quantities such as density, pressure or entropy. We further show that it is possible to choose the four ingredients to recover existing methods such as a provably entropy dissipative subcell shock-capturing approach or a sparse invariant domain preserving approach. We test the versatility of the presented strategies and mix and match the four ingredients to solve challenging simulation setups, such as the KPP problem (a hyperbolic conservation law with non-convex flux function), turbulent and hypersonic Euler simulations, and MHD problems featuring shocks and turbulence.

97 MATHEMATICS AND COMPUTING↗

A simple-structured benzo[1,2-b:4,5-b’]dithiophene-based molecule as a third component enables organic solar cells with over 18.2% efficiency

In this study, research has shown that ternary strategy is an effective tool for improving the performance of organic solar cells (OSCs) and simple-structured small molecules with low-cost and wide bandgap which broadens the scope of the choice in third components recently attracted the attention of researchers. Herein, a small molecule with benzo[1,2-b:4,5-b’]dithiophene (BDT) core named BR-C12 which mixed well with PM6 has been designed and synthesized. Incorporating BR-C12 as third component into PM6:Y6, BR-C12 exhibits a complementary absorption with PM6 and Y6, promoted exciton dissociation and charge collection, reduced bimolecular recombination, resulting in more ordered molecular packing. The PM6:Y6:BR-C12 (1:1.2:0.1, weight ratio) OSC shows an improved power conversion efficiency (PCE) of 17.02% in comparison to host binary devises (16.42%). In addition, by employing BR-C12 as third component into PM6:L8-BO, the ternary device based on PM6:L8-BO:BR-C12 also achieves significant enhancement in device photovoltaic performance with a PCE of 18.26%. In addition, BR-C12-added ternary devices exhibit improved stability than binary ones. This work presents a method for constructing high-performance ternary OSCs and enriches the abundance of materials for the third components.

14 SOLAR ENERGY↗

Phase transitions of BaTi{sub 0.9}Rh{sub 0.1}O{sub 3±δ} perovskite-type oxides under reducing environments

Highlights: • Solid solution formation BaTi{sub 0.9}Rh{sub 0.1}O{sub 3±δ} with a new wet chemical synthesis method. • Rhodium in the BaTiO{sub 3} perovskite stabilizes the hexagonal structure. • New Rh segregation mechanism for hexagonal BaTi{sub 0.9}Rh{sub 0.1}O{sub 3±δ} upon reduction. - Abstract: Perovskite-type oxides of composition BaTi{sub 0.9}Rh{sub 0.1}O{sub 3±δ} were prepared following a new chemical route that avoids the formation of hydroxyl species and precipitation, and allows the homogeneous distribution of Rh in the final mixed metal oxide. The high dispersion of Rh and the formation of the solid solution between Rh and the BaTiO{sub 3} perovskite is confirmed by means of X-ray diffraction (XRD) and extended X-ray absorption fine structure spectroscopy (EXAFS). The presence of Rh stabilized the hexagonal BaTi{sub 0.9}Rh{sub 0.1}O{sub 3±δ} phase, which decomposes into barium orthotitanate (BaTi{sub 2}O{sub 4}) and metallic Rh° in reducing environment. This phase transition starts already at 700 °C and is only partially completed at 900 °C suggesting that part of the Rh present in the perovskite lattice might not be easily reduced by hydrogen. These aspects and further open questions are discussed.

36 MATERIALS SCIENCE↗

Structural and functional impact of troponin C-mediated Ca 2+ sensitization on myofilament lattice spacing and cross-bridge mechanics in mouse cardiac muscle

Acto-myosin cross-bridge kinetics are important for beat-to-beat regulation of cardiac contractility; however, physiological and pathophysiological mechanisms for regulation of contractile kinetics are incompletely understood. Here we explored whether thin filament-mediated Ca 2+ sensitization influences cross-bridge kinetics in permeabilized, osmotically compressed cardiac muscle preparations. We used a murine model of hypertrophic cardiomyopathy (HCM) harboring a cardiac troponin C (cTnC) Ca 2+ -sensitizing mutation, Ala8Val in the regulatory N-domain. We also treated wild-type murine muscle with bepridil, a cTnC-targeting Ca 2+ sensitizer. Our findings suggest that both methods of increasing myofilament Ca 2+ sensitivity increase cross-bridge cycling rate measured by the rate of tension redevelopment (k TR ); force per cross-bridge was also enhanced as measured by sinusoidal stiffness and I 1,1 -1/I 1,0 ratio from X-ray diffraction. Computational modeling suggests that Ca 2+ sensitization through this cTnC mutation or bepridil accelerates k TR primarily by promoting faster cross-bridge detachment. To elucidate if myofilament structural rearrangements are associated with changes in k(TR), we used small angle X-ray diffraction to simultaneously measure myofilament lattice spacing and isometric force during steady-state Ca 2+ activations. Within in vivo lattice dimensions, lattice spacing and steady-state isometric force increased significantly at submaximal activation. We conclude that the cTnC N-domain controls force by modulating both the number and rate of cycling cross-bridges, and that the both methods of Ca 2+ sensitization may act through stabilization of cTnC's D-helix. Furthermore, we propose that the transient expansion of the myofilament lattice during Ca 2+ activation may be an additional factor that could increase the rate of cross-bridge cycling in cardiac muscle. These findings may have implications for the pathophysiology of HCM.

60 APPLIED LIFE SCIENCES↗

Phase Imaging: A Compressive Sensing Approach

Since Wolfgang Pauli posed the question in 1933, whether the probability densities |Ψ(r)|² (real-space image) and |Ψ(q)|² (reciprocal space image) uniquely determine the wave function Ψ(r) [1], the so called Pauli Problem sparked numerous methods in all fields of microscopy [2, 3]. Reconstructing the complete wave function Ψ(r) = a(r)e-iφ(r) with the amplitude a(r) and the phase φ(r) from the recorded intensity enables the possibility to directly study the electric and magnetic properties of the sample through the phase. In transmission electron microscopy (TEM), electron holography is by far the most established method for phase reconstruction [4]. Requiring a high stability of the microscope, next to the installation of a biprism in the TEM, holography cannot be applied to any microscope straightforwardly. Recently, a phase retrieval approach was proposed using conventional TEM electron diffractive imaging (EDI). Using the SAD aperture as reciprocal-space constraint, a localized sample structure can be reconstructed from its diffraction pattern and a real-space image using the hybrid input-output algorithm [5]. We present an alternative approach using compressive phase-retrieval [6]. Our approach does not require a real-space image. Instead, random complimentary pairs of checkerboard masks are cut into a 200 nm Pt foil covering a conventional TEM aperture (cf. Figure 1). Used as SAD aperture, subsequently diffraction patterns are recorded from the same sample area. Hereby every mask blocks different parts of gold particles on a carbon support (cf. Figure 2). The compressive sensing problem has the following formulation. First, we note that the complex-valued reciprocal-space wave-function is the Fourier transform of the (also complex-valued) real-space wave-function, Ψ(q) = F[Ψ(r)], and subsequently the diffraction pattern image is given by |Ψ(q)|2 = |F[Ψ(r)]|2. We want to find Ψ(r) given a few differently coded diffraction pattern measurements yn = |F[HnΨ(r)]|2, where the matrices Hn encode the mask structure of the aperture. This is a nonlinear inverse problem, but has been shown to be solvable even in the underdetermined case [6]. Since each diffraction pattern yn contains diffraction information from selected regions of the same sample, the differences in each pattern contain local phase information, which can be combined to form a full estimate of the real-space wave-function[7]. References: [1] W. Pauli in “Die allgemeinen Prinzipien der Wellenmechanik“, ed. H Geiger and W Scheel, (Julius Springer, Berlin). [2] A. Tonomura, Rev. Mod. Phys. 59 (1987), p. 639. [3] J. Miao et al, Nature 400 (1999), p. 342. [4] H. Lichte et al, Annu. Rev. Mater. Res. 37 (2007), p. 539. [5] J. Yamasaki et al, Appl. Phys. Lett. 101 (2012), 234105. [6] P Schniter and S Rangan. Signal Proc., IEEE Trans. on. 64(4), (2015), pp. 1043. [7] Supported by the Chemical Imaging, Signature Discovery, and Analytics in Motion initiatives at PNNL. PNNL is operated by Battelle Memorial Inst. for the US DOE; contract DE-AC05-76RL01830.

Schneider, Sebastian↗

Consistency-Enhanced Evolution for Variable Selection Can Identify Key Chemical Information from Spectroscopic Data

In the last few decades, spectroscopic techniques such as near-infrared (NIR) spectroscopy have gained wide applications in several industries, such as the pharmaceutical, agricultural, oil, and gas industries. As a result, various soft sensors have been developed to predict sample properties from spectroscopic readings. Because the spectroscopic readings at different wavelengths, especially at the adjacent wavelengths, are highly correlated, it has been shown that variable selection could significantly improve a soft sensor’s prediction performance while reducing the model complexity. To improve the prediction performance, most variable selection methods focus on identifying the variables (i.e., wavelengths or wavelength segments) that are strongly correlated with the dependent variable. Although many successful applications have been reported, these variable selection methods do have their limitations. Specifically, the selected wavelengths sometimes show little connection to the chemical bounds or functional groups presenting in the sample. In addition, the selected variables can be quite sensitive to the choice of the training samples. In this work, we address these limitations from a different perspective: if a variable selection algorithm can identify the truly relevant input variables, it should consistently identify the same subset of variables regardless of the choice of the training samples. Therefore, we propose a variable selection method that aims to improve the consistency of variable selection resulting from different training samples. Furthermore, the new algorithm is termed consistency-enhanced evolution for variable selection (CEEVS). To demonstrate the performance and robustness of CEEVS, we compare the proposed method with three representative variable selection methods using five published NIR data sets. These case studies clearly demonstrate that by improving the variable selection consistency, we can not only achieve improved prediction performance, but also identify key chemical information from spectroscopic data.

42 ENGINEERING↗

Inelastic Neutron Scattering from Thin Film Biaxially Oriented Polyethylene Terephthalate

Recent interest in emerging processes for polymer manufacturing and bio-based chemistries for direct chemical recycling/upcycling has motivated new research focused on a deeper understanding of atomic-scale polymer properties and how they influence macroscopic phenomena. Uncovering the fundamental properties of polymers that give rise to macroscopic behavior could enable new pathways for improved recyclability or utilization of alternative “greener” polymer analogues. In this study, the neutron vibrational spectrum was measured for a film of biaxially oriented polyethylene terephthalate (BoPET) using inelastic neutron scattering (INS), to investigate the relationship between the structure and dynamics of a widely used polymer. Compared to conventional spectroscopic techniques, the use of INS is advantageous for polymeric materials due to the absence of selection rules (i.e., all transitions are allowed), broad-band energy range, and considerable sensitivity to hydrogen modes. In order to distinguish the vibrational modes caused by trans and gauche rotational isomerism, the normal modes of vibration were calculated from a density functional theory-optimized structure of crystalline PET (cPET), representative of the all-trans state, and compared with INS from “highly crystalline” PET powder. Although in- and out-of-plane wagging of hydrogens on the ring structure exhibit significant contribution to both BoPET and cPET spectra, the wagging, rocking, and twisting modes of hydrogen on the ethylene glycol group are, in most cases, conformation-specific. These results were further rationalized by investigating the role of hyperconjugation in stabilizing both conformations using the natural bond order method. Here, through comparison of experimental and calculated INS results, this work provides the fundamental basis for discovering the role of structure and dynamics in shaping the macroscopic properties of PET and polymer analogues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved 140 Nd Production for the 140 Nd/ 140 Pr In Vivo Generator through Target Recycling and Radiochemical Optimization

Theranostic strategies that utilize f-block therapeutic radionuclides, including 161 Tb, 177 Lu, 225 Ac, and 227 Th, suffer from a shortage of positron emission tomography (PET) imaging counterparts in the same chemical space and often rely on 68 Ga as a surrogate. The 140 Nd/ 140 Pr in vivo PET generator, which belongs to the f-block, may address this issue and can be produced via the 141 Pr(p,2n) 140 Nd production route by using medium-energy cyclotrons. However, impurities in the target material, including stable Nd, and the inherent difficulty of adjacent lanthanide separations limit the achievable radionuclidic and chemical purity of 140 Nd. In this work, we address these challenges through the purification and recycling of praseodymium target material and optimization of Nd/Pr separation. The resulting purified 140Nd was evaluated using DOTA and Macropa chelators via radiolabeling and in vitro stability studies. A target material purification and recycling method was developed for the monoisotopic 141 Pr starting material to remove stable Nd impurities, yielding 90.3 ± 4.7% (n = 3) recovery. The purified 141 Pr was isolated as Pr 6 O 11 and irradiated with 24 MeV protons (20.07 MeV at the target surface) at 20 μA for 4 h, which produced 1417.0 ± 83.4 MBq (38.3 ± 2.2 mCi) of 140 Nd at the end of bombardment (EOB). The produced 140 Nd was purified through an optimized DGA normal method to recover 71.6 ± 6.3% pure 140 Nd. The amount of stable Nd reduced progressively in each target purification cycle from >340 ppm without purification to <250 ppb after three cycles, while other measured metallic impurities were below 30 ppb. This improvement in target purity was reflected in the direct increase of apparent molar activity (AMA), when purified 140 Nd was evaluated with DOTA and Macropa chelators. AMA of [ 140 Nd]Nd-DOTA and [ 140 Nd]Nd-Macropa increased from 70.3 MBq/μmol (1.9 mCi/μmol) and 74 MBq/μmol (2.0 mCi/μmol) to 8025.3 MBq/μmol (216.9 mCi/μmol) and 8473.0 MBq/μmol (229.0 mCi/μmol), respectively, after the third target purification cycle. Further evaluation of chelator-labeled 140 Nd showed that [ 140 Nd]Nd-DOTA was stable in phosphate-buffered saline (PBS), saline, human serum, and mouse serum, whereas [140Nd]Nd-Macropa was stable in all except human serum. This work established a practical methodological advance for the production of 140 Nd/ 140 Pr in vivo PET generators, combining optimized target recycling and radiochemical separation to enable scaled-up and high-molar activity 140 Nd suitable for preclinical imaging. These advances support broader development of 140 Nd/ 140 Pr as a robust PET analogue, especially for f-block therapeutics.

Irradiation↗

Synthetic fossilization of soft biological tissues and their shape-preserving transformation into silica or electron-conductive replicas

Structural preservation of complex biological systems from the subcellular to whole organism level in robust forms, enabling dissection and imaging while preserving 3D context, represents an enduring grand challenge in biology. Here we report a simple immersion method for structurally preserving intact organisms via conformal stabilization within silica. This self-limiting process, which we refer to as silica bioreplication, occurs by condensation of water-soluble silicic acid proximally to biomolecular interfaces throughout the organism. Conformal nanoscopic silicification of all biomolecular features imparts structural rigidity enabling the preservation of shape and nano-to-macroscale dimensional features upon drying to form a biocomposite and further high temperature oxidative calcination to form silica replicas or reductive pyrolysis to form electrically conductive carbon replicas of complete organisms. Ultimately, the simplicity and generalizability of this approach should facilitate efforts in biological preservation and analysis and could enable the development of new classes of biomimetic composite materials.

59 BASIC BIOLOGICAL SCIENCES↗

Nanocluster and single-atom catalysts for thermocatalytic conversion of CO and CO 2

In this Minireview, we discuss recent advances in understanding nanocluster and single-atom catalysts for CO and CO 2 emissions control applications. Through analyzing thermocatalytic CO oxidation and CO 2 reduction, two fundamentally and industrially important reactions, we compare representative nanocluster and single-atom catalytic systems from perspectives of intrinsic chemistry and reaction engineering. Generally, nanoclusters and single atoms display different catalytic performances (i.e., activity, selectivity, and stability) for these reactions depending on the synthesis methods, support materials, and reaction conditions. Key observations of activity and selectivity tradeoffs between nanoclusters and single-atom catalysts are highlighted. Here, the dynamic structural responses of these catalytic species under CO oxidation or CO 2 reduction reaction conditions are also discussed. Synthetic control and detailed experimental and computational characterization of single-atom and nanocluster catalysts for CO and CO 2 conversion have led to exciting progress over the past decade. Still, more efforts are needed to understand and develop catalysts that meet the environmental, energy, and technical requirements to power a sustainable global economy.

03 NATURAL GAS↗

Condition‐based exciter model complexity reduction for improved transient stability simulation

A reduced model approach to decrease computational complexity in power system transient stability simulation is presented. This study investigates conditions that make fast modes active or inactive using practical power system examples. Modes in the original system where fast dynamics do not appear can be neglected, allowing simulation steps to be increased without numerical stability issues. During a transient simulation, the proposed method switches dynamically between the original exciter model and the reduced one, depending on the switching criterion presented. A high level of accuracy is achieved and computational efficiency is increased. Case studies with the Western Electricity Coordinating Council system are provided to validate the performance of the proposed method.

Kim, Soobae↗

Prediction of the phase equilibria of methane hydrates using the direct phase coexistence methodology

The direct phase coexistence method is used for the determination of the three-phase coexistence line of sI methane hydrates. Molecular dynamics (MD) simulations are carried out in the isothermal–isobaric ensemble in order to determine the coexistence temperature (T{sub 3}) at four different pressures, namely, 40, 100, 400, and 600 bar. Methane bubble formation that results in supersaturation of water with methane is generally avoided. The observed stochasticity of the hydrate growth and dissociation processes, which can be misleading in the determination of T{sub 3}, is treated with long simulations in the range of 1000–4000 ns and a relatively large number of independent runs. Statistical averaging of 25 runs per pressure results in T{sub 3} predictions that are found to deviate systematically by approximately 3.5 K from the experimental values. This is in good agreement with the deviation of 3.15 K between the prediction of TIP4P/Ice water force field used and the experimental melting temperature of ice Ih. The current results offer the most consistent and accurate predictions from MD simulation for the determination of T{sub 3} of methane hydrates. Methane solubility values are also calculated at the predicted equilibrium conditions and are found in good agreement with continuum-scale models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design of semiconductor surface pits for fabrication of regular arrays of quantum dots and nanorings

We report a systematic computational study on the formation of complex nanostructures consisting of quantum dots and nanorings on surfaces of coherently strained thin films grown epitaxially on pit-patterned substrates. The analysis is based on self-consistent dynamical simulations according to a film surface evolution model that has been validated experimentally by comparison of its predictions with experimental observations on Ge films grown on Si pit-patterned substrates and is supported by linear stability theory that explains the film surface nanopattern formation as the outcome of a Stranski-Krastanow instability. Emphasis is placed on the design of conical pit patterns and the effects on the resulting film surface nanopattern of varying geometrical design parameters including film thickness, pit-pattern period, pit depth, pit opening diameter, and pit wall inclination. We demonstrate that varying the pit opening diameter and the pit wall slope leads to formation of complex nanostructures inside the pits of a regular pit pattern on the film surface, which include quantum dots, as well as single nanorings and multiple concentric nanorings that may or may not surround a central quantum dot inside each pit. Our simulation predictions demonstrate that the ordered nanostructure patterns forming on the film surface can be controlled precisely by tuning the geometrical parameters of the pits on the pit-patterned substrate. Our findings have important implications for designing optimal semiconductor surface patterns toward enabling future nanofabrication technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fabrication of pristine 2D heterostructures for scanning probe microscopy

Material-by-design has been a long-standing aspiration that has recently become a reality. Such designer materials have been repeatedly demonstrated using the top-down approach of mechanical exfoliation and stacking, leading to a variety of artificial 2D heterostructures with new properties that are otherwise unattainable. Consequently, tremendous research frontiers in physics, chemistry, engineering, and life science have been created. While thousands of layered crystals exist in nature, only a few dozen of them with manageable chemical-stability have been made into heterostructures using this method. Moreover, experimental investigations of materials that have received limited exploration in the 2D realm, such as cuprates, halides, and perovskites, along with their heterostructures, have been fundamentally hindered by their rapid chemical degradation. Another critical challenge imposed by exfoliating and stacking 2D layers in ambient environment is the absorption of itinerant gas molecules that further contaminate sensitive 2D interfaces in the heterostructures. Such contamination and compromised material properties significantly hinder surface-sensitive local probes—scanning probe microscopy (SPM)—that often require nanometer to atomic scale surface cleanliness. In this article, we aim to provide a technical review of recent development toward 2D materials and heterostructure fabrication in more controlled environments that are suitable for SPM characterizations. These include the development of more efficient mechanical exfoliation and dry-transfer techniques, as well as the incorporation of 2D material exfoliation and transfer in inert gas, low vacuum, and, eventually, ultra-high vacuum environments. Finally, we provide an outlook on the remaining challenges and opportunities in ultra-clean 2D material fabrication techniques.

Materials Science↗

A high-order spectral deferred correction strategy for low Mach number flow with complex chemistry

This article presents a fourth-order finite-volume algorithm in space and time for low Mach number reacting flow with detailed kinetics and transport. Our temporal integration scheme is based on a Multi-Implicit Spectral Deferred Correction (MISDC) strategy that iteratively couples advection, diffusion, and reactions evolving subject to a constraint. Our new approach overcomes a stability limitation of our previous second-order method encountered when trying to incorporate higher-order polynomial representations of the solution in time to increase accuracy. We have developed a new iterative scheme that naturally fits within our MISDC framework and allows us to conserve mass and energy while simultaneously satisfying the equation of state. We analyse the conditions for which the iterative schemes are guaranteed to converge to the fixed point solution. Finally, we present numerical examples illustrating the performance of the new method on premixed hydrogen, methane, and dimethyl ether flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of broken KAM surfaces for particle orbits in toroidal confinement systems

Here, the destruction of Kolmogorov–Arnold–Moser surfaces in a Hamiltonian system is an important topic in nonlinear dynamics, and in particular in the theory of particle orbits in toroidal magnetic confinement systems. Analytic models for transport due to mode-particle resonances are not sufficiently correct to give the effect of these resonances on transport. In this paper we compare three different methods for the detection of the loss of stability of orbits in the dynamics of charged particles in a toroidal magnetic confinement device in the presence of time dependent magnetic perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient Simulation of Logical Magic State Preparation Protocols

Developing space- and time-efficient logical magic state preparation (MSP) protocols will likely be an essential step toward building a large-scale fault-tolerant quantum computer. Motivated by this need, we introduce a scalable method for simulating logical MSP protocols under the standard circuit-level noise model. When applied to protocols based on code-switching, magic state cultivation, and magic state distillation, our method yields a complexity polynomial in (i) the number of qubits and (ii) the nonstabilizerness, e.g., stabilizer rank or Pauli rank, of the target encoded magic state. The efficiency of our simulation method is rooted in a curious fact: every circuit-level Pauli error in these protocols propagates to a Clifford error at the end. This property is satisfied by a large family of protocols, including those that repeatedly measure a transversal Clifford that squares to a Pauli. We provide a proof-of-principle numerical simulation that prepares a magic state using such logical Clifford measurements. Our work enables practical simulation of logical MSP protocols without resorting to approximations or resource-intensive state-vector simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗