SEARCH · Engineering Papers

Results for “Matrix elements”

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.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2023

The Neutron Mean Life and Big Bang Nucleosynthesis

We explore the effect of neutron lifetime and its uncertainty on standard big bang nucleosynthesis (BBN). BBN describes the cosmic production of the light nuclides, 1H, D, 3H+3He, 4He, and 7Li+7Be, in the first minutes of cosmic time. The neutron mean life τn has two roles in modern BBN calculations: (1) it normalizes the matrix element for weak n↔p interconversions, and (2) it sets the rate of free neutron decay after the weak interactions freeze-out. We review the history of the interplay between τn measurements and BBN, and present a study of the sensitivity of the light element abundances to the modern neutron lifetime measurements. We find that τn uncertainties dominate the predicted 4He error budget, but these theory errors remain smaller than the uncertainties in 4He observations, even with the dispersion in recent neutron lifetime measurements. For the other light element predictions, τn contributes negligibly to their error budget. Turning the problem around, we combine present BBN and cosmic microwave background (CMB) determinations of the cosmic baryon density to predict a “cosmologically preferred” mean life of τn(BBN+CMB)=870±16s, which is consistent with experimental mean life determinations. We show that if future astronomical and cosmological helium observations can reach an uncertainty of σobs(Yp)=0.001 in the 4He mass fraction Yp, this could begin to discriminate between the mean life determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Strong Interaction Dynamics and Fermi $β$ Decay in the Nucleon and the Nucleus

Nuclear super-allowed β decay has been used to obtain tight limits on the value of the CKM matrix element V ud that is important for unitarity tests and, therefore, for tests of the standard model. Current requirements on precision are so intense that effects formerly thought too small to matter are now considered relevant. This article is a brief review of personal efforts to include the effects of strong interactions on Fermi β decay. First, I examine the role of isospin violation in the decay of the neutron. The size of the necessary correction depends upon detailed strong-interaction dynamics. The isospin violating parts of the nucleon wave function, important at the low energy of β decay, can be constrained by data taken at much higher energies, via measurements, for example, of ed → e'π ± + X reactions at Jefferson Laboratory. The next point of focus is on the role of nuclear short-ranged correlations, which affect the value of the correction needed to account for isospin violation in extracting the value of V ud . The net result is that effects previously considered as irrelevant are now considered relevant for both neutron and nuclear β decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Improving the Precision of First-Principles Calculation of Parton Physics from Lattice Quantum Chromodynamics

Large momentum effective theory (LaMET) provides a general framework for computing the multi-dimensional partonic structure of the proton from first principles using lattice quantum chromodynamics (QCD). In this effective field theory approach, LaMET predicts parton distributions through a power expansion and perturbative matching of a class of Euclidean observables—quasi-distributions—evaluated at large proton momenta. Recent advances in lattice renormalization, such as the hybrid scheme with leading renormalon resummation, together with improved matching kernel that incorporates higher-loop corrections and resummations, have enhanced both the perturbative and power accuracy of LaMET, enabling a reliable quantification of theoretical uncertainties. Moreover, the Coulomb-gauge correlator approach further simplifies lattice analyses and improves the precision of transverse-momentum-dependent structures, particularly in the non-perturbative region. State-of-the-art LaMET calculations have already yielded certain parton observables with important phenomenological impact. In addition, the recently proposed kinematically enhanced lattice interpolation operators promise access to unprecedented proton momenta with greatly improved signal-to-noise ratios, which will extend the range of LaMET prediction and further suppress the power corrections. The remaining challenges, such as controlling excited-state contamination in lattice matrix elements and extracting gluonic distributions, are expected to benefit from emerging lattice techniques for ground-state isolation and noise reduction. Thus, lattice QCD studies of parton physics have entered an exciting stage of precision control and systematic improvement, which will have a broader impact for nuclear and particle experiments.

Zhao, Yong [Argonne National Laboratory (ANL), Arg↗

DOE OSTI2020

An Investigation of Poorly Studied Open Cluster NGC 4337 Using Multicolor Photometric and Gaia DR2 Astrometric Data

We present a comprehensive analysis (photometric and kinematical) of the poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia DR2 databases. By determining the membership probabilities of stars, we identified the 624 most probable members with membership probability higher than 50% by using proper motion and parallax data taken from Gaia DR2. The mean proper motion of the cluster is obtained as μ{sub x}=−8.83±0.01 and μ {sub y} = 1.49 ± 0.006 mas yr{sup −1}. We find the normal interstellar extinction toward the cluster region. The radial distribution of members provides a cluster radius of 7.′75 (5.63 pc). The estimated age of 1600 ± 180 Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as 1.46±0.18 within the mass range 0.75−2.0 M{sub ⊙}, which is in fair agreement with Salpeter’s value (x = 1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates (A, D) by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are (A {sub conv}, D {sub conv}) = (96.°27 ± 0.°10, 13.°14 ± 0.°27) and (A {sub ◦}, D {sub ◦}) = (100.°282 ± 0.°10, 9.°577 ± 0.°323) respectively. We also computed the Velocity Ellipsoid Parameters, matrix elements (μ {sub ij}), direction cosines (l {sub j}, m {sub j}, n {sub j}), and the Galactic longitude of the vertex (l {sub 2}).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Reexamining Circular Dichroism in Photoemission From a Topological Insulator

The orbital angular momentum (OAM) of electron states is an essential ingredient for topological and quantum geometric quantities in solids. For example, Dirac surface states with helical spin- and orbital-angular momenta are a hallmark of a 3D topological insulator. Angle-resolved photoemission spectroscopy (ARPES) with variable circular light polarization, known as circular dichroism (CD), has been assumed to be a direct probe of OAM and, by proxy, of the Berry curvature of electronic bands in energy- and momentum-space. Indeed, topological surface states have been shown to exhibit angle-dependent CD (CDAD), and more broadly, CD is often interpreted as evidence of spin-orbit coupling. Meanwhile, it is well-established that CD originates from the photoemission matrix elements, which can have extrinsic contributions related to the experimental geometry and the inherently broken inversion symmetry at the sample surface. Therefore, it is important to broadly examine CD-ARPES to determine the scenarios in which it provides a robust probe of intrinsic material physics. We performed CD-ARPES on the canonical topological insulator $\mathrm{Bi}_2\mathrm{Se}_3$ over a wide range of incident photon energies. Not only do we observe angle-dependent CD in the surface states, as expected, but we also find CD of a similar magnitude in virtually all bulk bands. Since OAM is forbidden by inversion symmetry in the bulk, we conclude this originates from symmetry-breaking in the photoemission process. Comparison with theoretical calculations supports this view and suggests that $\textit{hidden}$ OAM - localized to atomic sites within each unit cell - contributes significantly. Additional effects, including inter-atomic interference and final-state resonances, are responsible for the rapid variation of the CDAD signal with photon energy.

FOS: Physical sciences↗

DOE OSTI2025

Gluonic Structure from Instantons

The instanton vacuum provides an effective description of chiral symmetry breaking by local topological fluctuations of the gauge fields, as observed in lattice QCD simulations. The resulting effective dynamics at momenta below 1/$\bar{ρ}$ ≈ 0.6 GeV explains the basic features of light-quark correlation functions and is used extensively in studies of hadron structure. The instanton fields also make definite contributions to the gluonic structure of light hadrons, as expressed in the matrix elements of composite quark– gluon or gluon operators. The article reviews the gluonic structure of light hadrons (nucleon, pion) induced by instantons. This includes: (i) twist-2 parton distributions and momentum sum rule; (ii) twist-3 angular momentum and spin-orbit interactions; (iii) twist-3 and twist-4 quark–gluon correlations and power corrections; (iv) trace anomaly and hadron mass decomposition; (v) scalar gluon form factors and mechanical properties; (vi) axial anomaly and pseudoscalar gluon form factors. It also discusses possible further applications of the methods and recent developments including gauge field configurations beyond instantons.

Weiss, C. [Thomas Jefferson National Accelerator F↗

DOE OSTI2024

Structure of Low-spin States in $^{45}$Sc Studied via Coulomb Excitation

The electromagnetic structure of 45 Sc at low excitation energy was investigated via low-energy Coulomb excitation at the Heavy Ion Laboratory (HIL) of the University of Warsaw and at the Inter-University Accelerator Centre (IUAC) in New Delhi. A set of reduced E2, E3, and M1 matrix elements was extracted from the collected data using the GOSIA code. The reduced transition probability B(E2; 11/2 – → 7/2 – ) has been determined, allowing us to deduce the lifetime of the 11/2 – state at 1237 keV. In addition, the upper limit on the reduced transition probability B(E3; 7/2 – → 5/2 + ) has been determined for the first time. Finally, new large-scale shellmodel and beyond-mean-field calculations were performed to interpret the structure of this nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Real time evolution for ultracompact Hamiltonian eigenstates on quantum hardware

In this work we present a detailed analysis of variational quantum phase estimation (VQPE), a method based on real-time evolution for ground and excited state estimation on near-term hardware. We derive the theoretical ground on which the approach stands, and demonstrate that it provides one of the most compact variational expansions to date for solving strongly correlated Hamiltonians. At the center of VQPE lies a set of equations, with a simple geometrical interpretation, which provides conditions for the time evolution grid in order to decouple eigenstates out of the set of time evolved expansion states, and connects the method to the classical filter diagonalization algorithm. Further, we introduce what we call the unitary formulation of VQPE, in which the number of matrix elements that need to be measured scales linearly with the number of expansion states, and we provide an analysis of the effects of noise which substantially improves previous considerations. The unitary formulation allows for a direct comparison to iterative phase estimation. Our results mark VQPE as both a natural and highly efficient quantum algorithm for ground and excited state calculations of general many-body systems. We demonstrate a hardware implementation of VQPE for the transverse field Ising model. Further, we illustrate its power on a paradigmatic example of strong correlation (Cr2 in the SVP basis set), and show that it is possible to reach chemical accuracy with as few as ~50 timesteps.

Klymko, Katherine↗

DOE OSTI2022

Neumann Series in MGS-GMRES and Inner-Outer Iterations: Preprint

A low-synchronization MGS-GMRES Krylov solver employing a truncated Neumann series for the inverse compact WY MGS correction matrix T is presented. A corollary to the backward stability result of Paige et al. [1] establishes that T = I - Lk is sufficient for convergence of GMRES when kLkp F = O("p)_p F (B), where the strictly lower triangular matrix L is defined by the inner products of Krylov vectors V T 1:k-2 vk-1. The preconditioner is the classical Ruge-Stuben AMG algorithm with compatible relaxation and inner-outer Gauss-Seidel smoother. This smoother may also be expressed as a truncated Neumann series. Drop tolerances are applied to the lower triangular matrices arising in the smoother in order to reduce the number of non-zeros and accelerate the time to solution. The number of small matrix elements are found to increase from fine to coarse levels and thus the effciency gains are greater for large problems with many levels in the V -cycle. The solver is applied to the pressure continuity equation for the incompressible Navier-Stokes equations. Unlike the inner-outer iteration, the solver convergence rate with the standard Gauss-Seidel smoother deteriorates with dropping. The solver compute time is reduced by up to 50% without a change in the convergence rate.

Gauss-Seidel smoother↗

DOE OSTI2022

Assessing the theory-data tension in neutrino-induced charged pion production: the effect of final-state nucleon distortion

Pion production on nuclei constitutes a significant part of the total cross section in experiments involving few-GeV neutrinos. Combined analyses of data on deuterium and heavier nuclei points to tensions between the bubble chamber data and the data of the MINER$\nu$A experiment, which are often ascribed to unspecified nuclear effects. To understand the origin of these tensions, a microscopic quantum mechanical framework is needed to compute nuclear matrix elements. We use the local approximation to the relativistic distorted wave impulse approximation (RDWIA) to assess the role of final-state nucleon distortion. To perform this comparison under conditions relevant to neutrino experiments, we compute cross sections for the MINER$\nu$A and T2K charged pion production datasets. The inclusion of nucleon distortion leads to a reduction of the cross section up to 10%, but to no significant change in shape of the flux-averaged cross sections. Results with and without distortion compare favorably to experimental data, with the exception of the low-$Q^2$ MINER$\nu$A $\pi^+$ data. We point out that hydrogen target data from BEBC is also overpredicted at low-$Q^2$, and that the discrepancy is similar in shape and magnitude to what is found in comparison to MINER$\nu$A data. Including nucleon distortion alone cannot explain the overprediction of low-$Q^2$ cross sections measured by MINER$\nu$A. The similar overprediction of BEBC data on hydrogen means that it is impossible to ascribe this discrepancy solely to a nuclear effect. Axial couplings and their $Q^2$ dependence should ideally be derived from more precise data on hydrogen and deuterium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Hadron spectroscopy of pseudoscalar and vector meson photoproduction using linearly polarised photons at CLAS

The field of hadron spectroscopy is informed by experiments which study excited states of the nucleon by detecting the decay products following excitation with a hadronic, leptonic or electromagnetic probe. Experiments using a polarised beam of photons provide an important contribution to the world data. This thesis describes the analysis of data from one such experiment, the CEBAF Large Acceptance Spectrometer at Jefferson Laboratory, Virginia, USA. In the experiment, an incident beam of electrons passes through a diamond radiator producing a linearly polarised photon beam of energy from 1.1 to 2.1 GeV on a liquid hydrogen target. The analysis method and results for two reaction channels are presented. In the first analysis, the polarisation observables {?, P, T, Ox, Oz} for the reaction ??? p ? K0 S ?+ are extracted simultaneously using likelihood sampling with Markov-Chain Monte-Carlo. The results are presented as angular distributions (cos ?K0S) for bins in photon energy E? for a total of 21 bins. The values extracted for T , Ox and Oz are a first measurement of these quantities for this reaction and the beam asymmetry (?) measurement supplements the one previous measurement for this quantity and extends the energy range over which it has been extracted. The measurement of the recoil polarisation, P, is consistent with previous measurements. In a second analysis, the spin density matrix elements (SDMEs) for ? p ? p ? are extracted using similar analysis techniques. The results are presented as angular distributions (cos ?? ) for bins in photon energy E? for a total of 24 bins. This work is the first to extract all nine SDMEs for the full angular range of the ? meson production. A comparison is made to the predictions of the Vector Meson Dominance (VMD) model and the results provide evidence for contributing processes other than VMD, particularly at more backward angles of the ? meson production. The data also indicate non-helicity conserving process in both the s-channel and the t-channel. The results from the first analysis have been passed to the Juelich-Bonn theory group and the effect of including the new data within their dynamical coupled channel model is described. Several nucleon resonances are affected and the impact on the pole positions, widths and photon decay amplitudes of the affected resonances is described.

Clark, Louise↗

DOE OSTI2024

Photoproduction of K0K0-bar with the GlueX Experiment

In this dissertation, we study photoproduction of K0 ?K 0 through the KSKLp and KSKSp ?nal states with the GlueX Phase-I (GlueX-I) data set. We measure the spin-density matrix elements (SDMEs) and di?erential cross section of ?(1020) ! KSKL to better understand photoproduction of light vector mesons. A mass independent Partial Wave Analysis (PWA) of the KSKL and KSKS meson spectrum below 2 GeV is also undertaken to study the meson spectrum with JPC = even++ and odd??. The ?(1020) di?erential cross section is measured at E = 8:2 ? 8:8 GeV and ?t = 0:15 ? 1:00 GeV2. The di?erential cross section is well described by an exponential decay with slope 4:44?0:01 GeV?2 and the integrated cross section is determined to be 295:7?0:4 nb, only statistical uncertainties are quoted. Both measurements are consistent and far more precise than the previous measurement by Ballam et al. [1]. The ?(1020) SDMEs were measured in nine bins of ?t in the same range. At low ?t, we ?nd the data were consistent with s-channel helicity conservation, SCHC, i. e. the only non-zero SDMEs were ?1 1?1 = ?Im(?2 1?1) = 1=2. At higher ?t, the SDMEs deviate from SCHC and the measurements indicate this is due to natural parity exchange since we observe that the SDMEs are consistent with zero unnatural exchange. We also ?find that the contribution from helicity double-flip amplitudes is consistent with zero. The measured SDMEs are in poor agreement with theoretical predictions put forward by JPAC [35]. These measurements will serve as input to re?ne models of production processes, which will be essential for the interpretation of possible signals of exotic mesons in GlueX. Analysis of the spectrum above the ?(1020) indicates the presence of at least three particles at ? 1:50, ? 1:75, and 2:20 GeV. We employ simple parametrizations of the resonance line shape based on relativistic Breit-Wigner functions considering cases with and without interference. The resonance at ? 1:50 GeV is not identifi?ed as any specifi?c resonance, but could be due to interference between the ?(1450) and !(1420). Our models favor the resonance at ~1:75 GeV as the X(1750) rather than the ?(1680). We ?find that introducing a third Breit-Winger into the models improves the ?t quality with the ?2/ndf going from 1.75 to 1.59 and 1.41 to 1.18 for models without and with interference respectively. A PWA of the KSKL spectrum indicates that the spectrum is consistent with exclusively spin-1 contribution up to ? 1:6 GeV. Above ? 1:6 GeV we ?find evidence for a small spin-3 contribution. We also ?nd that the spectrum predominantly positive reflectivity up to ? 1:6 GeV, after which the spectrum becomes a nearly equal mix of both. The KSKS system shows a rich spectrum with multiple resonance structures but is more statistically limited than the KSKL system. Modelling the line shape we ?find evidence at greater than 10? level in favor of a resonance at ? 1:75 GeV. The model parameters suggest this resonance is the f0(1710). The PWA of the KSKS system was inconclusive. However, using a small set of amplitudes suggests that the spin-2 contribution to the spectrum is primarily in the range 1:2?1:6 GeV, where f2(1270), a2(1320) and f02 (1525) are expected.

Linera, Gabriel Rodriguez↗

DOE OSTI2024

Gravitational form factors of glueballs in Yang-Mills theory

This work presents preliminary results of the first determination of the energy-momentum tensor form factors of the scalar glueball, referred to as gravitational form factors (GFFs). The calculation has been carried out in lattice Yang-Mills theory at a single lattice spacing. Using variationally optimized operators, the matrix elements are extracted from ratios of three-point functions to two-point functions. The glueball GFFs and their kinematic dependence are compared to those of other hadrons from previous calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

A Thermal Relic Encyclopedia: Dark Matter Candidates Coupled to Quarks

Thermal freeze-out is a compelling framework for naturally generating the dark matter abundance. We systematically study a broad range of dark matter and mediator particle combinations that can viably realize thermal freeze-out, focusing on models in which the mediator couples to Standard Model quarks. In each case, we calculate the relic density and consider existing constraints from accelerators, cosmology, direct detection, and indirect detection over the full range of dark matter and mediator masses. We present an encyclopedic catalog of matrix elements, cross sections, and decay rates which can be used as a reference for future studies of dark matter phenomenology.

Hooper, Dan [Fermilab; Chicago U., EFI; Chicago U.↗

DOE OSTI2026

Extracting $B_s\to D_s^*\ellν_\ell$ form factors

Semileptonic $B_{(s)}$ decays are of great phenomenological interest because they allow to determine e.g. CKM matrix elements or test lepton flavor universality. Taking advantage of already existing lattice data, we demonstrate the analysis steps to extract the four form factors describing exclusive semileptonic $B_s\to D_s^*\ellν_\ell$ decays using the narrow width approximation. Our data are based on RBC/UKQCD's set of 2+1 flavor gauge field ensembles with Shamir domain-wall fermion and Iwasaki gauge field action featuring inverse lattice spacings of $a^{-1}=$1.785, 2.383, and 2.785 GeV as well as pion masses between 268 and 433 MeV. Light, strange and charm quarks are simulated using domain-wall fermions, whereas bottom quarks are generated with the relativistic heavy quark (RHQ) action.

Boushmelev, Anastasia [Siegen U.] (ORCID:000000033↗

DOE OSTI2025

Internal-Gelation Production of Uranium Oxide Sol-Gel Particles for Forensic Applications

The purpose of this project was to develop and demonstrate a novel method for the production of uranium oxide microsphere particles with tunable chemical compositions via a sol-gel process using a 3D-printer setup. These particles can serve several purposes in research and development as a forensic training tool or as standard reference materials. A key component of the project was to demonstrate the ability to control physical and chemical parameters of the particles created. First, we demonstrated the ability to employ an internal gelation sol-gel process to create individual uranium oxide particles. The particles were successfully dispensed using a unique 3D-printing setup onto a substrate to react and then were collected and thermally processed. A series of temperatures for the annealing process was tested on individual samples to investigate the effect on the sol-gel chemical composition and physical integrity. Next, we demonstrated the ability to control matrix composition of the particles by separately incorporating fission product isotopes as well as Np-237 into the sol-gel solution. It was shown by gamma-ray spectroscopy that these matrix elements were successfully retained during the gelation process. We studied the retention of the elements across a series of annealing temperatures. Additionally, we demonstrated the ability to quantitatively control the isotopic composition of the particles by altering the U-237/U-238 ratio to a controlled value. Finally, X-ray diffraction analysis (XRD) was used to investigate the oxidation state of the sol-gel after annealing at different temperatures.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

In this work, we provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint-based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. Time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

97 MATHEMATICS AND COMPUTING↗