Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conjugate spaces”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quinonoid versus Aromatic π-Conjugated Oligomers and Polymers and Their Diradical Characters

Quinonoid elements in π-conjugated oligomers and polymers have promising characteristics by providing unique structures and electronic properties, including varying degrees of ground state open shell diradical character and small band gaps. So far, quinonoid ground state polymers were rarely identified, and one of the goals of this work was to expand the portfolio of such systems in which over 17 different repeat units were considered, about half preferred a quinonoid and half preferred aromatic structures plus several with the alternating A-Q composition were also included. The dependency of the diradical character as a function of size and composition opens a rich space to discover trends affecting the aromaticity and diradical character. Herein, we show that the diradical character is a useful parameter to classify a wide variety of polymers based on their ground states in a systematic study of a large number of homo- and heterosystems. Homo-polymers with quinonoid ground states show a high diradical character, while those with the aromatic ground state behave as closed-shell systems with no or a very small diradical character. The diradical character of quinonoid systems depends strongly on the size, which increases with the size of the oligomer. Here, we found a correlation between the diradical character and the interring C–C distance as well as the singlet–triplet energy gaps. In alternating copolymers mixing quinonoid and aromatic subunits, we found open-shell hetero-polymers with moderate to high diradical characters.

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

Specific heats for quantum BTZ black holes in extended thermodynamics

It was shown recently that extended black hole thermodynamics, where the cosmological constant is a dynamical variable, giving rise to a pressure p and its conjugate volume V , can be given a natural setting in the context of braneworld models. We study the specific heat capacities C p ( T ) and C V ( T ) of the quantum version of the Bañados, Teitelboim, and Zanelli (BTZ) black hole that lives in the induced gravity theory on the brane. There are multiple branches of solutions, and we explore and characterize key features of the possible behavior. We identify and study a critical point in the space of solutions where both specific heats diverge. In the regime of weak backreaction where we are close to an ordinary theory of gravity, the black hole is “subentropic,” but as backreaction is increased we note that there are parts of parameter space that has regions where it is “superentropic.” While a study of the sign of the specific heats does not always show a corresponding instability (conjectured in the literature), the presence of strong backreaction makes interpretation unclear. Published by the American Physical Society 2024

Johnson, Clifford V. (ORCID:0000000189645830)↗
DOE OSTI2026

A Measurement of the Pion-Energy Dependence of Muon Neutrino Charged-Current Scattering to Final States With One Charged Pion in NOvA

The study of neutrino oscillations is a main priority for particle physics as the most immediately tractable lever on physics beyond the Standard Model. In particular, more insight into violation of the combined symmetry of charge-conjugation plus parity could yield clues to the origin of matter-antimatter asymmetry, and in some theoretical frameworks the neutrino mass could give insights into dark matter. NOvA is a long-baseline accelerator neutrino experiment with both a near and far detector that seeks to measure several of the parameters of the neutrino mixing matrix, as well as carry out a broad program of additional physics. NOvA has innovated and developed a variety of techniques in the space of neutrino physics, including expanding the use of machine learning techniques in reconstruction. NOvA has also set the stage for the next generation US-based long baseline experiment, DUNE. This dissertation details the creation of a new neutrino interaction vertex reconstruction package for NOvA, which offers enormous improvements in accuracy above the previous vertexer it replaces. This is accomplished using a Convolutional Visual Network trained on large datasets of simulated events. The vertexer is validated thoroughly against data. This dissertation also presents a cross-section measurement for $\nu_\mu + N \rightarrow \mu^- + 1\pi^\pm + X$ (where X does not include additional charged pions) binned in pion kinetic energy, a challenging measurement in general and for NOvA in particular. Measurements like this provide crucial inputs for neutrino interaction models utilized by neutrino oscillation experiments. This measurement in particular is important to the energy ranges relevant for DUNE in a band not well-covered by any other experiment.

Ewart, Erin [Indiana U., Bloomington (main)]↗
DOE OSTI2024

Null Raychaudhuri: canonical structure and the dressing time

Abstract We initiate a study of gravity focusing on generic null hypersurfaces, non-perturbatively in the Newton coupling. We present an off-shell account of the extended phase space of the theory, which includes the expected spin-2 data as well as spin-0, spin-1 and arbitrary matter degrees of freedom. We construct the charges and the corresponding kinematic Poisson brackets, employing a Beltrami parameterization of the spin-2 modes. We explicitly show that the constraint algebra closes, the details of which depend on the non-perturbative mixing between spin-0 and spin-2 modes. Finally we show that the spin zero sector encodes a notion of a clock, called dressing time, which is dynamical and conjugate to the constraint. It is well-known that the null Raychaudhuri equation describes how the geometric data of a null hypersurface evolve in null time in response to gravitational radiation and external matter. Our analysis leads to three complementary viewpoints on this equation. First, it can be understood as a Carrollian stress tensor conservation equation. Second, we construct spin-0, spin-2 and matter stress tensors that act as generators of null time reparametrizations for each sector. This leads to the perspective that the null Raychaudhuri equation can be understood as imposing that the sum of CFT-like stress tensors vanishes. Third, we solve the Raychaudhuri constraint non-perturbatively. The solution relates the dressing time to the spin-2 and matter boost charge operators. Finally we establish that the corner charge corresponding to the boost operator in the dressing time frame is monotonic. These results show that the notion of an observer can be thought of as emerging from the gravitational degrees of freedom themselves. We briefly mention that the construction offers new insights into focusing conjectures.

Physics↗
DOE OSTI2024

DNA-mediated assembly of Au bipyramids into anisotropic light emitting kagome superlattices

Colloidal crystal engineering with DNA allows one to design diverse superlattices with tunable lattice symmetry, composition, and spacing. Most of these structures follow the complementary contact model, maximizing DNA hybridization on building blocks and producing relatively close-packed lattices. Here, low-symmetry kagome superlattices are assembled from DNA-modified gold bipyramids that can engage only in partial DNA surface matching. The bipyramid dimensions and DNA length can be engineered for two different superlattices with rhombohedral unit cells, including one composed of a periodic stacking of kagome lattices. Enabled by the partial facet alignment, the kagome lattices exhibit lattice distortion, bipyramid twisting, and planar chirality. When conjugated with Cy-5 dyes, the kagome lattices serve as cavities with high-density optical states and large Purcell factors along lateral directions, leading to strong dipole radiation along thezaxis and facet-dependent light emission. Such complex optical properties make these materials attractive for lasers, displays, and quantum sensing constructs.

Science & Technology - Other Topics↗
DOE OSTI2025

Orbital Engineering Band Degeneracy in a Dual-Square Carbon-Oxide Framework

Electron band degeneracies in momentum space give rise to exotic quantum phenomena that have sparked intense interest in condensed matter physics and materials science. Nodal-lines─isolines in k-space formed by the incidental touching of two bands that share the same energy but belong to discrete eigenstates─arise in the presence of symmetries that preclude effective hybridization. Despite recent advances in the design, bottom-up assembly, and engineering of exotic electronic states in graphene nanomaterials, the extension of this approach to access synthetic two-dimensional (2D) quantum materials derived from metal- or covalent-organic frameworks (COFs) has lagged behind. Here we present a molecular orbital engineering approach for designing and fabricating an edge-centered dual square lattice within a π-conjugated 2D-tetraoxa[8]circulene (2D-TOC) COF. First-principles calculations and scanning tunnelling spectroscopy reveal the emergence of Frontier states at the center of a 3 × 3 lattice that give rise to Dirac nodal-lines in 2D-TOC. Our findings not only provide a general guide for the design of conjugated COFs with custom tailored electronic properties from molecular fragments but enable the exploration of emergent topological phenomena in synthetic 2D materials with potential application for high-speed, low-power data processing, transmission, and storage.

Dirac nodal line semimetal↗
DOE OSTI2024

Using Molecular Structure to Tune Intrachain and Interchain Charge Transport in Indacenodithiophene-Based Copolymers

Here, in this work, we compare two structurally near-amorphous rigid-rod polymers-poly(indacenodithiophene-co-benzothiadiazole), p(IDT-BT), and poly(indacenodithiophene-co-benzopyrollodione), p(IDT-BPD)-with orders of magnitude different mobilities to understand the effect charge carrier intrachain delocalization has on electronic transport. Quantum chemical calculations show that p(IDT-BPD) has a barrier to torsion that is significantly lower than that of p(IDT-BT) and is thus more likely to have reduced conjugation lengths. We utilize absorption and photoluminescence spectroscopy to characterize energetic disorder and show that p(IDT-BPD) has higher energetic disorder. Charge modulation spectroscopy (CMS) and model calculations are used to show that charge carriers are substantially delocalized in p(IDT-BT) and occupy near-uniform energetic environments. We find that mobility activated hopping barriers are similar in these two materials. Electronic structure calculations show that both intrachain and interchain couplings of monomer units are poor enough in p(IDT-BPD) that charge carriers collapse to single IDT units and transport via a through-space tunneling mechanism. This work highlights the remarkable charge transport properties of p(IDT-BT) by showing that high mobilities are achievable on device-relevant length scales with only 1D carrier delocalization.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantum simulation of QED in Coulomb gauge

A recent work considered quantum simulation of Quantum Electrodynamics on a lattice in the Coulomb gauge with gauge degrees of freedom represented in the occupation basis in momentum space. In this work, we consider the more efficient representation of the gauge degrees of freedom in field basis in position space and develop a quantum algorithm for real-time simulation. We show that the continuum Coulomb gauge Hamiltonian is equivalent to the temporal gauge Hamiltonian when acting on physical states consisting of fermion and transverse gauge fields. The Coulomb gauge Hamiltonian is discretized by using the Green’s function of the discrete Laplacian operator under the Dirichlet boundary conditions. Both the continuum Coulomb gauge Hamiltonian and the discretized one proposed here guarantee that the unphysical longitudinal gauge fields are decoupled and commute with the corresponding Hamiltonian. Thus there is no need to impose any constraint. The local gauge field basis and the canonically conjugate variable basis are swapped efficiently using the quantum Fourier transform. We prove that the qubit cost to represent physical states and the gate count for real-time simulation scale polynomially with the lattice size, energy, time, accuracy, and Hamiltonian parameters in lattice units. The gate cost here for implementing the time evolution of the gauge field is reduced at least by a factor on the order of 10 8 for modest lattice size and accuracy level compared with the previous work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Comparing the Effects of Side Chain Dipole–Dipole Interactions and Hydrogen Bonding on the Mechanical and Electrical Properties of Poly(3-hexylthiophene)

Poly(3-alkylthiophenes) are simple conjugated polymers with good electrical properties, but tend to be brittle, limiting their application. Here, we investigate how side chain modifications can improve their toughness without compromising their charge transport. Specifically, we compare a weak dipole–dipole interaction of ester groups with a stronger hydrogen-bonding interaction of the hydroxyl groups. Two copolymers with targeted 5 and 10 mol % ester-functionalized side chains were synthesized and hydrolyzed in the solid state to introduce hydroxyl groups. Both hydrolyzed polymers became insoluble in most organic solvents, providing a path to layer-by-layer solution processing. All polymers’ optical, thermal, structural, electrical, and mechanical properties were investigated. Optical spectroscopy revealed that functionalization at these low levels did not significantly affect aggregation or π–π stacking. Differential scanning calorimetry indicated reduced crystallinity in the functionalized copolymers compared to that of P3HT. Grazing-incidence wide-angle X-ray diffraction data showed that functionalization led to an increase in lamellar spacing without altering molecular orientation or π–π stacking. Mechanical testing highlighted that at the lower functionalization level, both ester and hydroxyl groups similarly enhanced toughness, suggesting that disorder, rather than the specific nature of the functional group, was responsible for these improvements. At the higher functionalization level, hydroxyl groups significantly increased the elastic modulus while also negatively impacting charge carrier mobility. This study suggests that with stronger hydrogen-bonding groups, only very small amounts can improve P3HT’s toughness. For a wider tunable range, it may be a better strategy to use weak interactions such as dipole–dipole interactions of ester groups.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Parametric Study of Panel PCM–Air Heat Exchanger Designs

Heat exchangers, devices for the transfer of heat between two or more working fluids, are extensively used in cooling applications and heating applications. Heat exchangers in buildings are typically components of space-conditioning systems, as well as of water-heating applications. Heat exchangers are also sometimes used in applications that require storage and release of energy at specific times. Phase change materials (PCMs) enhance these heat-exchange processes, given their ability to melt and solidify at a fixed range of temperatures, absorbing or releasing significant amounts of latent heat. Five different configurations of PCM–air heat exchangers for thermal control in buildings are analyzed in this work. The heat exchangers were fitted with PCM encapsulated in plastic and composite pouches of various shapes, and packaged in stackable panel layers. Three-dimensional computational fluid dynamics (CFD) modeling of coupled incompressible fluid and conjugate heat transfer were performed on the designs. The phase change process was numerically modelled using the apparent heat capacity method. Steady-state CFD simulations provided quantification of pressure drop as a function of air flow velocity. Transient simulation results describe the thermal evolution of PCM in the pouches, helping to determine the best performing configuration with respect to total thermal charging time.

10 SYNTHETIC FUELS↗
DOE OSTI2025

Taming the virtual space for incremental full configuration interaction

Incremental full configuration interaction (iFCI) closely approximates the FCI limit with polynomial cost through a many-body expansion of the correlation energy, providing highly accurate total energies within a given basis set. To extend iFCI beyond previous basis set limitations, this work introduces a novel natural orbital (NO) screening approach, incremental NO full configuration interaction (iNO-FCI). By consideration of the importance of virtual orbital selection in the convergence of iFCI, iNO-FCI maximizes the consistency between orbitals selected for each correlated body. iNO-FCI employs a principle of cancellation of errors and ensures that the same set of virtual NOs is used for interdependent terms. Here, this strategy significantly reduces computational cost without compromising precision. Computational savings of up to 95% are demonstrated, allowing access to larger basis sets that were previously computationally prohibitive. iNO-FCI is herein introduced and benchmarked for several difficult test cases involving double-bond dissociation, biradical systems, conjugated π systems, and the spin gap of a Cu-based transition metal complex.

Correlation energy↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

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

Higher-rank tensor field theory of non-abelian fracton and embeddon

We formulate a new class of tensor gauge field theories in any dimension that is a hybrid class between symmetric higher-rank tensor gauge theory (i.e., higher-spin gauge theory) and anti-symmetric tensor topological field theory. Our theory describes a mixed unitary phase interplaying between gapless and gapped topological order phases (which can live with or without Euclidean, Poincaré, or anisotropic symmetry, at least in ultraviolet high or intermediate energy field theory, but not yet to a lattice cutoff scale). The “gauge structure” can be compact, continuous, abelian or non-abelian. Our theory sits outside the paradigm of Maxwell electromagnetic theory in 1865 and Yang–Mills isospin/color theory in 1954. We discuss its local gauge transformation in terms of the ungauged vector-like or tensor-like higher-moment global symmetry. The non-abelian gauge structure is caused by gauging the non-commutative symmetries: a higher-moment symmetry and a charge conjugation (particle–hole) symmetry. Vector global symmetries along time direction may exhibit time crystals. We explore the relation of these long-range entangled matters to a non-abelian generalization of Fracton order in condensed matter, a field theory formulation of foliation, the spacetime embedding and Embeddon that we newly introduce, and possible fundamental physics applications to dark matter or dark energy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Photo-Induced Charge Transfer of Fullerene and Non-Fullerene Conjugated Polymer Blends via Ab Initio Excited-State Dynamics

Organic conjugated polymers (CPs) are promising candidates for organic photovoltaic (OPV) devices due to their unique tunable mechanical and optoelectronic performance. Over the last decade, optoelectronic properties of narrow band gap CPs as a blend with acceptor units are largely optimized, which leads to noticeable progress in OPV technology. However, their power conversion efficiency is still lower than their organic counterparts (i.e., silicon), limiting their practical usage. In this study, we employ ab initio molecular dynamics to explore photo-induced charge transfer (CT) of the diketopyrrolopyrrole-based polymer as a blend with non-fullerene (i.e., ITIC) and fullerene (i.e., PCBM) acceptor units. The results of charge carrier dynamics induced by selected photoexcitation show that hole density redistribution in space is much faster than electron relaxation. We track the relaxation rates of charge carriers over time, where the derivative of the difference between the rate of electron and hole implies the current density at zero voltage. Furthermore, this can be utilized to characterize the CT performance of CPs blended with different acceptor units. Relaxation rate results indicate that CP blend with ITIC promises a better PV performance, illustrating that the current computational approach opens the door to determine bulk heterojunctions’ electronic performance for OPV devices and narrowing down the list of potential donor–acceptor candidates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Potent Inhibition of E. coli DXP Synthase by a gem -Diaryl Bisubstrate Analog

New antimicrobial strategies are needed to address pathogen resistance to currently used antibiotics. Bacterial central metabolism is a promising target space for the development of agents that selectively target bacterial pathogens. 1-Deoxy- D -xylulose 5-phosphate synthase (DXPS) converts pyruvate and d-glyceraldehyde 3-phosphate ( D -GAP) to DXP, which is required for synthesis of essential vitamins and isoprenoids in bacterial pathogens. Thus, DXPS is a promising antimicrobial target. Toward this goal, our lab has demonstrated selective inhibition of Escherichia coli DXPS by alkyl acetylphosphonate (alkylAP)-based bisubstrate analogs that exploit the requirement for ternary complex formation in the DXPS mechanism. Here, we present the first DXPS structure with a bisubstrate analog bound in the active site. Insights gained from this cocrystal structure guided structure–activity relationship studies of the bisubstrate scaffold. A low nanomolar inhibitor (compound 8) bearing a gem-dibenzyl glycine moiety conjugated to the acetylphosphonate pyruvate mimic via a triazole-based linker emerged from this study. Compound 8 was found to exhibit slow, tight-binding inhibition, with contacts to E. coli DXPS residues R99 and R478 demonstrated to be important for this behavior. This work has discovered the most potent DXPS inhibitor to date and highlights a new role of R99 that can be exploited in future inhibitor designs toward the development of a novel class of antimicrobial agents.

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

Discovery of Self-Assembling π-Conjugated Peptides by Active Learning-Directed Coarse-Grained Molecular Simulation

Electronically active organic molecules have demonstrated great promise as novel soft materials for energy harvesting and transport. Self-assembled nanoaggregates formed from π-conjugated oligopeptides composed of an aromatic core flanked by oligopeptide wings offer emergent optoelectronic properties within a water-soluble and biocompatible substrate. Nanoaggregate properties can be controlled by tuning core chemistry and peptide composition, but the sequence–structure–function relations remain poorly characterized. Here, we employ coarse-grained molecular dynamics simulations within an active learning protocol employing deep representational learning and Bayesian optimization to efficiently identify molecules capable of assembling pseudo-1D nanoaggregates with good stacking of the electronically active π-cores. We consider the DXXX-OPV3-XXXD oligopeptide family, where D is an Asp residue and OPV3 is an oligophenylenevinylene oligomer (1,4-distyrylbenzene), to identify the top performing XXX tripeptides within all 20 3 = 8000 possible sequences. By direct simulation of only 2.3% of this space, we identify molecules predicted to exhibit superior assembly relative to those reported in prior work. Spectral clustering of the top candidates reveals new design rules governing assembly. This work establishes new understanding of DXXX-OPV3-XXXD assembly, identifies promising new candidates for experimental testing, and presents a computational design platform that can be generically extended to other peptide-based and peptide-like systems.

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

A Scalable Interior‐Point Gauss–Newton Method for PDE‐Constrained Optimization With Bound Constraints

Here, we present a scalable approach to solve a class of partial differential equation (PDE)‐constrained optimization problems with bound constraints. This approach utilizes a robust full‐space interior‐point (IP)‐Gauss–Newton optimization method. To cope with the poorly‐conditioned IP‐Gauss–Newton saddle‐point linear systems that need to be solved approximately, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE‐constrained optimization problems. A block Gauss–Seidel preconditioner is proposed for the GMRES‐based solution of the IP‐Gauss–Newton linear systems. It is shown, for a large‐class of PDE‐ and bound‐constrained optimization problems, that the spectrum of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix is asymptotically independent of discretization and is not impacted by the ill‐conditioning that notoriously plagues interior‐point methods. We exploit symmetry of the IP‐Gauss–Newton linear systems and propose a regularization and log‐barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)‐based solution of the equivalent IP‐Gauss–Newton–Schur complement linear systems. The eigenvalues of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix, that are not equal to one, are identical to the eigenvalues of the regularization and log‐barrier Hessian preconditioned Schur complement matrix. The scalability of the approach is demonstrated on two example problems. The numerical solution of these optimization problems is shown to require a discretization independent number of IP‐Gauss–Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill‐conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of the preconditioner, achieved via algebraic multigrid component solvers when applicable, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE‐ and bound‐constrained problems.

PDE-constrained optimization↗