Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “overlapping”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Response functions and giant monopole resonances for light to medium-mass nuclei from the ab initio symmetry-adapted no-core–shell model

Using the ab initio symmetry-adapted no-core–shell model, we compute sum rules and response functions for light to medium-mass nuclei, starting from interactions that are derived in the chiral effective field theory. Specifically, we investigate electromagnetic transitions of monopole, dipole and quadrupole nature for 4 He, and explore dominant features of giant monopole resonances in symmetric nuclei such as the closed-shell 4 He and 16 O light nuclei, the intermediate-mass open-shell 20 Ne and the medium-mass closed-shell 40 Ca. Furthermore, for the NNLO opt chiral potential, we determine parameter-free monopole sum rules, which can provide information on the incompressibility of symmetric nuclear matter. Here, we report 213(10) MeV as an estimate for the compression modulus for infinite nuclear matter, which overlaps with the lower range of values often used in current astrophysical applications.

ab initio nuclear structure↗
DOE OSTI2021

Low-temperature high-frequency dynamic magnetic susceptibility of classical spin-ice Dy 2 Ti 2 O 7

In this study, radio-frequency (14.6 MHz) AC magnetic susceptibility, , of Dy 2 Ti 2 O 7 was measured using self-oscillating tunnel-diode resonator. Measurements were made with the excitation AC field parallel to the superimposed DC magnetic field up to 5 T in a wide temperature range from 50 mK to 100 K. At 14.6 MHz, a known broad peak of from kHz—range audio-frequency measurements around 15 K for both [111] and [110] directions shifts to 45 K, continuing the Arrhenius activated behavior with the same activation energy barrier of E a ≈ 230 K. Magnetic field dependence of along [111] reproduces previously reported low-temperature two-in-two-out to three-in-one-out spin configuration transition at about 1 T, and an intermediate phase between 1 and 1.5 T. The boundaries of the intermediate phase show reasonable overlap with the literature data and connect at a critical endpoint of the first order transition line, suggesting that these features are frequency independent. An unusual upturn of the magnetic susceptibility at T → 0 was observed in magnetic fields between 1.5 T and 2 T for both magnetic field directions, before fully polarized configuration sets in above 2 T.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Two methods to study inelastic neutron-scattering measurements based on ω n (q) versus S(q, ω) applied to the magnetic open honeycomb lattice Tb 2 Ir 3 Ga 9

This work describes two methods to fit the inelastic neutron-scattering spectrum S(q, ω) with wavevector q and frequency ω. The common and well-established method extracts the experimental spin-wave branches ω n (q) from the measured spectra S(q, ω) and then minimizes the difference between the observed and predicted frequencies. When n branches of frequencies are predicted but the measured frequencies overlap to produce only m < n branches, the weighted average of the predicted frequencies must be compared to the observed frequencies. A penalty is then exacted when the width of the predicted frequencies exceeds the width of the observed frequencies. The second method directly compares the measured and predicted intensities S(q, ω) over a grid {q i , ω j } in wavevector and frequency space. After subtracting background noise from the observed intensities, the theoretical intensities are scaled by a simple wavevector-dependent function that reflects the instrumental resolution. Furthermore, the advantages and disadvantages of each approach are demonstrated by studying the open honeycomb material Tb 2 Ir 3 Ga 9 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ground and excited states of even-numbered Hubbard ring at half-filling: comparison of the extended Gutzwiller approach with exact diagonalization

It remains a great challenge in condensed matter physics to develop a method to treat strongly correlated many-body systems with balanced accuracy and efficiency. We introduce an extended Gutzwiller (EG) method incorporating a manifold technique, which builds an effective manifold of the many-body Hilbert space, to describe the ground- and excited-state properties of strongly correlated electrons. Here we systematically apply an EG projector onto the ground and excited states of a non-interacting system. Diagonalization of the true Hamiltonian within the manifold formed by the resulting EG wavefunctions gives the approximate ground and excited states of the correlated system. To validate this technique, we implement it on even-numbered fermionic Hubbard rings at half-filling with periodic boundary conditions, and compare the results with the exact diagonalization (ED) method. The EG method is capable of generating high-quality ground and low-lying excited state wavefunctions, as evidenced by the high overlaps of wavefunctions between the EG and ED methods. Favorable comparisons are also achieved for other quantities including the total energy, the double occupancy, the total spin and the staggered magnetization. With the capability of accessing the excited states, the EG method can capture the essential features of the one-electron removal spectral function that contains contributions from states deep in the excited spectrum. Finally, we provide an outlook on the application of this method on large extended systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The role of defect charge, crystal chemistry, and crystal structure on positron lifetimes of vacancies in oxides

Density functional theory based positron lifetime (PL) calculations for cation and oxygen monovacancies in a range of oxides—hematite, magnetite, hercynite, and alumina—have been conducted to compare the impact of defect chemistry and crystal structure on the predicted lifetimes. The role of defect charge state has also been examined. A comparison across the same type of crystalline structure but different composition shows that oxygen vacancies only induce a slight increase in the positron-electron overlap and thus barely modify the PL as compared to the bulk. A much more substantial increase of PL is observed for cation monovacancies, regardless of crystal structure or the elemental nature of the vacancy, which we ascribe to an enhanced localization of charge density around the vacant site. The structural and compositional richness of the oxide leads to longer defect PLs, with defected hercynite exhibiting the longest PLs. The charge state of cation monovacancies modifies only by a small percentage the positron localization, relegating to secondary importance the metal defect's oxidation state in modifying the lifetime of positrons within vacancy traps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

First principles band structure of interacting phosphorus and boron/aluminum δ -doped layers in silicon

Silicon can be heavily doped with phosphorus in a single atomic layer (a δ layer), significantly altering the electronic structure of the conduction bands within the material. Recent progress has also made it possible to further dope silicon with acceptor-based δ layers using either boron or aluminum, making it feasible to create devices with interacting δ layers with opposite polarity. Using density functional theory, we calculate the electronic structure of a phosphorus-based δ layer interacting with a boron or aluminum δ layer, varying the distances between the δ layers. At separations 1 nm and smaller, the dopant potentials overlap and largely cancel each other out, leading to an electronic structure closely mimicking intrinsic silicon. At separations greater than 1 nm, the two δ layers behave independently of one another, with an equivalent electronic structure to a p–n diode with an intrinsic layer taking the place of the depletion region. One mechanism for charge transfer between δ layers at larger distances could be tunneling, where we see a tunneling probability exceeding what would be seen for a standard silicon 1.1 eV triangular barrier, indicating that the interaction between delta layers may enhance tunneling compared to a traditional junction.

delta layers↗
DOE OSTI2026

Dual density wave state in liquid Ga

Liquid Ga is considered to be in a mixed state of metallic and covalent local states, and shows complex behavior, similar to those of other polyvalent Group III–V elements on the Periodic Table, such as Si, Ge, Sn, and Bi. We studied liquid Ga by inelastic neutron scattering to obtain the pair-distribution function. We found that the structure of liquid Ga has the medium-range order characterized by two overlapping density waves (DWs), one originating from ionic repulsion and the other due to electronic driving force to create charge DWs. We suggest that the dual DW state is naturally elucidated by the DW theory, rather than the widely used Ornstein–Zernike theory.

atomic structure↗
DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Synergy of coupled kink and tearing modes in fast ion transport

One aspect of the interaction between fast ions and magnetohydrodynamic (MHD) instabilities is the fast ion transport. Coupled kink and tearing MHD instabilities have also been reported to cause fast ion transport. Recently, the 'kick' model has been developed to compute the evolution of the fast ion distribution from neutral beam injection using instabilities as phase-space resonance sources. Overall, the goal of this paper is to utilize the kick model to understand the physics of fast ion transport caused by the coupled kink and tearing modes. Soft x-ray diagnostics are used to identify the mode parameters in the National Spherical Torus Experiment. The comparison of neutron rates measured and computed from time-dependent TRANSP simulation with the kick model shows that the coupling of kink and tearing mode is important in determination of the fast ion transport. The numerical scan of the mode parameters shows that the relative phase of the kink and tearing modes and the overlapping of kink and tearing mode resonances in the phase space can affect the fast ion transport, suggesting that the synergy of the coupled modes may be causing the fast ion transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Assessment of filtered cameras for quantitative 2D analysis of divertor conditions during detachment in JET L-mode plasmas

We report that estimates for 2D distributions of electron temperature, T e , electron density, n e , and atomic deuterium density, n o , in the JET divertor volume have been inferred from deuterium Balmer line intensity ratios obtained from tomographic reconstructions of divertor camera measurements. This enables also investigation of ionization, S ion , and recombination, S rec , rates. The analysis shows a decrease of T e to 0.5–1.0 eV throughout the outer divertor during detachment in low-confinement (L-mode) plasmas. Simultaneously, the high-n e region and the n 0 distribution in the outer divertor are observed to elongate and shift from the outer strike point towards the X-point. The observations are in qualitative agreement and follow the same sequence with modelling predictions of EDGE2D-EIRENE simulations of a density scan. While the method was found to provide good representation of the evolution of volumetric recombination during detachment, in agreement with the simulations, the movement of the ionization front upstream could not be followed due to lack of spatial overlap between the ionization region and the necessary emission distributions. Consequently, the representation of the ionization conditions and the particle balance in the detached outer divertor are compromised.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of the neutral beam injector operational regime on the Alfven eigenmode saturation phase in DIII-D plasma

The aim of this study is to analyze the effect of the neutral beam injector (NBI) operation regime on the saturation phase of the Alfven Eigenmodes (AEs) in DIII-D plasma. Here, the analysis is done using the linear and nonlinear versions of the gyro-fluid code FAR3d. A set of parametric analyses are performed modifying the nonlinear simulation EP $β$ (NBI injection power), EP energy (NBI voltage) and the radial location of the EP density profile gradient (NBI radial deposition). The analysis indicates a transition from the soft (local plasma relaxation) to the hard MHD (global plasma relaxation) limit if the simulation EP $β$ $\unicode{x2A7E}$ $0.02$, leading to bursting MHD activity caused by radial AEs overlapping. MHD bursts cause an enhancement of the EP transport showing ballistic-like features as avalanche-like events. Simulations in the soft MHD limit show an increment of the EP density gradient as the EP $β$ increases. On the other hand, there is a gradient upper limit in the hard MHD limit, consistent with the critical-gradient behavior. AEs induce shear flows and zonal current leading to the deformation of the flux surfaces and the safety factor profile, respectively, particularly strong for the simulation in the hard MHD limit. Simulations in the hard MHD regime show a decrease of the AE frequency in the saturation phase; this is caused by the destabilization of a transitional mode between a $9/3-10/3$ TAE and a $9/3$ RSAE that may explain the AE frequency down-sweeping observed in some DIII-D discharges. Reducing the EP energy in the nonlinear simulations leads to a weakening of the plasma perturbation. On the other hand, increasing the EP energy causes the opposite effect. Nonlinear simulations of off-axis NBI profiles indicate a lower plasma perturbation as the EP density gradient is located further away from the magnetic axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Diagnosing electric and magnetic fields in laser-driven coil targets

Laser-driven capacitor coils are widely used to generate intense magnetic fields for various applications in high-energy-density (HED) physics research. Accurate measurement of the magnetic fields is essential but challenging, due to the overlapping contributions from magnetic and electric fields in proton radiography, which is the primary tool diagnosing the field generation around the coils. In this study, we systematically analyze proton radiographs obtained from laser-driven capacitor-coil targets along two orthogonal axes under various electromagnetic field conditions, including magnetic field only, electric field only, and combined electromagnetic fields. By analyzing key features in the radiographs, we distinguish and characterize the respective contributions from magnetic and electric fields. Using detailed simulations validated by experimental benchmarks, methods to isolate and quantify the magnetic field and electric field are given. The methods are successfully applied to determine the electric current and charge distribution in a double coil configuration. Our findings provide insights into improving the diagnostic capability of proton radiography, potentially leading to more accurate measurements of electromagnetic fields and enhancing the utility of laser-driven capacitor coils in HED experiments.

electromagnetic field diagnostics↗
DOE OSTI2022

2D imaging of absolute methyl concentrations in nanosecond pulsed plasma by photo-fragmentation laser-induced fluorescence

The methyl radical plays a central role in plasma-assisted hydrocarbon chemistry but is challenging to detect due to its high reactivity and strongly pre-dissociative electronically excited states. In this work, we report the development of a photo-fragmentation laser-induced fluorescence (PF-LIF) diagnostic for quantitative 2D imaging of methyl profiles in a plasma. This technique provides temporally and spatially resolved measurements of local methyl distributions, including in near-surface regions that are important for plasma-surface interactions such as plasma-assisted catalysis. The technique relies on photo-dissociation of methyl by the fifth harmonic of a Nd:YAG laser at 212.8 nm to produce CH fragments. These photofragments are then detected with LIF imaging by exciting a transition in the B-X(0, 0) band of CH with a second laser at 390 nm. Fluorescence from the overlapping A-X(0, 0), A-X(1, 1), and B-X(0, 1) bands of CH is detected near 430 nm with the A-state populated by collisional B-A electronic energy transfer. This non-resonant detection scheme enables interrogation close to a surface. The PF-LIF diagnostic is calibrated by producing a known amount of methyl through photo-dissociation of acetone vapor in a calibration gas mixture. We demonstrate PF-LIF imaging of methyl production in methane-containing nanosecond pulsed plasmas impinging on dielectric surfaces. Absolute calibration of the diagnostic is demonstrated in a diffuse, plane-to-plane discharge. Measured profiles show a relatively uniform distribution of up to 30 ppm of methyl. Relative methyl measurements in a filamentary plane-to-plane discharge and a plasma jet reveal highly localized intense production of methyl. The utility of the PF-LIF technique is further demonstrated by combining methyl measurements with formaldehyde LIF imaging to capture spatiotemporal correlations between methyl and formaldehyde, which is an important intermediate species in plasma-assisted oxidative coupling of methane.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ammonia generation in Ns pulse and Ns pulse/RF discharges over a catalytic surface

Plasma-catalytic ammonia synthesis in a ns pulse discharge and a 'hybrid' ns pulse/RF discharge in plane-to-plane geometry is studied by Fourier Transform infrared absorption spectroscopy. The data are taken in a preheated H 2 –N 2 mixture, with and without Ni/γ-Al 2 O 3 or Co/γ-Al 2 O 3 catalyst placed in the discharge section. The measurement results are taken using two different approaches. The first is a 'single-stage' process, where a ns pulse discharge in the H 2 –N 2 mixture is sustained continuously. In this case, the ammonia yield increases slowly, over a period of tens of minutes. The second is a 'two-stage' process, where the catalyst is first activated by the ns pulse discharge sustained in pure nitrogen, and then the activated catalyst is exposed to the H 2 –N 2 flow, with or without the discharge. In this case, a strong overshoot of the NH 3 number density at the reactor exit is detected, by over a factor of two compared to the single-stage process. After the initial overshoot, the ammonia yield gradually decreases to the 'single stage' value (with the discharge on), or to near zero (with the discharge off). The results demonstrate that the ammonia yield in the plasma-catalytic reactor is controlled by the N atom accumulation on the catalyst surface, which reacts with H atoms thermally dissociated on the catalyst or generated in the plasma. The results also show that the plasma-catalytic ammonia yield is significantly higher compared to that in the ns pulse discharge without the catalyst. The accumulation of H atoms on the catalyst, with their subsequent reactions with N atoms generated in the plasma, is of relatively minor importance at the present conditions. An additional series of measurements was made with a sub-breakdown RF waveform overlapped with the ns pulse discharge train, to enhance the vibrational excitation of nitrogen. The ammonia yield measured with the RF waveform added is approximately 20% higher compared to that at the baseline ns pulse discharge conditions, both with and without the catalyst. This effect is weaker compared to that of the catalyst activation by N atoms. Additional data are necessary to isolate the possible effect of the vibrationally excited N 2 molecules on the ammonia synthesis in the plasma catalytic reactions.

surface reactions↗
DOE OSTI2023

Plasma catalysis: separating plasma and surface contributions for an Ar/N 2 /O 2 atmospheric discharge interacting with a Pt catalyst

Atmospheric pressure non-equilibrium plasmas can form nitrogen oxide (NO x ) compounds directly from nitrogen and oxygen without a catalyst, and at lower catalyst temperatures than would be possible without plasma. In this work, the oxidation of plasma-produced NO from an Ar/N 2 /O 2 non-equilibrium atmospheric-pressure plasma-jet (APPJ) over a platinum-on-alumina powder catalyst was investigated with in-situ infrared spectroscopy. Products downstream of the catalyst bed were analyzed along with catalyst surface species. The catalyst was exposed to plasma at both constant temperature and a cyclic temperature ramp in order to study long-lasting and transient surface changes. Primary incident reactive species to the catalyst were assessed to be NO and O 3 . Pt-Al 2 O 3 at 350 °C increased oxidation of NO relative to Al 2 O 3 or an empty chamber. The surface state of Pt-Al 2 O 3 evolves during plasma-effluent exposure and requires upwards of 20 min exposure for stabilization compared to Al 2 O 3 . Once stable surface conditions are achieved, thermal cycling reveals a repeatable hysteresis pattern in downstream products. At low temperature, oxygen and NO x accumulate on the catalyst surface and react at elevated temperatures to form NO 2 . Increasing plasma power and O 2 :N 2 ratio increases the hysteresis of the heating relative to the cooling curves in the pattern of NO 2 formation. The limitation on NO oxidation at high temperatures was assessed to be Pt-O which is depleted as the catalyst is heated. Once stored species have been depleted, NO oxidation rates are determined by incoming reactants. Two overlapping NO oxidation patterns are identified, one determined by surface reactants formed at low temperature, and the other by reactants arriving at the surface at high temperature. As a result, the plasma is responsible for providing the reactants to the catalyst surface, while the catalyst enables reaction at high temperature or storage at low temperature for subsequent reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Improving wafer-scale Josephson junction resistance variation in superconducting quantum coherent circuits

Quantum bits, or qubits, are an example of coherent circuits envisioned for next-generation computers and detectors. A robust superconducting qubit with a coherent lifetime of O (100 µs) is the transmon: a Josephson junction functioning as a non-linear inductor shunted with a capacitor to form an anharmonic oscillator. In a complex device with many such transmons, precise control over each qubit frequency is often required, and thus variations of the junction area and tunnel barrier thickness must be sufficiently minimized to achieve optimal performance while avoiding spectral overlap between neighboring circuits. Simply transplanting our recipe optimized for single, stand-alone devices to wafer-scale (producing 64, 1x1 cm dies from a 150 mm wafer) initially resulted in global drifts in room-temperature tunneling resistance of ± 30%. Inferring a critical current variation from this resistance distribution, we present an optimized process developed from a systematic 38 wafer study that results in < 3.5% relative standard deviation (RSD) in critical current ( ) for 3000 Josephson junctions (both single-junctions and asymmetric SQUIDs) across an area of 49 cm 2 . Looking within a 1x1 cm moving window across the substrate gives an estimate of the variation characteristic of a given qubit chip. Our best process, utilizing ultrasonically assisted development, uniform ashing, and dynamic oxidation has shown = 1.8% within 1x1 cm, on average, with a few 1x1 cm areas having < 1.0% (equivalent to < 0.5%). Such stability would drastically improve the yield of multi-junction chips with strict critical current requirements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Investigations in the tape-to-tape contact resistance and contact composition in superconducting CORC ® wires

Conductor on Round Core (CORC®) wires and cables, constructed from multiple layers of helically wound REBa 2 Cu 3 O 7–δ tapes, are a promising cable technology for high field magnet applications. An important feature of high-temperature superconductor cables is the ability to share current between conductors, allowing current to bypass drops in I c and minimizing the risk of hot spot formation, which could lead to potential burnout in the superconductor. In contrast to stacked-tape cables, which have continuous contact between tapes, in CORC® the transfer points occur at discrete tape crossovers. The tape-to-tape contact resistance, R c , plays a critical role in the current sharing capabilities and current distribution in cables. For the work reported here, special CORC® wires were manufactured using different winding parameters to investigate variations in R c . Variations comprised inclusion of a lubricant, different lubricant conductivity, inclusion of pre-tinning, and heating briefly to melt the solder. Cables were first tested as straight lengths, followed by bending to a 10 cm diameter. In straight cables R c values ranged from 1 to over 1000 μΩ cm 2 , depending on cabling parameters, with the highest values being found for cables made by the present 'standard' process. Bending the cables to a 10 cm diameter decreased R c by a factor 2–5. Tinning with PbSn decreased R c by three orders of magnitude compared to standard CORC® wires, and heat treating wires with tinned conductor resulted in only a small further decrease in R c . Based on the measured R c at an electric field of 1 μV cm –1 the resulting current transfer length between layers can range from a few millimeters to a tens of centimeters. Examination of contacts with a laser confocal microscope showed plastic deformation of the copper at the edges of the contact overlap area, apparently caused by thicker plating at tape edges digging into the copper of neighboring layers. These images reveal that only a fraction of the total contact surface may actually be touching when there is nothing to compensate for height differential. Images of the PbSn coated tapes indicated that application of solder produces a much more uniform contact surface and higher contact area. Furthermore, imaging of CORC® cross-sections confirmed that in the non-tinned cables there are many regions where tapes are not in contact, while in contrast the PbSn cable shows significantly more contact between the tapes. These different imaging techniques reveal that tape surface morphology is a significant parameter in determining R c .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A non-orthogonal variational quantum eigensolver

Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. This strategy doubles the size of the state preparation circuits but not their depth, while adding a small number of additional two-qubit gates relative to standard variational quantum eigensolver. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H4 and the π-system of hexatriene (C6H8).

97 MATHEMATICS AND COMPUTING↗