Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator study”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Imaging the kinetics of anisotropic dissolution of bimetallic core–shell nanocubes using graphene liquid cells

Chemical design of multicomponent nanocrystals requires atomic-level understanding of reaction kinetics. Here, we apply single-particle imaging coupled with atomistic simulation to study reaction pathways and rates of Pd@Au and Cu@Au core-shell nanocubes undergoing oxidative dissolution. Quantitative analysis of etching kinetics using in situ transmission electron microscopy (TEM) imaging reveals that the dissolution mechanism changes from predominantly edge-selective to layer-by-layer removal of Au atoms as the reaction progresses. Dissolution of the Au shell slows down when both metals are exposed, which we attribute to galvanic corrosion protection. Morphological transformations are determined by intrinsic anisotropy due to coordination-number-dependent atom removal rates and extrinsic anisotropy induced by the graphene window. Our work demonstrates that bimetallic coreshell nanocrystals are excellent probes for the local physicochemical conditions inside TEM liquid cells. Furthermore, single-particle TEM imaging and atomistic simulation of reaction trajectories can inform future design strategies for compositionally and architecturally sophisticated nanocrystals.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Data-Centric Development of Lignin Structure–Solubility Relationships in Deep Eutectic Solvents Using Molecular Simulations

Lignin is a natural source of aromatic chemicals with significant potential as an abundant, renewable feedstock for value-added products. Deep eutectic solvents (DES)–solvents composed of a hydrogen bond donor (HBD) and acceptor (HBA) in varying ratios–have emerged as a highly tunable class of solvents for lignin solubilization. However, the variety of possible DES compositions and limited molecular-scale understanding of lignin solubility makes solvent selection a challenge without laborious trial-and-error experimentation. To address these challenges, we use classical molecular dynamics (MD) simulations to study the interactions of lignin model compounds with various DES–water systems. Quantitative parameters (descriptors) were calculated by postprocessing the MD results and used to train a regression model that predicts experimentally determined solubilities of lignin model compounds. This approach revealed that the most important descriptors of solubility are the system temperature, solute hydrophilicity, and metrics quantifying hydrogen bonding. Maximizing the interactions between solute–HBD (hydrophobic group), water–HBD (hydrophilic group), and water–HBA molecules led to the highest model compound solubility. Our results support a hydrotropic mechanism in which extensive DES–water hydrogen bonding and favorable HBD interactions with the solute promote high solubility. We applied the regression model derived using model compounds to predict the solubility of representative lignin oligomers. The model predicted lignin oligomers’ solubilities in good agreement with experiments, indicating that the simulations of model compounds can be extended to predict the solubility of larger lignin compounds across a range of solvent compositions and temperatures. Furthermore, these findings provide new molecular-scale insight into lignin solubilization mechanisms and a new method for computationally screening potential solvent systems for lignin valorization.

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

Effect of Sulfonation Level on the Percolated Morphology and Proton Conductivity of Hydrated Fluorine-Free Copolymers: Experiments and Simulations

Using all-atom molecular dynamics simulations and a variety of experimental methods, we previously reported on a linear polyethylene with pendant phenyl sulfonated groups precisely on every fifth carbon along the backbone. With increasing relative humidity this fluorine-free polymer self-assembled to form nanoscale water channels and exhibited exceptional proton conductivity. Expanding upon those findings, here we explore partially sulfonated random copolymers, referred to as p 5PhSH-Y. Using either acetyl sulfate or sulfuric acid, a wide range of sulfonation levels were prepared ( Y = 34−98%) corresponding to ion-exchange capacities (IEC) of 2.0−4.4 mmol/g. Combining experimental techniques and all-atom molecular dynamics simulations, we study the effect of Y on water uptake, nanoscale morphology, and the proton/water transport properties of p5PhSH- Y . The proton conductivity of p 5PhSH- Y increases with relative humidity and with Y and achieves values in excess of 0.1 S/cm. These high conductivities are attributed to high IEC and welldeveloped nanoscale percolated hydrophilic domains made possible by the flexible backbone. We quantitatively describe the nature of the water channels using the characteristic distance, channel width distribution, the area per sulfonate group at the hydrophilic/ hydrophobic interface, and the fractal dimension. Notably, the channel widths and the areas per sulfonate group are nominally independent of the level of sulfonation, while depending significantly on the level of hydration. The fractal dimension of the water channels correlates strongly with the water diffusion coefficients calculated from the molecular dynamics (MD) simulations. These findings demonstrate that the p 5PhSH- Y hydrocarbon copolymers can be modified to tune properties, particularly proton conductivity.

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

The Future Electron-Ion Collider

Generalized Parton Distributions include rich information and became a powerful tool for studying hadron structure. Deeply Virtual Compton Scattering is the golden channel to access GPDs. The development of high luminosity and high-acceptance detectors (Electron-Ion Collider) allows physicists to overcome the difficulty of DVCS measurements. The outstanding performance of the accelerator performs the eA collisions at the center of mass energy from 20 to 140 GeV, while the luminosity at the scale of 1034 cm?2s?1. The new type of high-density crystal, PWO-II, produced by CRYTUR in 2×2×20 cm3 will be utilized in the calorimeter. The transparency of PWO-II crystals is > 70% at 620 nm, > 60% at 420 nm, and >35% at 360 nm. The crystals also show good radiation hardness under 30 Gy radiation exposure. Ten PWO-II crystals have a uniform light yield at 30 p.e./MeV, providing sufficient light yield to reduce the fluctuation of the energy measurements. Near 3000 crystals were constructed inside the 12-sided polygon supporting structure in the simulation. The island clustering algorithm was used for reconstructed the energy. The energy resolution study by the particle gun shows the stochastic term and constant term are 1.8% and 1.2%, respectively. The spatial resolution of NEEMC varies from 5% to 15% of the crystal?s width depending on the particle?s incident angle. The pion rejection of NEEMC can reach 103 with electron efficiency > 85% when the particle?s energy is larger than 1GeV. The Pi0- identified efficiency study can be interpreted as finding the local maxima in the single cluster caused by two high energy close photons. The study results show that efficiency is nearly 100% for pi0 energy is smaller than 10 GeV, and efficiency drops to 30% with 20 GeV pi0. The new type of 3x3 pixelated AC-LGAD were wire bonded to ALTIROC for performance test. The cross-talk between the adjacent channels is about 20% for VPA and 10% for TZ. A series of the TDC characteristic measurements show that the jitter for both preamplifiers is about 20 ps, and the time-walk effect is mild for inject charge > 12 pF. Furthermore, the beta source radiation results quantify the sharing scale of the 3x3 pixels AC-LGAD (? 20%). The electronics simulation study results show no significant changes in the spatial resolution, whether ADC resolution is 8, 10, or 12 bits. The 8-bit ADC is decided to use in the EICROC as it has a smaller size and power consumption than the 10-bit and 12-bit ADC. The ECCE is one of the full detector proposals for EIC. The new type of DVCS generator, called TOPEG, is used to generate the high acceptance beam configurations of 18×110 GeV2 electron and 4He. The acceptance study shows a -1.8 < ? < -1.4 gap between the BEMC and EEMC. The 10?x,y geometry cut is applied on the Roman Pots, so the acceptance of Roman Pots quickly drops to 0% when the polar angle of the recoiled 4He < 2 mrad. This primary ECCE simulation study suggests extending the acceptance of BEMC longitudinally as the structure limits the radial size of EEMC. The acceptance of the Roman Pots is still challenging. In 2023, The overall design of ePIC was finalized by merging two full detector proposals (ECCE and ATHENA) after a series of intensive simulation studies.

Wang, Pu-Kai↗
DOE OSTI2024

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The Southern hemisphere midlatitude circulation response to rapid adjustments and sea surface temperature driven feedbacks

Rapid adjustments – the response of meteorology to external forcing while sea surface temperatures (SST) and sea ice are held fixed – can affect the midlatitude circulation and contribute to long-term forced circulation responses in climate simulations. This study examines rapid adjustments in the Southern hemisphere (SH) circulation using nine models from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), which perform fixed SST and coupled ocean experiments for five perturbations: a doubling of carbon dioxide (2xCO2), a tripling of methane (3xCH4), a five-fold increase in sulphate aerosol (5xSO4), a ten-fold increase in black carbon aerosol (10xBC) and a 2% increase in solar constant (2%Sol). In the coupled experiments, the SH eddy-driven jet shifts poleward and strengthens for forcings that produce global warming (and vice versa for 5xSO4), with the strongest response found in austral summer. In austral winter, the responses project more strongly onto a change in jet strength. For 10xBC, which induces strong shortwave absorption, the multi-model mean (MMM) rapid adjustment in DJF jet latitude is ~75% of the change in the coupled simulations. For the other forcings, which induce larger SST changes, the effect of SST-mediated feedbacks on the SH circulation is larger than the rapid adjustment. Nevertheless, for these perturbations the magnitude of the MMM jet shift due to the rapid adjustment is still around 25-33% of that in the coupled experiments. The results demonstrate the need to understand the mechanisms for rapid adjustments in the midlatitude circulation, in addition to the effect of changing SSTs.

Wood, Thomas↗
DOE OSTI2022

Impact of Changes in ENDF/B-VII.1 and ENDF/B-VIII.0 235U Nuclear Data Indicated by NDSE Studies on LLNL Pulsed Sphere Simulations

A recent journal article by J.A. Gomez et al. measured gamma-ray die-off curves of three subcritical static highly-enriched uranium (HEU) assemblies driven by an external neutron source with a new detector system. This new detector system was developed for being used in dynamically driven subcritical assemblies as part of the Neutron Diagnosed Subcritical Experiments (NDSE) program. Simulations of these die-off curves with various nuclear data and comparison to experimental data indicated that a decrease of the ENDF/B-VII.1 235 U(n,inl) cross section by a factor 0.8 and 0.85 for ENDF/B-VIII.0 would lead to better predictions of experimental data. Here, we test the proposed changes with another type of measurement response, namely neutron-leakage spectra emitted in LLNL pulsed sphere measurements. These spheres were pulsed by 14-MeV neutrons produced via the D+T reaction in their center. The proposed changes in nuclear data have a distinctly smaller impact on predicting pulsed-sphere neutron-leakage spectra than for the die-off curves; they lead to a worsened prediction of the inelastic valley of LLNL pulsed-sphere neutron spectra indicating that the proposed change could constitute a compensating error. Changes in the inelastic angular distributions along with the cross section might be worthwhile to study

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Evaluating indoor thermal resilience of passive and low-power cooling shelters for outdoor workers in India

India’s 231 million outdoor workers are exposed to deadly heat during heat waves. We report results from a simulation-based study that focuses on passive design of cooling shelters for the hot and dry climate zone of India. The goal is to design a cooler-than-outdoor shelter for outdoor workers’ intermittent rest and recovery from heat stress during their arduous outdoor work, to avoid heat-related morbidities and mortality. Expected indoor thermal conditions in low-cost cooling shelters deploying a range of passive designs and low-power active measures are investigated, consistent with the current Indian standards and practices, located in the city of Jodhpur in western India. Heat resilience of the shelter was measured with the predicted reduction of hazardous hours inside the shelter based on the wet-bulb globe temperature (WBGT) and extended heat index. These measures include improving the building envelope, installing cool roofs, using internal thermal mass, and natural ventilation. Additionally, a low-power active measure of using ventilation fans and ceiling fans is explored. EnergyPlus-based simulation results show that simple, commonly available measures, particularly natural ventilation and ceiling fans, achieve the most significant reductions in overheating danger hours, by 22% and 21%, respectively, while more complex or costly passive strategies yield only marginal additional benefits. Moreover, combined measure packages are identified that can reduce the indoor WBGT by 8 to 10 °C, eliminating overheating danger hours. These findings highlight the critical importance of scalable, low-cost, and easily deployable cooling solutions for developing resilient cooling shelters for vulnerable outdoor workers during heat waves in India.

building performance↗
DOE OSTI2025

Micro-Mechanically Guided High-Throughput Alloy Design Exploration Towards Metastability-Induced H Embrittlement Resistance

We develop a high-throughput approach for studying H embrittlement (HE)-resistance in alloys, which is based on combinatorial compositional screening of metastability effects by in situ scanning electron microscopy H-analyses. The project objective included: (i) technique development of high-throughput screening (HTS) for HE-resistance, (ii) discovery of new metallic materials with superior HE-resistance, (iii) Multiscale verification of HE-resistance of the new alloys and H-barrier layers, from atomic scale to an engineering scale. The investigation focused on a model system that enables exploration of different metastable states in Fe-based complex-concentrated alloys. Composition spread islands with hundreds of varying compositions are fabricated on a single substrate using a combinatorial co-sputtering technique. To screen these alloys, we developed and employed a home-built SEM-based integrated analysis system capable of characterizing H-permeability, H-trapping, H-influence on mechanical properties, and other material properties, as well as atomistic to continuum simulations to study the underlying physics.

08 HYDROGEN↗
DOE OSTI2025

Hybrid quantum simulations with qubits and qumodes on trapped-ion platforms

We explore the feasibility of gate-based hybrid quantum computing using both discrete (qubit) and continuous (qumode) variables on trapped-ion platforms. Trapped-ion systems have demonstrated record one- and two-qubit gate fidelities and long qubit coherence times, while qumodes, which can be represented by the collective vibrational modes of the ion chain, have remained relatively unex- plored for their use in computing. Using numerical simulations, we show that high-fidelity hybrid gates and measurement operations can be achieved for existing trapped-ion quantum platforms. As an exemplary application, we consider quantum simulations of the Jaynes-Cummings-Hubbard model, which is given by a one-dimensional chain of interacting spin and boson degrees of free- dom. Using classical simulations, we study its real-time evolution and develop a suitable variational quantum algorithm for ground state preparation. Furthermore, our results motivate further studies of hybrid quantum computing in this context, which may lead to direct applications in condensed matter and fundamental particle and nuclear physics.

Lower-dimensional field theories↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2024

Evaluation of precipitation across the contiguous United States, Alaska, and Puerto Rico in multi-decadal convection-permitting simulations

Abstract This study is an early effort to generate a multi-decadal convection-permitting regional climate dataset that covers nearly the entire North American continent. We assessed a 20 year dynamically downscaled regional climate simulation at a 4 km spatial resolution with explicit convection across the contiguous United States (CONUS), Alaska, and Puerto Rico. Specifically, we evaluated the model’s performance in representing mean, 95th percentile, and extreme precipitation across regions. Our findings indicate that when compared with ERA5 reanalysis, the forcing data, convection-permitting simulation improves representations of seasonal, 95th percentile, and extreme precipitation over a large portion of the CONUS, Alaska, and Puerto Rico, particularly in areas where precipitation is heaviest. The simulation adds value over its forcing data (ERA5) in up to 53% of all grid cells in the CONUS, 68.8% in Alaska, and 84.0% in Puerto Rico. It is important to note that, however, despite improvements, model errors in Puerto Rico remain large. Similar improvements are observed in extreme indices, including consecutive dry days, maximum 5 days precipitation, and extreme precipitation. Analysis of the diurnal cycle of mean hourly precipitation suggests that representations of convective processes—including onset, dissipation, suppression, downstream propagation, and local circulation—improved overall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sequence-defined Pareto frontier of a copolymer structure

The correlations between the sequence of monomers in a macromolecule and its three-dimensional (3D) structure is a grand challenge in polymer science. The properties and functions of macromolecules depend on their 3D shape that has appeared to be dictated by their monomer sequence. However, the progress towards understanding the sequence–structure-property correlations and their utilization in materials engineering are slow because it is almost impossible to characterize an astronomically large number of possible sequences of a copolymer using traditional experimental and simulation methods. To address this problem, here, we combine evolutionary computing and coarse-grained molecular dynamics (CGMD) simulation and study the sequence-structure correlations of a model AB-type copolymer in a solution and assess the impact of sequence on the packing density in its bulk phase. The CGMD-based evolutionary algorithm (EA) screens the sequence space of a single chain copolymer efficiently and identifies a wide range of single-molecule structures including extremal radii of gyration. The data are utilized to estimate the Pareto front of the structure-space of a binary copolymer as a function of its composition. The monomer packings in single-molecule solution phase and multimolecular bulk phase are found to be identical. Finally, this work highlights the opportunities of sequence-specific control of macromolecular structure for designing target materials.

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

Reversible Surface Energy Storage in Molecular-Scale Porous Materials

Forcible wetting of hydrophobic pores represents a viable method for energy storage in the form of interfacial energy. The energy used to fill the pores can be recovered as pressure–volume work upon decompression. For efficient recovery, the expulsion pressure should not be significantly lower than the pressure required for infiltration. Hysteresis of the wetting/drying cycle associated with the kinetic barrier to liquid expulsion results in energy dissipation and reduced storage efficiency. In the present work, we use open ensemble (Grand Canonical) Monte Carlo simulations to study the improvement of energy recovery with decreasing diameters of planar pores. Near-complete reversibility is achieved at pore widths barely accommodating a monolayer of the liquid, thus minimizing the area of the liquid/gas interface during the cavitation process. At the same time, these conditions lead to a steep increase in the infiltration pressure required to overcome steric wall/water repulsion in a tight confinement and a considerable reduction in the translational entropy of confined molecules. In principle, similar effects can be expected when increasing the size of the liquid particles without altering the absorbent porosity. While the latter approach is easier to follow in laboratory work, we discuss the advantages of reducing the pore diameter, which reduces the cycling hysteresis while simultaneously improving the stored-energy density in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Rotational transition, domain formation, dislocations, and defects in vortex systems with combined sixfold and twelvefold anisotropic interactions

We introduce a phenomenological model for a pairwise repulsive interaction potential of vortices in a type-II superconductor, consisting of superimposed sixfold and twelvefold anisotropies. Using numerical simulations we study how the vortex lattice configuration varies as the magnitudes of the two anisotropic interaction terms change. A triangular lattice appears for all values, and rotates through 30 ° as the ratio of the sixfold and twlevefold anisotropy amplitudes is varied, in agreement with experimental results. The transition causes the vortex lattice to split into domains that have rotated clockwise or counterclockwise, with grain boundaries that are “decorated” by dislocations consisting of fivefold and sevenfold coordinated vortices. We also find intradomain dislocations and defects, and characterize them in terms of their energy cost. We discuss how this model could be generalized to other particle-based systems with anisotropic interactions, such as colloids, and consider the limit of very large anisotropy where it is possible to create cluster crystal states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

New interaction potentials for alkaline earth silicate and borate glasses

New interaction potentials were developed for molecular dynamics simulations to study the role of Mg and Ca in modifying the structure and properties of alkaline earth silicates and borates. Competition between the depolymerization of the silica network and the formation of new bonds between oxygen atoms and modifiers leads to the enhancement of the elastic moduli with increasing modifier content in alkaline earth silicate glasses. Compared with calcium silicate, the higher elastic moduli of magnesium silicate result from a higher connectivity of the overall glass network due to the incorporation of fourfold coordinated magnesium and a more rigid connection between oxygen atoms and modifiers. In contrast to the silicates, the effect of modifier on the elastic moduli of alkaline earth borates is dominated by the formation of fourfold coordinated boron (N 4 ). Calcium borate with higher N 4 shows a more rigid network structure and higher elastic moduli.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tidal Stripping of a White Dwarf by an Intermediate-mass Black Hole

During the inspiralling of a white dwarf (WD) into an intermediate-mass black hole (~10 2-5 M ⊙ ), both gravitational waves (GWs) and electromagnetic (EM) radiation are emitted. Once the eccentric orbit's pericenter radius approaches the tidal radius, the WD would be tidally stripped upon each pericenter passage. The accretion of this stripped mass would produce EM radiation. It is suspected that the recently discovered new types of transients, namely the quasiperiodic eruptions and the fast ultraluminous x-ray bursts, might originate from such systems. Modeling these flares requires a prediction of the amount of stripped mass from the WD and the details of the mass supply to the accretion disk. We run hydrodynamical simulations to study the orbital parameter dependence of the stripped mass. We find that our results match the analytical estimate that the stripped mass is proportional to z 5/2 , where z is the excess depth by which the WD overfills its instantaneous Roche lobe at the pericenter. The corresponding fallback rate of the stripped mass is calculated, which may be useful in interpreting the individual flaring light curve in candidate EM sources. We further calculate the long-term mass-loss evolution of a WD during its inspiral and the detectability of the GW and EM signals. The EM signal from the mass-loss stage can be easily detected: the limiting distance is ~320(M h /10 4 M ⊙ ) Mpc for the Einstein Probe. The GW signal, for space-borne detectors such as Laser Interferometer Space Antenna or TianQin, can be detected only within the Local Supercluster (~33 Mpc).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Adaptive variational ground state preparation for spin-1 models on qubit-based architectures

Here, we apply the adaptive variational quantum imaginary time evolution (AVQITE) method to prepare ground states of one-dimensional spin S = 1 models. We compare different spin-to-qubit encodings (standard binary, Gray, unary, and multiplet) with regard to the performance and quantum resource cost of the algorithm. Using state-vector simulations, we study two well-known spin-1 models: the Blume-Capel model of transverse-field Ising spins with single-ion anisotropy, and the XXZ model with single-ion anisotropy. We consider system sizes of up to 20 qubits, which corresponds to spin-1 chains up to length 10. We determine the dependence of the number of CNOT gates in the AVQITE state preparation circuit on the encoding, the initial state, and the choice of operator pool in the adaptive method. Independent of the choice of encoding, we find that the CNOT gate count scales cubically with the number of spins for the Blume-Capel model and quartically for the anisotropic XXZ model. However, the multiplet and Gray encodings present smaller prefactors in the scaling relations. These results provide useful insights for the implementation of AVQITE on quantum hardware.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗