Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained optimization”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Recent efforts of vapour-phase strategies for EUV resist toward high- and hyper-NA extreme ultraviolet lithography

Extreme ultraviolet lithography (EUVL, λ = 13.5 nm) is critical for sub-1 nm technology nodes but remains constrained by inherent trade-offs among resolution, line-edge roughness (LER), and sensitivity. Stochastic effects originating from photon shot noise, low-energy secondary electron blur, and the random distribution of resist components further limit its advancement toward high numerical-aperture (NA, NA = 0.55) and hyper-NA (NA ≥0.75) EUVL. While the optimization of spin-on chemically amplified resists (CARs) continues, metal-oxide resists (MORs) have emerged as strong candidates for next-generation EUVL by incorporating metals with high EUV absorption coefficients, which enhances both resist sensitivity and etch resistance during pattern transfer. Besides spin-coating, recent advances in vapour-phase techniques, such as vapour-phase infiltration (VPI), chemical vapour deposition (CVD), and molecular atomic layer deposition (MALD), offer promising pathways to achieve new resist platforms, such as dry resists, that satisfy the stringent thickness and uniformity requirements of next-generation EUVL. These methods enable the direct incorporation of metal species into existing resist matrices or the formation of hybrid inorganic–organic resist platforms, thereby improving film uniformity, etch durability, and pattern fidelity while mitigating stochastic defects. This review highlights the latest advancements in vapour-phase-synthesized EUV resists, emphasizing material design, lithographic performance, and the underlying exposure mechanisms. Although still emerging, vapour-phase strategies are paving the way for an all-dry integration framework that could improve EUV patterning workflows and meet the demands of future technology nodes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beam breakup instability studies of powerful energy recovery linac for experiments

The maximum achievable beam current in an energy recovery linac (ERL) is often constrained by beam breakup (BBU) instability. Our previous research highlighted that filling patterns have a substantial impact on BBU instabilities in multipass ERLs. In this study, we extend our investigation to the eight-cavity model of the Powerful ERL for Experiment (PERLE). We evaluate its requirements for damping cavity higher order modes (HOMs) and propose optimal filling patterns and bunch timing strategies. Our findings reveal a significant new insight: while filling patterns are crucial, the timing of bunches also plays a critical role in mitigating HOM beam loading and BBU instability. This previously underestimated factor is essential for effective BBU control. We estimated the PERLE threshold current using both analytical and numerical models, incorporating the designed PERLE HOM dampers. During manufacturing, HOM frequencies are expected to vary slightly. Our study found no significant difference in BBU suppression for relative rms frequency jitters of 0.001, 0.002, and 0.005 for the same HOM. Introducing a jitter of 0.001 into our models, we found that the dampers effectively suppressed BBU instability, achieving a threshold current an order of magnitude higher than the design requirement. Our results offer new insights into ERL BBU beam dynamics and have important implications for the design of future ERLs. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Estimate of the Mass and Radial Profile of the Orphan–Chenab Stream's Dwarf-galaxy Progenitor Using MilkyWay@home

We fit the mass and radial profile of the Orphan–Chenab Stream's (OCS) dwarf-galaxy progenitor by using turnoff stars in the Sloan Digital Sky Survey and the Dark Energy Camera to constrain N-body simulations of the OCS progenitor falling into the Milky Way on the 1.5 PetaFLOPS MilkyWay@home distributed supercomputer. We infer the internal structure of the OCS's progenitor under the assumption that it was a spherically symmetric dwarf galaxy composed of a stellar system embedded in an extended dark matter halo. We optimize the evolution time, the baryonic and dark matter scale radii, and the baryonic and dark matter masses of the progenitor using a differential evolution algorithm. The likelihood score for each set of parameters is determined by comparing the simulated tidal stream to the angular distribution of OCS stars observed in the sky. We fit the total mass of the OCS's progenitor to (2.0 ± 0.3) × 10 7 M ⊙ with a mass-to-light ratio of γ = 73.5 ± 10.6 and (1.1 ± 0.2) × 10 6 M ⊙ within 300 pc of its center. Within the progenitor's half-light radius, we estimate a total mass of (4.0 ± 1.0) × 10 5 M ⊙ . We also fit the current sky position of the progenitor's remnant to be (α, δ) = ((166.0 ± 0.9)°, (–11.1 ± 2.5)°) and show that it is gravitationally unbound at the present time. The measured progenitor mass is on the low end of previous measurements and, if confirmed, lowers the mass range of ultrafaint dwarf galaxies. Our optimization assumes a fixed Milky Way potential, OCS orbit, and radial profile for the progenitor, ignoring the impact of the Large Magellanic Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Overview on Shielding Analyses for the VENUS Instrument at SNS

VENUS, a world‐class versatile neutron imaging instrument, is under construction and is expected to be completed and ready to start commissioning in 2023. The range of cold to epithermal neutrons at SNS will give users of VENUS access to novel imaging methods, as well as to significantly improved existing methods. The instrument is being built on beam line 10 at Spallation Neutron Source (SNS) First Target Station (FTS) facing a decoupled poisoned hydrogen moderator. Instrument design, which includes optics, Front-End components and instrument enclosure components and shape started over 10 years ago and during this time had significant changes. All changes were supported by neutronics analyses to provide adequate shielding for both Front-End and instrument enclosure.Final optics design will provide the field of view (FOV) at the detector position (where the image is formed) to be as high as 0.20 by 0.20 m. Both the beam, which contains a large fraction of high-energy neutrons, and the desire for a large footprint on the detector are challenges for shielding design and budget, because the driving cost for the instrument is the beam line and enclosure shielding. For cost reason VENUS baseline design is optimized with respect to both the instrument cave footprint and its wall thickness, and Front-End shielding is tailored along the beam. Analyses are performed with the Monte Carlo particle transport code MCNPX version 2.7.0 to make choices on materials, thicknesses and configurations of the shielding. Numerous calculations are performed to meet space and coast constrain and to satisfy instrument physics needs and to comply with radiation protection requirements.

Gallmeier, Franz X.↗
DOE OSTI2020

Uncertainty-based weighted least squares density integration for background-oriented schlieren

We propose an improved density integration methodology for Background-Oriented Schlieren (BOS) measurements that overcomes the noise sensitivity of the commonly used Poisson solver. Here, the method employs a weighted least-squares (WLS) optimization of the 2D integration of the density gradient field by solving an over-determined system of equations. Weights are assigned to the grid points based on density gradient uncertainties to ensure that a less reliable measurement point has less effect on the integration procedure. Synthetic image analysis with a Gaussian density field shows that WLS constrains the propagation of random error and reduces it by 80% in comparison to Poisson for the highest noise level. Using WLS with experimental BOS measurements of flow induced by a spark plasma discharge shows a 30% reduction in density uncertainty in comparison to Poisson, thereby increasing the overall precision of the BOS density measurements.

42 ENGINEERING↗
DOE OSTI2023

Nucleation-mediated reshaping of facetted metallic nanocrystals: Breakdown of the classical free energy picture

Shape stability is key to avoiding degradation of performance for metallic nanocrystals synthesized with facetted non-equilibrium shapes to optimize properties for catalysis, plasmonics, and so on. Reshaping of facetted nanocrystals is controlled by the surface diffusion-mediated nucleation and growth of new outer layers of atoms. Here, Kinetic Monte Carlo (KMC) simulation of a realistic stochastic atomistic-level model is applied to precisely track the reshaping of Pd octahedra and nanocubes. Unexpectedly, separate constrained equilibrium Monte Carlo analysis of the free energy profile during reshaping reveals a fundamental failure of the classical nucleation theory (CNT) prediction for the reshaping barrier and rate. Why? Nucleation barriers can be relatively low for these processes, so the system is not locally equilibrated before crossing the barrier, as assumed in CNT. This claim is supported by an analysis of a first-passage problem for reshaping within a master equation framework for the model that reasonably captures the behavior in KMC simulations.

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

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Characterizing Stress Orientations in Southern Kansas

We repot that Induced seismicity predominantly occurs along faults that are optimally oriented to the local principal compressive stress direction, and the characterization of these stress orientations is an important component of understanding seismic hazards. The seismicity rate in southern Kansas rapidly increased in 2013 primarily due to the disposal of large volumes of wastewater into the Arbuckle Group. Previously, local stress orientations in this area were poorly constrained, which limited our understanding of the complex faulting and diverse earthquake mechanisms in this region. We use shear-wave splitting and focal mechanism inversion techniques to create multiple, independent estimates of maximum horizontal stress directions (S Hmax ) and their spatial variation across the study area. We then create an integrated model of stress orientations for southern Kansas and northern Oklahoma using our local results in conjunction with previous, regional stress orientation estimates. We find that S Hmax in both southern Kansas and central Oklahoma exhibits an east-northeast (~N78°E) orientation, and these regions bound a northeast (~N59°E) rotation within a ~20km area in northern Oklahoma near the Nemaha ridge.

58 GEOSCIENCES↗
DOE OSTI2025

Baseflow Identification via Explainable AI With Kolmogorov‐Arnold Networks

Abstract Hydrological models often involve constitutive laws that may not be optimal in every application. We propose to replace such laws with the Kolmogorov‐Arnold networks (KANs), a class of neural networks designed to identify symbolic expressions. We demonstrate KAN's potential on the problem of baseflow identification, a notoriously challenging task plagued by significant uncertainty. KAN‐derived functional dependencies of the baseflow components on the aridity index outperform their original counterparts; they demonstrate that water availability, rather than potential evapotranspiration, drives baseflow by constraining actual evapotranspiration under arid conditions. On a test set, they increase the Nash‐Sutcliffe efficiency (NSE) by 65%, decrease the root mean squared error by 29%, and increase the Kling‐Gupta efficiency by 34%. This superior performance is achieved while reducing the number of fitting parameters from three to two. Next, we use data from 378 catchments across the continental United States to refine the water‐balance equation at the mean‐annual scale. The KAN‐derived equations based on the refined water balance outperform both the current aridity index model, with up to a 105% increase in NSE, and the KAN‐derived equations based on the original water balance. While the performance of our model and tree‐based machine learning methods is similar, KANs offer the advantage of simplicity and transparency and require no specific software or computational tools. This case study focuses on the aridity index formulation, but the approach is flexible and transferable to other hydrological processes. Plain Language Summary Equations used in hydrologic model are often suboptimal, resulting in reduced prediction accuracy and efficiency. We implemented Kolmogorov‐Arnold networks (KAN), a machine learning algorithm for deriving symbolic formulations, to estimate groundwater recharge and showed that it outperforms an existing state‐of‐the‐art semi‐empirical formulation. In hydrology, Nash‐Sutcliffe efficiency (NSE), root mean squared error (RMSE), and Kling‐Gupta efficiency (KGE) are commonly used to evaluate model performance. Higher NSE and KGE values indicate better performance, while lower RMSE values are preferable. Our results show that NSE increased by 71%, RMSE decreased by 32%, and KGE improved by 25%. In addition, KAN identifies an optimal functional form and can be used to derive new analytical formulas using the prior knowledge. The KAN‐inspired equation outperformed the original formulation and reduced the fitting parameters. Furthermore, we refined the water‐balance equation at the mean‐annual scale and showed that, based on the new water‐balance equation, KAN can derive new formulations that are superior to the original aridity index formulations (up to 105% increase in NSE) and KAN‐derived equations based on the original water balance. These findings highlight the significant potential of KAN to advance the scientific understanding of a wide range of hydrologic processes. Key Points Kolmogorov‐Arnold networks (KANs) enhance interpretability of machine‐learned hydrological models KAN‐derived symbolic formulations outperform state‐of‐the‐art semi‐empirical aridity indices KAN‐identified functional form yields an analytical index with fewer fitting parameters and improved performance

baseflow↗
DOE OSTI2022

Representativeness assessment of the pan-Arctic eddy covariance site network and optimized future enhancements

Abstract. Large changes in the Arctic carbon balance are expected as warming linked to climate change threatens to destabilize ancient permafrost carbon stocks. The eddy covariance (EC) method is an established technique to quantify net losses and gains of carbon between the biosphere and atmosphere at high spatiotemporal resolution. Over the past decades, a growing network of terrestrial EC tower sites has been established across the Arctic, but a comprehensive assessment of the network's representativeness within the heterogeneous Arctic region is still lacking. This creates additional uncertainties when integrating flux data across sites, for example when upscaling fluxes to constrain pan-Arctic carbon budgets and changes therein. This study provides an inventory of Arctic (here > = 60∘ N) EC sites, which has also been made available online (https://cosima.nceas.ucsb.edu/carbon-flux-sites/, last access: 25 January 2022). Our database currently comprises 120 EC sites, but only 83 are listed as active, and just 25 of these active sites remain operational throughout the winter. To map the representativeness of this EC network, we evaluated the similarity between environmental conditions observed at the tower locations and those within the larger Arctic study domain based on 18 bioclimatic and edaphic variables. This allows us to assess a general level of similarity between ecosystem conditions within the domain, while not necessarily reflecting changes in greenhouse gas flux rates directly. We define two metrics based on this representativeness score: one that measures whether a location is represented by an EC tower with similar characteristics (ER1) and a second for which we assess if a minimum level of representation for statistically rigorous extrapolation is met (ER4). We find that while half of the domain is represented by at least one tower, only a third has enough towers in similar locations to allow reliable extrapolation. When we consider methane measurements or year-round (including wintertime) measurements, the values drop to about 1/5 and 1/10 of the domain, respectively. With the majority of sites located in Fennoscandia and Alaska, these regions were assigned the highest level of network representativeness, while large parts of Siberia and patches of Canada were classified as underrepresented. Across the Arctic, mountainous regions were particularly poorly represented by the current EC observation network. We tested three different strategies to identify new site locations or upgrades of existing sites that optimally enhance the representativeness of the current EC network. While 15 new sites can improve the representativeness of the pan-Arctic network by 20 %, upgrading as few as 10 existing sites to capture methane fluxes or remain active during wintertime can improve their respective ER1 network coverage by 28 % to 33 %. This targeted network improvement could be shown to be clearly superior to an unguided selection of new sites, therefore leading to substantial improvements in network coverage based on relatively small investments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Thermal Hall conductivity of electron-doped cuprates

Measurements of the thermal Hall conductivity in hole-doped cuprates have shown that phonons acquire chirality in a magnetic field both in the pseudogap phase and in the Mott insulator state. The microscopic mechanism at play is still unclear. A number of theoretical proposals are being considered including skew scattering of phonons by various defects, the coupling of phonons to spins, and a state of loop-current order with the appropriate symmetries, but more experimental information is required to constrain theoretical scenarios. Here we present our study of the thermal Hall conductivity κ xy in the electron-doped cuprates Nd 2–x Ce x CuO 4 and Pr 2–x Ce x CuO 4 for dopings across the phase diagram, from x = 0 in the insulating antiferromagnetic phase up to x = 0.17 in the metallic phase above optimal doping. We observe a large negative thermal Hall conductivity at all dopings in both materials. Since heat conduction perpendicular to the CuO 2 planes is dominated by phonons, the large thermal Hall conductivity we observe in electron-doped cuprates for a heat current in that direction must also be due to phonons, as in hole-doped cuprates. However, the degree of chirality, measured as the ratio |κ xy /κ xx | where κ xx is the longitudinal thermal conductivity, is much larger in the electron-doped cuprates. Here, we discuss various factors that may be involved in the mechanism that confers chirality to phonons in cuprates, including short-range spin correlations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Physics-informed State-space Neural Networks for transport phenomena

This work introduces Physics -informed State -space neural network Models (PSMs), a novel solution to achieving real-time optimization, flexibility, and fault tolerance in autonomous systems, particularly in transportdominated systems such as chemical, biomedical, and power plants. Traditional data -driven methods fall short due to a lack of physical constraints like mass conservation; PSMs address this issue by training deep neural networks with sensor data and physics -informing using components' Partial Differential Equations (PDEs), resulting in a physics -constrained, end -to -end differentiable forward dynamics model. Further, through two in silico experiments - a heated channel and a cooling system loop - we demonstrate that PSMs offer a more accurate approach than a purely data -driven model. In the former experiment, PSMs demonstrated significantly lower average root -mean -square errors across test datasets compared to a purely data -driven neural network, with reductions of 44 %, 48 %, and 94 % in predicting pressure, velocity, and temperature, respectively. Beyond accuracy, PSMs demonstrate a compelling multitask capability, making them highly versatile. In this work, we showcase two: supervisory control of a nonlinear system through a sequentially updated state -space representation and the proposal of a diagnostic algorithm using residuals from each of the PDEs. The former demonstrates PSMs' ability to handle constant and time -dependent constraints, while the latter illustrates their value in system diagnostics and fault detection.

42 ENGINEERING↗
DOE OSTI2021

Resilient NdFeB magnet recycling under the impacts of COVID-19 pandemic: Stochastic programming and Benders decomposition

Neodymium-iron-boron (NdFeB) magnets are the most powerful magnets per unit volume sold in the commercial market. Despite the increasing demand for clean energy applications such as electric vehicles and wind turbines, disruptive events including the COVID-19 pandemic have caused significant uncertainties in the supply and demand for NdFeB magnets. Therefore, this study aims to alleviate the risk of supply shortage for NdFeB magnets and the containing critical materials, rare-earth elements (REEs), through the development of a resilient reverse supply chain and logistics network design. We develop scenarios to model the unique impact of the COVID-19 pandemic on the proposed business, incorporating both disruption intensity and recovery rate. We formulate a chance-constrained two-stage stochastic programming model to maximize the profit while guaranteeing the network resiliency against disruption risks. To solve the problem in large-scale instances, we develop an efficient Benders decomposition algorithm that reduces the computational time by 98.5% on average compared to the default CPLEX algorithm. When applied to the United States, the model suggests the optimal facility locations, processing capacities, inventory levels, and material flows for NdFeB magnet recyclers that could meet 99.7% of the demand. To the best of our knowledge, this study is the first to incorporate the impacts of the COVID-19 pandemic to design a resilient NdFeB magnet recycling supply chain and logistics network, leveraging risk-averse stochastic programming.

42 ENGINEERING↗
DOE OSTI2021

BEAM COLLIMATION IN THE PIP-II LINAC TO BOOSTER TRANSFER LINE

The new PIP-II superconducting linac will deliver a 2 mA average H⁻ beam to the existing Booster synchrotron. The injected beam is accumulated by charge exchange over approximately 300 turns; phase space painting is used to mitigate space charge effects. To limit the power load on the internal waste beam absorber from the transverse tails of the H⁻ distribution missing the foil, the beam will be collimated in both planes in the linac to Booster transfer line using compact collimators of a novel design. Both the number of parasitic hits and the fraction of the beam missing the foil are sensitive functions of the H⁻ beam centroid position with respect to the edge of the foil. The positioning of the collimation is constrained by the availability of suitable space in the transfer line lattice, by specifics of the collimator design, by the phase space orientation at the collimator, and by the betatron phase advance to the foil needed to achieve proper orientation of the spatial distribution at the injection point. In this contribution, we describe the procedure by which collimator positions were optimized. We then discuss the expected performance of the overall system.

Johnson, David E.↗
DOE OSTI2022

SpectralFly: Ramanujan Graphs as Flexible and Efficient Interconnection Networks

In recent years, graph theoretic considerations have become increasingly important in the design of HPC interconnection topologies. One approach is to seek optimal or near-optimal families of graphs with respect to a particular graph theoretic property, such as diameter. For example, the SlimFly topology is based on a construction of McKay, Miller, and \v{S}ir\'{a}\v{n} which produces a diameter two graph on a number of nodes approaching the Moore bound, i.e. the largest possible diameter two graph with a fixed radix. Motivated by recent work of Aksoy, Bruillard, Young, and Raugas, we consider topologies which optimize the spectral gap rather than the diameter. In particular, we introduce a novel HPC topology, SpectralFly, designed around the Ramanujan graph construction of Lubotzky, Phillips, and Sarnak (LPS). In this work, we show that the combinatorial properties, such as diameter, bisection bandwidth, average path length, and resilience to link failure, of SpectralFly topologies are better than, or comparable to, similarly constrained DragonFly, SlimFly, and BundleFly topologies. Additionally, we simulate the performance of SpectralFly topologies on a representative sample of physics-inspired HPC workloads using the Structure Simulation Toolkit Macroscale Element Library simulator and demonstrate considerable benefit to using LPS construction as the basis of the SpectralFly topology.

graphs and networks, network topology, interconnec↗
DOE OSTI2021

Gleipnir: toward practical error analysis for Quantum programs

Practical error analysis is essential for the design, optimization, and evaluation of Noisy Intermediate-Scale Quantum(NISQ) computing. However, bounding errors in quantum programs is a grand challenge, because the effects of quantum errors depend on exponentially large quantum states. In this work, we present Gleipnir, a novel methodology toward practically computing verified error bounds in quantum programs. Gleipnir introduces the (ρ,δ)-diamond norm, an error metric constrained by a quantum predicate consisting of the approximate state ρ and its distance δ to the ideal state ρ. This predicate (ρ,δ) can be computed adaptively using tensor networks based on the Matrix Product States. Gleipnir features a lightweight logic for reasoning about error bounds in noisy quantum programs, based on the (ρ,δ)-diamond norm metric. Furthermore, our experimental results show that Gleipnir is able to efficiently generate tight error bounds for real-world quantum programs with 10 to 100 qubits, and can be used to evaluate the error mitigation performance of quantum compiler transformations.

Tao, Runzhou↗
DOE OSTI2021

Bandwidth extension and conversion efficiency improvements beyond phase matching limitations using cavity-enhanced OPCPA

The conversion efficiency and phase matching bandwidth of ultrafast optical parametric amplification (OPA) are constrained by the dispersion and nonlinear coefficient of the employed crystal as well as pulse shaping effects. In our work we show that an enhancement cavity resonant with the pump seeded at the full repetition rate of the pump laser can automatically reshape the small-signal gain in optical parametric chirped-pulse amplification (OPCPA) to achieve close-to-optimal operation. This new method termed cavity-enhanced OPCPA or C-OPCPA significantly increases both the gain bandwidth and the conversion efficiency, in addition to boosting gain for high-repetition-rate amplification. The goal in C-OPCPA is to arrive at a condition of impedance matching at all temporal coordinates, such that, in the absence of linear losses, all the incident pump power is dissipated in the nonlinear loss element, i.e., converted to signal and idler. The use of a low finesse enhancement cavity resonant with a low average power (<1W) and a high repetition rate (78MHz) pump source is shown to achieve more than 50% conversion efficiency into signal and idler from the coupled pump in an optical parametric process, whereas an equivalent amount of pump power in a single-pass configuration leads to negligible conversion. Additionally, the gain bandwidth is extended by a factor of 3-4 beyond the phase-matching limit. Our empirical observations are corroborated by a numerical analysis of depletion optimizing the single-pass case, which assesses the underlying impedance matching that is responsible for the observed performance improvements.

47 OTHER INSTRUMENTATION↗