SEARCH · Engineering Papers

Results for “arbitrary choices”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

48 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Twisted bilayer graphene. II. Stable symmetry anomaly

We show that the entire continuous model of twisted bilayer graphene (TBG) (and not just the two active bands) with particle-hole symmetry is anomalous and hence incompatible with lattice models. Previous works [e.g., Z. Song, Z. Wang, W. Shi, G. Li, C. Fang, and B. A. Bernevig, Phys. Rev. Lett. 123, 036401 (2019); J. Ahn, S. Park, and B.-J. Yang, Phys. Rev. X 9, 021013 (2019); H. C. Po, L. Zou, T. Senthil, and A. Vishwanath, Phys. Rev. B 99, 195455 (2019); J. Kang and O. Vafek, Phys. Rev. X 8, 031088 (2018); M. Koshino, N. F. Q. Yuan, T. Koretsune, M. Ochi, K. Kuroki, and L. Fu, Phys. Rev. X 8, 031087 (2018); J. Liu, J. Liu, and X. Dai, Phys. Rev. B 99, 155415 (2019); L. Zou, H. C. Po, A. Vishwanath, and T. Senthil, Phys. Rev. B 98, 085435 (2018)] found that the two flatbands in TBG possess a fragile topology protected by the C 2z T symmetry. Song et al. [Z. Song, Z. Wang, W. Shi, G. Li, C. Fang, and B. A. Bernevig, Phys. Rev. Lett. 123, 036401 (2019)] also pointed out an approximate particle-hole symmetry (P) in the continuous model of TBG. In this paper, we numerically confirm that P is indeed a good approximation for TBG and show that the fragile topology of the two flatbands is enhanced to a P-protected stable topology. This stable topology implies 4l + 2 (l ϵ N) Dirac points between the middle two bands. The P-protected stable topology is robust against arbitrary gap closings between the middle two bands and the other bands. We also show that, remarkably, this P-protected stable topology, as well as the corresponding 4l + 2 Dirac points, cannot be realized in lattice models that preserve both C 2z T and P symmetries. In other words, the continuous model of TBG is anomalous and cannot be realized on lattices. Two other topology related topics, with consequences for the interacting TBG problem, i.e., the choice of Chern band basis in the two flatbands and the perfect metal phase of TBG in the so-called second chiral limit, are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Adaptive variational quantum eigensolvers for highly excited states

Highly excited states of quantum many-body systems are central objects in the study of quantum dynamics and thermalization that challenge classical computational methods due to their volume-law entanglement content. In this work, we explore the potential of variational quantum algorithms to approximate such states. We propose an adaptive variational quantum eigensolver (VQE) for excited states (X) that self-generates a variational ansatz for arbitrary eigenstates of a many-body Hamiltonian H by attempting to minimize the energy variance with respect to H. We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. Here, we also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Analytical and EZmock covariance validation for the DESI 2024 results

The estimation of uncertainties in cosmological parameters is an important challenge in Large-Scale-Structure (LSS) analyses. For standard analyses such as Baryon Acoustic Oscillations (BAO) and Full-Shape two approaches are usually considered. First: analytical estimates of the covariance matrix use Gaussian approximations and (nonlinear) clustering measurements to estimate the matrix, which allows a relatively fast and computationally cheap way to generate matrices that adapt to an arbitrary clustering measurement. On the other hand, sample covariances are an empirical estimate of the matrix based on an ensemble of clustering measurements from fast and approximate simulations. While more computationally expensive due to the large amount of simulations and volume required, these allow us to take into account systematics that are impossible to model analytically. In this work we compare these two approaches in order to enable DESI's key analyses. We find that the configuration space analytical estimate performs satisfactorily in BAO analyses and its flexibility in terms of input clustering makes it the fiducial choice for DESI's 2024 BAO analysis. On the contrary, the analytical computation of the covariance matrix in Fourier space does not reproduce the expected measurements in terms of Full-Shape analyses, which motivates the use of a corrected mock covariance for DESI's 2024 Full Shape analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Calorimetric Energies and Chemical Potentials of Metal Atoms in Catalytic Nanoparticles on Oxide and Carbon Supports: Improved Size Dependencies and Adhesion Energies

The chemical potential of metal atoms, μM, in supported metal nanoparticles is an important descriptor related to both the catalytic activity and the stability of the nanoparticles. Here, we derive an expression relating μM to the radius of the particle’s contact area with the support and the adhesion energy at the metal/support interface (E adh ) that assumes the particles have the shape of spherical caps but of arbitrary contact angle with the support (θ c ) and includes an empirical correction for the increase in metal surface energy and adhesion energy with decreasing radius of curvature. We then show that, at any assumed contact angle, we can simultaneously fit previously reported measurements of both calorimetric μ M (from heats of metal vapor adsorption during nanoparticle growth by vapor deposition) versus metal coverage data and the He + low-energy ion scattering (LEIS) intensities for the metal and/or support versus metal coverage (using our recently developed spherical cap model for quantitative LEIS intensities), to determine the particle size versus coverage and E adh . Only one choice of contact angle gives a pair of values for contact angle and E adh , which is consistent with the Young–Dupré equation for the equilibrium shape of a spherical particle. At this equilibrium shape, we then applied this spherical cap model (SCM) to reanalyze microcalorimetric metal chemical potentials and LEIS signals versus coverage data for nine metal/support combinations that were previously analyzed by assuming that the particles had the shape of hemispherical caps, i.e., with a contact angle of 90°. We show that this revised approach gives close agreement with the calorimetric and LEIS data; the best-fit contact angles vary from 64 to 84°, correcting the earlier assumption of 90°. Furthermore, these results provide significant accuracy improvements in particle size versus coverage, metal chemical potential versus size and coverage, metal/support adhesion energies and contact angles for Cu, Ag and Au on CeO 2 (111), Ni on MgO(100), Ag on Fe 3 O 4 (111) and TiO 2 (100), and Ag, Ni and Pd on Ni-supported graphene. This revised approach is much more broadly applicable than the earlier hemispherical cap model (HCM).

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Elimination of QCD Renormalization Scale and Scheme Ambiguities

The setting of the renormalization scale (μ r ) in the perturbative QCD (pQCD) is one of the crucial problems for achieving precise fixed-order pQCD predictions. The conventional prescription is to take its value as the typical momentum transfer Q in a given process, and theoretical uncertainties are then evaluated by varying it over an arbitrary range. The conventional scale-setting procedure introduces arbitrary scheme-and-scale ambiguities in fixed-order pQCD predictions. The principle of maximum conformality (PMC) provides a systematic way to eliminate the renormalization scheme-and-scale ambiguities. The PMC method has rigorous theoretical foundations; it satisfies the renormalization group invariance (RGI) and all of the self-consistency conditions derived from the renormalization group. The PMC has now been successfully applied to many physical processes. In this paper, we summarize recent PMC applications, including event shape observables and heavy quark pair production near the threshold region in e + e – annihilation and top-quark decay at hadronic colliders. In addition, estimating the contributions related to the uncalculated higher-order terms is also summarized. These results show that the major theoretical uncertainties caused by different choices of μ r are eliminated, and the improved pQCD predictions are thus obtained, demonstrating the generality and applicability of the PMC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Grating Polarizers at 170 GHz for ECRH Systems: Low Power Tests and Simulations

The elliptical polarization produced by a pair of sinusoidally profiled miter-bend grating polarizers has been measured experimentally and found to be in excellent agreement with theory. The polarizers were designed to be used at 170 GHz in the miter bends of a 63.5 mm corrugated metallic waveguide for creating the arbitrary elliptically polarized microwave beam needed for electron cyclotron resonance heating of plasma. Using a vector network analyzer to generate a linearly polarized HE11 incident mode, each polarizer was individually tested to measure the rotation α and ellipticity β of the elliptically polarized reflected microwave beam as a function of the grating rotation angle. The grating polarizers were then tested together to measure the elliptical polarization as a function of the angles of the combination of the two polarizers arranged sequentially on the transmission line. The map of the ellipticity versus grating angles agreed very well with numerical simulations using high-frequency structure simulator. Numerical simulations show that up to eight combinations of rotation angles of the two grating polarizers can provide the same ellipticity but with the varying ohmic loss, thus allowing the choice of settings to minimize the ohmic loss, a necessary consideration to minimize heating effects in megawatt power level systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Noise-aware circuit compilations for a continuously parameterized two-qubit gateset

State-of-the-art noisy-intermediate-scale quantum processors are currently implemented across a variety of hardware platforms, each with their own distinct gatesets. As such, circuit compilation should not only be aware of but also deeply connect to the native gateset and noise properties of each. Trapped-ion processors are one such platform that provides a gateset that can be continuously parameterized across both one- and two-qubit gates. Here we use the Quantum Scientific Computing Open User Testbed to study noise-aware compilations focused on continuously parameterized two-qubit 𝑍⁢𝑍 gates (based on the Mølmer-Sørensen interaction) using $\scriptsize{SUPERSTAQ}$, a quantum software platform for hardware-aware circuit compiler optimizations. We discuss the realization of 𝑍⁢𝑍 gates with arbitrary angle on the all-to-all connected trapped-ion system. Then we discuss a variety of different compiler optimizations that innately target these 𝑍⁢𝑍 gates and their noise properties. These optimizations include moving from a restricted maximally entangling gateset to a continuously parameterized one, swap mirroring to further reduce the total entangling angle of the operations, focusing the heaviest 𝑍⁢𝑍 angle participation on the best-performing gate pairs, and circuit approximation to remove the least impactful 𝑍⁢𝑍 gates. We demonstrate these compilation approaches on the hardware with randomized quantum volume circuits, observing the potential to realize a larger quantum volume as a result of these optimizations. Using differing yet complementary analysis techniques, we observe the distinct improvements in system performance provided by these noise-aware compilations and study the role of stochastic and coherent error channels for each compilation choice.

Noise↗

DOE OSTI2023

Coupled magnetostructural continuum model for multiferroic BiFeO 3

Here, a continuum approach to study magnetoelectric multiferroic BiFeO 3 (BFO) is proposed. Our modeling effort marries the ferroelectric (FE) phase field method and micromagnetic simulations to describe the entire multiferroic order parameter sector (polarization, oxygen antiphase tilts, strain, and magnetism) self-consistently on the same time- and length scale. In this paper, we discuss our choice of FE and magnetic energy terms and demonstrate benchmarks against known behavior. We parametrize the lowest order couplings of the structural distortions against previous predictions from density functional theory calculations giving access to simulations of the FE domain wall (DW) topology. This allows us to estimate the energetic hierarchy and thicknesses of the numerous structural DWs. We then extend the model to the canted antiferromagnetic order and demonstrate how the FE domain boundaries influence the resulting magnetic DWs. We also highlight some capabilities of this model by providing two examples relevant for applications. We demonstrate spin-wave transmission through the multiferroic domain boundaries which identify rectification in qualitative agreement with recent experimental observations. As a second example of application, we model fully dynamical magnetoelectric switching, where we find a sensitivity on the Gilbert damping with respect to switching pathways. We envision that this modeling effort will set the basis for further work on properties of arbitrary 3D nanostructures of BFO (and related multiferroics) at the mesoscale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

IRMA

IRMA (In)elastic Representation of Materials As S(α,β) evaluations IRMA turns one phonon model into three outputs that usually require three separate tool chains: an evaluated nuclear-data file, predicted neutron-scattering spectra, and scattering kernels for Monte Carlo transport. The three outputs draw on a single, consistent description of the material, so the evaluation, the spectroscopy that can validate it, and the transport that uses it always agree about the physics. Nuclear data. IRMA writes ENDF-6 File 7 thermal scattering evaluations on automatically constructed (α, β) grids. This part reimplements and generalizes NJOY's LEAPR: the classic kernels reproduce freshly generated NJOY2016 tapes digit for digit and published reference tapes to about 1e-4, and the generalized paths add the exact coherent one-phonon term, anisotropic Debye-Waller tensors, coherent elastic for arbitrary crystals, and a per-species partition for polyatomic materials. The tapes feed NJOY, AMPX, FUDGE, and every transport code downstream of them. Neutron spectroscopy. The irma.spectra forward model projects the same physics onto an instrument's kinematics and resolution: INS spectra for VISION and generic indirect geometries, and 2-D S(Q,E) powder maps for direct-geometry spectrometers, from a phonopy model or straight from a phonon DOS. It can be used to predict a proposed measurement before beam time; in analysis, it supplies the calculated single-scattering counterpart of a measured spectrum, from the same material description the evaluation was built from. Monte Carlo transport. The irma.ncrystal exporter writes per-temperature scattering kernels for the companion NCrystal plugin, so McStas, OpenMC, and other NCrystal-aware codes sample the same physics. The exported kernels carry the per-site anisotropic Debye-Waller tensors, keeping directional coherent-elastic physics that NCrystal's standard scalar treatment does not represent. With the same physics inside a transport code, an entire beamline becomes a virtual experiment: IRMA's end-to-end validation ran a custom McStas implementation of the ARCS spectrometer, assembled from the existing McVine and McStas models, against measured data. From a bare crystal structure. The irma mlip front end builds the phonon model itself: a structure file and a choice of potential are enough. Nine pretrained machine-learned interatomic potentials are supported, on a laptop CPU, with no first-principles calculation; an approximate phonon model for a new material costs minutes, not a DFT campaign, and the build emits prefilled inputs for all three outputs. The result is a good starting point rather than a finished evaluation: survey-quality physics with every parameter exposed for review. A converged atomistic calculation enters the same way, as a phonopy model, when higher fidelity is needed.

Ramic, Kemal [Oak Ridge National Laboratory (ORNL)↗

DOE OSTI2023

Generalized echo squeezing protocol with near-Heisenberg-limit sensitivity and strong robustness against detection noise and variation in squeezing parameter

We present a generalized echo squeezing protocol (GESP) as a generalization of the Schrödinger cat state protocol with the value of the squeezing parameter being an arbitrary number rather than π / 2 . We show analytically that over a very broad range of the squeezing parameter the sensitivity of the GESP reaches the Heisenberg limit within a factor of 2 . For a large number of atoms, N , this plateau interval is almost the whole range from zero to π / 2 , and the sensitivity is independent of the parity of N . Therefore, it is possible to operate a sensor over a wide interval of the squeezing parameter without changing the sensitivity. This is to be contrasted with the conventional echo squeezing protocol (CESP) which only works for a very small interval of the squeezing parameter. We also show that, in contrast to the CESP, the sensitivity of the GESP is close to the quantum Cramér-Rao bound over the whole range of the squeezing parameter, indicating that the phase-shift information contained in the quantum state is near-optimally extracted. We find that the enhancement in sensitivity in the case of the GESP is due to a combination of two parameters: the phase magnification factor (PMF) and the noise amplification factor (NAF). As the value of the squeezing parameter increases, both PMF and NAF increase, while keeping the ratio of the PMF to the NAF essentially constant, yielding a net enhancement of sensitivity at the Heisenberg limit within a factor of 2 over the whole plateau interval. An important consequence of this behavior is that the robustness of the GESP against detection noise easily exceeds that of the CESP for a broad range of values of the squeezing parameter. As such, in the context of an experimental study, it should be possible to achieve a net enhancement in sensitivity higher than that for the CESP, under typical conditions where detection noise exceeds the quantum projection noise of an unsqueezed state with the same number of atoms. Finally, we consider the fragility of the GESP against decoherence mechanisms, and show how a balance between the fragility against the decoherence mechanisms and the robustness against detection noise would in practice determine the optimal choice of parameters for the GESP.

Li, Jinyang↗