Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Configuration interaction”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effects of Coronal Magnetic Field Configuration on Particle Acceleration and Release during the Ground Level Enhancement Events in Solar Cycle 24

Ground level enhancements (GLEs) are extreme solar energetic particle (SEP) events that are of particular importance in space weather. In solar cycle 24, two GLEs were recorded on 2012 May 17 (GLE 71) and 2017 September 10 (GLE 72), respectively, using a range of advanced modern instruments. Here we conduct a comparative analysis of the two events by focusing on the effects of large-scale magnetic field configuration near active regions on particle acceleration and release. Although the active regions are both located near the western limb, temporal variations of SEP intensities and energy spectra measured in situ display different behaviors at early stages. By combining a potential field model, we find the coronal mass ejection (CME) in GLE 71 originated below the streamer belt, while in GLE 72 it originated near the edge of the streamer belt. We reconstruct the CME shock fronts with an ellipsoid model based on nearly simultaneous coronagraph images from multiple viewpoints and further derive the 3D shock geometry at the GLE onset. The highest-energy particles are primarily accelerated in the shock–streamer interaction regions, i.e., likely at the nose of the shock in GLE 71 and the eastern flank in GLE 72, due to quasi-perpendicular shock geometry and confinement of closed fields. Subsequently, they are released to the field lines connecting to near-Earth spacecraft when the shocks move through the streamer cusp region. This suggests that magnetic structures in the corona, especially shock–streamer interactions, may have played an important role in the acceleration and release of the highest-energy particles in the two events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2025

Optimization of spray breakup model parameters for predicting fuel spray and film characteristics in gasoline direct injection engines

This study investigated the behavior of gasoline direct injection (GDI) sprays using computational fluid dynamics (CFD). The authors developed an approach to identify optimal spray breakup model parameters by evaluating an error function across numerous simulations, with the goal of minimizing discrepancies from experimental data. Using the optimal setup, the simulated spray matched well with projected liquid volume distributions, liquid penetration, and spray width measured in a constant-pressure continuous-flow chamber. To further validate the approach, the same setup was tested across various fuels, injectors, and operating conditions. Subsequently, the optimal setup, along with a recently developed spray-wall interaction model, were applied to a direct-injection spark-ignited engine under late-injection conditions to predict and evaluate fuel film formation and evolution at varying engine coolant temperatures. Here, with the centrally mounted injector directing the spray toward the piston, simulations indicated that the spray tends to impinge on the piston surface. The proposed simulation framework also accurately captured the aggregate film area on the piston surface, aligning with previously published experimental results. Moreover, simulations showed that increasing the coolant temperature from cold start conditions (333 K) to warm conditions (363 K) reduced the fuel mass deposited on the piston by roughly 50%. Furthermore, for the spray-guided engine configuration studied in this work, the CFD model predicted minimal film deposition on the spark plug electrodes regardless of the coolant temperatures due to a relatively weak in-cylinder flow during the compression phase.

Computational fluid dynamics (CFD)↗
DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗
DOE OSTI2025

Mitigating urban climate-energy feedback with citywide building-integrated photovoltaics implementation

Urban heat islands (UHIs) intensify cooling demands, whose waste heat further exacerbates UHI, creating a self-reinforcing feedback loop. Building-integrated photovoltaics (BIPV) can break this cycle by generating electricity and altering urban climate-energy interactions, yet its net impact remains unresolved. Using a cross-scale modeling framework for 118,521 buildings in Hong Kong, we demonstrate that BIPV glazing reduces cooling demand through improved insulation and lower solar heat gain, while inducing diurnal thermal asymmetry of daytime warming (up to +1.9°C) and nighttime cooling (up to −2.5°C). This microclimate regulation indirectly contributes to a 0.6% reduction in citywide building energy use. Cumulative savings from electricity generation and reduced energy use reach 4.7% (2,060.2 GWh) and 10.1% (4,401.3 GWh) under low- and high-coverage deployment, respectively, with building-level savings of −20.1% ± 11.7% (mean ± SD) at high coverage. We establish BIPV as a configurable climate-energy regulator that enables net-zero, heat-resilient planning and quantifies deployment trade-offs in cooling-dominated cities worldwide.

Building-integrated photovoltaics (BIPV)↗
DOE OSTI2021

Valence shell photoelectron angular distributions and vibrationally resolved spectra of imidazole: A combined experimental–theoretical study

Linearly polarized synchrotron radiation has been used to record polarization dependent valence shell photoelectron spectra of imidazole in the photon energy range 21–100 eV. These have allowed the photoelectron angular distributions, as characterized by the anisotropy parameter β, and the electronic state intensity branching ratios to be determined. Complementing these experimental data, theoretical photoionization partial cross sections and β-parameters have been calculated for the outer valence shell orbitals. The assignment of the structure appearing in the experimental photoelectron spectra has been guided by vertical ionization energies and spectral intensities calculated by various theoretical methods that incorporate electron correlation and orbital relaxation. Strong orbital relaxation effects have been found for the 15a', nitrogen lone-pair orbital. The calculations also predict that configuration mixing leads to the formation of several low-lying satellite states. The vibrational structure associated with ionization out of a particular orbital has been simulated within the Franck–Condon model using harmonic vibrational modes. The adiabatic approximation appears to be valid for the X 2A" state, with the β-parameter for this state being independent of the level of vibrational excitation. However, for all the other outer valence ionic states, a disparity occurs between the observed and the simulated vibrational structure, and the measured β-parameters are at variance with the behavior expected at the level of the Franck–Condon approximation. Furthermore, these inconsistencies suggest that the excited electronic states may be interacting vibronically such that the nuclear dynamics occur over coupled potential energy surfaces.

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

The Liquid Jet Endstation for Hard X-ray Scattering and Spectroscopy at the Linac Coherent Light Source

The ability to study chemical dynamics on ultrafast timescales has greatly advanced with the introduction of X-ray free electron lasers (XFELs) providing short pulses of intense X-rays tailored to probe atomic structure and electronic configuration. Fully exploiting the full potential of XFELs requires specialized experimental endstations along with the development of techniques and methods to successfully carry out experiments. The liquid jet endstation (LJE) at the Linac Coherent Light Source (LCLS) has been developed to study photochemistry and biochemistry in solution systems using a combination of X-ray solution scattering (XSS), X-ray absorption spectroscopy (XAS), and X-ray emission spectroscopy (XES). The pump–probe setup utilizes an optical laser to excite the sample, which is subsequently probed by a hard X-ray pulse to resolve structural and electronic dynamics at their intrinsic femtosecond timescales. The LJE ensures reliable sample delivery to the X-ray interaction point via various liquid jets, enabling rapid replenishment of thin samples with millimolar concentrations and low sample volumes at the 120 Hz repetition rate of the LCLS beam. This paper provides a detailed description of the LJE design and of the techniques it enables, with an emphasis on the diagnostics required for real-time monitoring of the liquid jet and on the spatiotemporal overlap methods used to optimize the signal. Additionally, various scientific examples are discussed, highlighting the versatility of the LJE.

EXAFS↗
DOE OSTI2021

Electron-Induced Proton Transfer

The pathway of activationless proton transfer induced by an electron-transfer reaction is studied theoretically. Long-range electron transfer produces highly nonequilibrium medium polarization that can drive proton transfer through an activationless transition during the process of thermalization, dynamically altering the screening of the electron–proton Coulomb interaction by the medium. The cross electron–proton reorganization energy is the main energy parameter of the theory, which exceeds in magnitude the proton-transfer reorganization energy roughly by the ratio of the electron-transfer to proton-transfer distance. This parameter, which can be either positive or negative, is related to the difference in pK a values in two electron-transfer states. The relaxation time of the medium is on the (sub)picosecond time scale, which establishes the characteristic time for activationless proton transfer. Microscopic calculations predict substantial retardation of the collective relaxation dynamics compared to the continuum estimates due to the phenomenology analogous to de Gennes narrowing. As a result, nonequilibrium medium configuration promoting proton transfer can be induced by either thermal or photoinduced charge transfer.

14 SOLAR ENERGY↗
DOE OSTI2021

Numerical investigation of active control of tearing mode by magnetic coils and the role of Δ'

Magnetic feedback stabilization of the tearing mode (TM) is numerically investigated, utilizing the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681) for toroidal tokamak equilibria. With control coil configurations assumed in this study, magnetic feedback partially or fully stabilizes the TM, with either vanishing or finite equilibrium pressure. The best control is achieved by the combination of internal active coils and internal poloidal sensors. The internal and external tearing indices are evaluated for the close-loop system, based on the MARS-F computed mode eigenvalue and eigenfunction, respectively. In the absence of the favorable curvature effect, these two indices are real-valued and quantitatively agree well with each other. For the equilibrium with finite pressure gradient at the mode rational surface, the favorable average curvature effect becomes important and the close-loop tearing index also becomes complex-valued, partly due to interaction of the feedback system with the dissipative wall eddy current response. Isolating the inner layer and outer region response to magnetic feedback, with either proportional or proportional-derivative actions, allows to establish that feedback stabilization of the TM occurs mainly due to modification of the behavior of the external ideal solution, further confirming the analytic result reported in He et al 2021 Phys. Plasmas 28 012504.

Physics↗
DOE OSTI2024

Optimizing Grid-interactive Efficient Building Designs with Stacked Value Streams

Grid-interactive efficient buildings (GEBs) are those characterized by the combination of energy efficiency and demand flexibility with smart technologies and communications to not only deliver greater affordability and comfort to buildings, but also help utilities manage grid operations and lower system costs. This paper presents an innovative techno-economic assessment framework to effectively examine different GEB design options, explore various use cases, define technically achievable benefits, and thereby assist in informed decision-making. In particular, building load flexibility, thermal storage, and battery energy storage are considered. Advanced optimal dispatch problem is formulated to maximize the stacked value streams from multiple, competing use cases, subject to the physical capabilities and operational flexibility associated with different designs and configurations. Comprehensive case studies were performed for a real-world building to evaluate the cost-effectiveness of different GEB designs and offer in-depth insights. It was found that the proposed assessment method could effectively capture the costs and benefits linked to each GEB design option. Furthermore, the study revealed that outage mitigation and demand response are the two most significant sources of benefits for GEBs.

Ma, Xu↗
DOE OSTI2023

The Structural Characterisation and DFT-Aided Interpretation of Vibrational Spectra for Cyclo(l-Cys-d-Cys) Cyclic Dipeptide in a Solid State

Cyclic dipeptides with two intramolecular peptide bonds forming a six-membered 2,5-diketopiperazine ring are gaining significant attention due to their biological and chemical properties. Small changes in the local geometry of such molecules (from cis to trans) can lead to significant structural differences. This work presents the results of a study of cyclo(l-Cys-d-Cys), a dipeptide comprising two cysteine molecules in opposite chiral configurations, with the functional groups situated at both sides of the diketopiperazine ring. X-ray diffraction (XRD) experiment revealed that the molecule crystallises in the $P$-1 space group, which includes the centre of inversion. The IR and Raman vibrational spectra of the molecule were acquired and interpreted in terms of the potential energy distribution (PED) according to the results of density functional theory (DFT) calculations. The DFT-assisted analysis of energy frameworks for the hydrogen bond network within molecular crystals was performed to support the interpretation of X-ray structural data. The optimisation of the computational model based on three-molecule geometry sections from the crystallographic structure, selected to appropriately reflect the intermolecular interactions responsible for the formation of 1D molecular tapes in cyclo(l-Cys-d-Cys) crystal, allowed for better correspondence between theoretical and experimental vibrational spectra. This work can be considered the first complete structural characterisation of cyclo(l-Cys-d-Cys), complemented via vibrational spectroscopy results with full band assignment aided with the use of the DFT method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Unraveling the Atomic Mechanism of the Crystalline Phase‐Dependent Structural Features and Special Spectral Design of α‐, β‐, and Ɛ‐Ga₂O₃

Atomic‐scale phase transformations profoundly influence the functional properties of Ga₂O₃ polymorphs. By combining irradiation experiments with microstructure characterization and theoretical approaches, phase‐specific energy‐dissipation pathways in α‐, β‐, and ε‐Ga₂O₃ are uncovered and strategies for targeted property design are outlined. Competing antiphase boundaries (APBs) and twin domain boundaries (TDBs) promote irreversible α→ε interconversion through domain fragmentation. In β‐Ga₂O₃, defect‐induced stress gradients drive two distinct local transformations: surface Ga‐aggregated β→δ that stabilizes transient states, and latent‐track‐confined β→κ phase transition with recoverable distortions via cation reordering. Under electronic excitation, β‐Ga₂O₃ forms nanohillocks via robust GaO₆ octahedra (high density/strong Ga─O bonds), while α/ε‐Ga₂O₃ generates nanopores from tetrahedral Ga looseness (low bonding energy), highlighting phase‐dependent surface dynamics shaped by atomic packing and bonding anisotropy. Defect‐regulated recombination suppresses visible photoluminescence in α/β‐Ga₂O₃, whereas in ε‐Ga₂O₃ bandgap narrowing of ΔE: 0.30 eV is observed, enhancing emission. Linking phase‐dependent defect‐carrier interactions and metastable‐phase engineering in Ga₂O₃ enables property optimization for power‐electronics and optoelectronics devices.

electronic state configuration↗
DOE OSTI2022

Water Network Shape-Dependence of Local Interactions with the Microhydrated –NO 2 – and –CO 2 – Anionic Head Groups by Cold Ion Vibrational Spectroscopy

We report the structural evolutions of water networks and solvatochromic response of the CH 3 NO 2 – radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH 3 NO 2 – ·(H 2 O) n=1–6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO 2 – ·(H 2 O) n=5–7 (R = CD 3 CD 2 ), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a “book” isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor–acceptor–donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient –XO 2 H–OH–(H 2 O) 2 intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH 3 NO 2 – ·(H 2 O) n=1–6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. Furthermore, this behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

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

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2024

IDAES-PSE 2.6.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.

AS↗
DOE OSTI2023

PERLE: A novel facility for ERL development and applications in multi-turn configuration and high-power regime

The development of ERLs has been recognized as one of the five main pillars of accelerators R&D in support of the European Strategy for Particle Physics (ESPP). The international panel in charge of the ERL Roadmap definition recognized PERLE project as “a central part of the roadmap for the development of energy-recovery linacs”, with milestones to be achieved by the next ESPP in 2026. PERLE project is aiming at the construction of a novel ERL facility for the development and application of the energy recovery technique in multi-turn configuration, high current and large energy regime. It will operate in a 3-turns mode, first at 250 MeV, then upgraded to 500 MeV with 20mA beam current. Such challenging parameters make PERLE a unique multi-turn ERL facility operating at an unexplored operational power regime (10MW), studying and validating a broad range of accelerator phenomena, paving the way for the future larger scale ERLs. PERLE will be the necessary demonstrator for the future HEP machine (LHeC / FCC-eh), with which it shares the same technological choices and beam parameters. Furthermore, PERLE opens a new frontier for the physics of “the electromagnetic probe”. It will be the first ERL dedicated to Nuclear Physics for studying the eN interaction with radioactive nuclei. Here we will report on the project status, introduce the main ongoing achievements and describe the staged strategy we will adopt toward the construction of PERLE machine at its nominal performances.

Kaabi, Walid↗
DOE OSTI2025

Observation of Topological Chirality Switching Induced Freezing of a Skyrmion Crystal

Magnetic skyrmions are topologically protected quasi-particles with a well-defined chirality. Control over their chirality is proposed as an additional feature for encoding data bits or as qubits in quantum computing due to their high efficiency and stability against achiral magnetic textures. Here it is shown that an in-plane magnetic field can be utilized to reshape the energy barriers between different skyrmionic bubbles (e.g., Bloch type, type-II) enabling spontaneous chirality fluctuations with a frequency that increases with the strength of the in-plane field. The insulating van der Waals ferromagnet CrBr3 is used as an archetypal system for low damping, reduced energy dissipation and a high number of magnetic phases to capture the chirality dynamics in real time through cryo-Lorentz transmission electron microscopy. It is observed that the interplay between the intrinsic Dzyaloshinskii–Moriya interaction and out-of-plane field biased the chirality dynamics, favoring one handedness over the other. A remarkable consequence of the spontaneous chirality switching mechanism is that it induces a freezing (or crystallization) process in the skyrmion lattice. As the bubbles fluctuate between Bloch and type-II they elongate and shrink parallel to the in-plane field. Subsequently, the overall lattice crystallizes along the in-plane field direction, inducing a phase transition from a disordered liquid state to a hexatic phase where skyrmions are highly ordered resembling that of a solid. The results indicate chirality as an active element in the creation of topologically protected skyrmion crystals unveiling pathways toward chiral spintronic device platforms with tunable embedded configuration.

2D Phase Transition↗