Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

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 1,459 records · Page 81

Nanoparticle self-assemblies with modern complexity

Thanks to decades of tireless efforts, nanoparticle assemblies have reached at an extremely high level of controllability, sophistication, and complexity, with new insights provided by integration with graph theory, cutting-edge characterization, and machine learning (ML)-based computation and modeling, as well as with ever-diversifying applications in energy, catalysis, biomedicine, optics, electronics, magnetics, organic biosynthesis, and quantum technology. Nanoparticle assemblies can be crystalline, known as superlattices or supracrystals. Their assembly entails a transition from disorder—dispersed nanoparticles—to order, which can be achieved through classical nucleation pathways or nonclassical pathways via prenucleation precursors or particle aggregation. Further, the periodic lattices allow facile manipulations of electrons, phonons, photons, and even spins, leading to advanced device components and metamaterials. Meanwhile, aperiodic assemblies out of nanoparticles, such as gels, networks, and amorphous solids, also start to attract attentions. Despite the loss of periodicity, symmetry-lowering or symmetry-breaking three-dimensional (3D) structures emerge with unique properties, such as chiroptical activity, topological mechanical strength, and quantum entanglement. Real-space imaging such as electron microscopy and X-ray based tomography methods are utilized to characterize these complex structures, while mathematical tools such as graph theories are in need to describe such complex structures. This issue aims to provide a timely review of the efforts in this greatly broadened materials design space including experiment, simulation, theory, and applications. Nine top experts (and their teams) from four countries deliver six review papers, summarizing fundamental mechanistic understandings of nanoparticle assemblies, highlighted with the developments of state-of-the-art in situ characterization tools and ML-assisted reverse engineering, and newly emergent applications of nanoarchitectures.

36 MATERIALS SCIENCE↗

FINAL REPORT

Hydrogen storage systems based on the readily reversible adsorption of H{sub 2} in porous materials have a number of very attractive properties with the potential to provide superior performance among candidate materials currently being investigated were it not for the fact that the interaction of H{sub 2} with the host material is too weak to permit viable operation at room temperature. Our study has delineated in quantitative detail the structural elements which we believe to be the essential ingredients for the future synthesis of porous materials, where guest-host interactions are intermediate between those found in the carbons and the metal hydrides, i.e. between physisorption and chemisorption, which will result in H{sub 2} binding energies required for room temperature operation. The ability to produce porous materials with much improved hydrogen binding energies depends critically on detailed molecular level analysis of hydrogen binding in such materials. However, characterization of H{sub 2} sorption is almost exclusively carried by thermodynamic measurements, which give average properties for all the sites occupied by H{sub 2} molecules at a particular loading. We have therefore extensively utilized the most powerful of the few molecular level experimental probes available to probe the interactions of hydrogen with porous materials, namely inelastic neutron scattering (INS) spectroscopy of the hindered rotations of the hydrogen molecules adsorbed at various sites, which in turn can be interpreted in a very direct way in by computational studies. This technique can relate spectral signatures of various H{sub 2} molecules adsorbed at binding sites with different degrees of interaction. In the course of this project we have synthesized a rather large number of entirely new hybrid materials, which include structural modifications for improved interactions with adsorbed hydrogen. The results of our systematic studies on many porous materials provide detailed information on the effects on hydrogen binding from framework modifications, including charged frameworks and extraframework cations, from reduction in pore sizes, functionalization of the organic linking group, and most importantly, that of the various types of metal sites. We provided a clear demonstration that metal sites are most effective if the metal is highly undercoordinated, open and completely accessible to the H{sub 2} molecule, a condition which is not currently met in MOFs with intra-framework metals. The results obtained from this project therefore will give detailed direction to efforts in the synthesis of new materials that can reach the goal of a practical sorption based hydrogen storage material.

08 HYDROGEN↗

Reactivation of Kamb Ice Stream tributaries triggers century-scale reorganization of Siple Coast ice flow in West Antarctica

Ongoing, centennial-scale flow variability within the Ross ice streams of West Antarctica suggests that the present-day positive mass balance in this region may reverse in the future. Here we use a three-dimensional ice sheet model to simulate ice flow in this region over 250 years. The flow responds to changing basal properties, as a subglacial till layer interacts with water transported in an active subglacial hydrological system. We show that a persistent weak bed beneath the tributaries of the dormant Kamb Ice Stream is a source of internal ice flow instability, which reorganizes all ice streams in this region, leading to a reduced (positive) mass balance within decades and a net loss of ice within two centuries. This hitherto unaccounted for flow variability could raise sea level by 5 mm this century. Furthermore, better constraints on future sea level change from this region will require improved estimates of geothermal heat flux and subglacial water transport.

58 GEOSCIENCES↗

Origin and evolution of HIV-1 subtype A6

Background: HIV outbreaks in the Former Soviet Union (FSU) countries were characterized by repeated transmission of the HIV variant AFSU, which is now classified as a distinct subtype A sub-subtype called A6. The current study used phylogenetic/phylodynamic and signature mutation analyses to determine likely evolutionary relationship between subtype A6 and other subtype A sub-subtypes. Methods: For this study, an initial Maximum Likelihood phylogenetic analysis was performed using a total of 553 full-length, publicly available, reverse transcriptase sequences, from A1, A2, A3, A4, A5, and A6 sub-subtypes of subtype A. For phylogenetic clustering and signature mutation analysis, a total of 5961 and 3959 pol and env sequences, respectively, were used. Results: Phylogenetic and signature mutation analysis showed that HIV-1 sub-subtype A6 likely originated from sub-subtype A1 of African origin. A6 and A1 pol and env genes shared several signature mutations that indicate genetic similarity between the two subtypes. For A6, tMRCA dated to 1975, 15 years later than that of A1. Conclusion: The current study provides insights into the evolution and diversification of A6 in the backdrop of FSU countries and indicates that A6 in FSU countries evolved from A1 of African origin and is getting bridged outside the FSU region.

60 APPLIED LIFE SCIENCES↗

Understanding Proton Movement in [Fe-Fe] Hydrogenases

Nature uses specialized metalloenzymes to carry out small molecule activation reactions, including CO 2 fixation, O 2 activation, and proton reduction, with unparalleled efficiency, rates, and selectivity. The latter reactions are performed by hydrogenases, protein metallo-complexes that interconverts H 2 to protons and electrons (H 2 oxidation) and the reverse reaction (H 2 production) with incredibly low energy input and amazingly fast kinetics. Reproducing both the activity and efficiency of metalloenzymes in sustainable anthropogenic systems remains one of the “holy grails” of inorganic chemistry. However, identifying the precise molecular components responsible for these desirable properties has been challenging in the natural metalloenzymes, hindering efforts to develop analogous processes in synthetic compounds. Considering the inherent complexity of a metalloenzyme and the many interactions, both strong and weak, that contribute to the function of an enzyme, we have elected to model natural metalloenzymes on a biochemical platform. Towards this end, we have a developed structural, functional, and mechanistic mimic of the [Ni-Fe] hydrogenases within a robust protein scaffold, rubredoxin, to understand the influence of the secondary coordination environment on the metal center. This involved preparing a series of three rubredoxin constructs containing a single point mutation at Val positions and making the NMR chemical shift assignments for the paramagnetic (nickel-substituted) and non-paramagnetic (zinc-substituted) form. These physical studies were complemented with in silico molecular dynamic studies on proteins with [Fe-S] clusters to engineer new proton channels to test in vitro . To assist our search for new natural metalloprotein scaffolds within the vast number of sequenced genomes, we created a neural network-based program to identify proteins with specific metal-binding sites. These computational and physical studies with metalloenzymes provide direct insight into the fundamental chemical principles driving the natural systems and offer design principles for developing catalysts that utilize analogous principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deformation of Lie–Poisson algebras and chirality

Linearization of a Hamiltonian system around an equilibrium point yields a set of Hamiltonian symmetric spectra: If λ is an eigenvalue of the linearized generator, -λ and $\barλ$ (hence, -$\barλ$) are also eigenvalues—the former implies a time-reversal symmetry, while the latter guarantees the reality of the solution. However, linearization around a singular equilibrium point (which commonly exists in noncanonical Hamiltonian systems) works out differently, resulting in breaking of the Hamiltonian symmetry of spectra; time-reversal asymmetry causes chirality. This interesting phenomenon was first found in analyzing the chiral motion of the rattleback, a boat-shaped top having misaligned axes of inertia and geometry [Z. Yoshida et al., Phys. Lett. A 381, 2772–2777 (2017)]. To elucidate how chiral spectra are generated, we study the three-dimensional Lie–Poisson systems and classify the prototypes of singularities that cause symmetry breaking. The central idea is the deformation of the underlying Lie algebra; invoking Bianchi’s list of all three-dimensional Lie algebras, we show that the so-called class-B algebras, which are produced by asymmetric deformations of the simple algebra s o ( 3 ) , yield chiral spectra when linearized around their singularities. The theory of deformation is generalized to higher dimensions, including the infinite-dimensional Poisson manifolds relevant to fluid mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian Inference with Latent Hamiltonian Neural Networks (L-HNNs)

When sampling for Bayesian inference, one popular approach is to use Hamiltonian Monte Carlo (HMC) and the No-U-Turn Sampler (NUTS). However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice. We propose Hamiltonian neural networks (HNNs) with HMC and NUTS for solving Bayesian inference problems [1, 2]. Once trained, HNNs do not require gradients of the target density while sampling. Moreover, they satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well suited for use within HMC and NUTS because stationarity can be shown. We also propose an HNN extension called latent HNNs (L-HNNs), which predict latent variable outputs. Compared to HNNs, L-HNNs offer improved expressivity and a reduction in integration errors. Finally, we propose employing L-HNNs in NUTS with an online error monitoring scheme to prevent degeneracy of the sampling in regions of low probability density. We demonstrate L-HNNs in NUTS with online error monitoring by using several example cases involving complex, heavy-tailed, and high local curvature probability densities. Overall, L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. Compared to traditional NUTS, L-HNNs in NUTS with online error monitoring improved the effective sample size (ESS) per gradient by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗

Vortex Diode Analysis and Testing for Fluoride Salt-Cooled High-Temperature Reactors

Fluidic diodes are presently being considered for use in several fluoride salt-cooled high-temperature reactor designs. A fluidic diode is a passive device that acts as a leaky check valve. These devices are installed in emergency heat removal systems that are designed to passively remove reactor decay heat using natural circulation. The direct reactor auxiliary cooling system (DRACS) uses DRACS salt-to-salt heat exchangers (DHXs) that operate in a path parallel to the core flow. Because of this geometry, under normal operating conditions some flow bypasses the core and flows through the DHX. A flow diode, operating in reverse direction, is-used to minimize this flow when the primary coolant pumps are in operation, while allowing forward flow through the DHX under natural circulation conditions. The DRACSs reject the core decay heat to the environment under loss-of-flow accident conditions and as such are a reactor safety feature. Fluidic diodes have not previously been used in an operating reactor system, and therefore their characteristics must be quantified to ensure successful operation. This report parametrically examines multiple design parameters of a vortex-type fluidic diode to determine the size of diode needed to reject a particular amount of decay heat. Additional calculations were performed to size a scaled diode that could be tested in the Oak Ridge National Laboratory Liquid Salt Flow Loop. These parametric studies have shown that a 152.4 mm diode could be used as a test article in that facility. A design for this diode is developed, and changes to the loop that will be necessary to test the diode are discussed. Initial testing of a scaled flow diode has been carried out in a water loop. The 150 mm diode design discussed above was modified to improve performance, and the final design tested was a 171.45 mm diameter vortex diode. The results of this testing indicate that diodicities of about 20 can be obtained for diodes of this size. Experimental results show similar trends as the computational fluid dynamics (CFD) results presented in this report; however, some differences exist that will need to be assessed in future studies. The results of this testing will be used to improve the diode design to be tested in the liquid salt loop system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Fungal indole alkaloid biogenesis through evolution of a bifunctional reductase/Diels–Alderase

Prenylated indole alkaloids such as the calmodulin-inhibitory malbrancheamides and anthelmintic paraherquamides possess great structural diversity and pharmaceutical utility. Here, we report complete elucidation of the malbrancheamide biosynthetic pathway accomplished through complementary approaches. These include a biomimetic total synthesis to access the natural alkaloid and biosynthetic intermediates in racemic form and in vitro enzymatic reconstitution to provide access to the natural antipode (+)-malbrancheamide. Reductive cleavage of an l-Pro–l-Trp dipeptide from the MalG non-ribosomal peptide synthetase (NRPS) followed by reverse prenylation and a cascade of post-NRPS reactions culminates in an intramolecular [4+2] hetero-Diels–Alder (IMDA) cyclization to furnish the bicyclo[2.2.2]diazaoctane scaffold. Enzymatic assembly of optically pure (+)-premalbrancheamide involves an unexpected zwitterionic intermediate where MalC catalyses enantioselective cycloaddition as a bifunctional NADPH-dependent reductase/Diels–Alderase. Furthermore, the crystal structures of substrate and product complexes together with site-directed mutagenesis and molecular dynamics simulations demonstrate how MalC and PhqE (its homologue from the paraherquamide pathway) catalyse diastereo- and enantioselective cyclization in the construction of this important class of secondary metabolites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On-Site Pilot Study - Removal of Uranium, Radium-226 and Arsenic from Impacted Leachate by Reverse Osmosis - 13155

Conestoga-Rovers and Associates (CRA-LTD) performed an on-site pilot study at the Welcome Waste Management Facility in Port Hope, Ontario, Canada, to evaluate the effectiveness of a unique leachate treatment process for the removal of radioactive contaminants from leachate impacted by low-level radioactive waste. Results from the study also provided the parameters needed for the design of the CRA-LTD full scale leachate treatment process design. The final effluent water quality discharged from the process to meet the local surface water discharge criteria. A statistical software package was utilized to obtain the analysis of variance (ANOVA) for the results from design of experiment applied to determine the effect of the evaluated factors on the measured responses. The factors considered in the study were: percent of reverse osmosis permeate water recovery, influent coagulant dosage, and influent total dissolved solids (TDS) dosage. The measured responses evaluated were: operating time, average specific flux, and rejection of radioactive contaminants along with other elements. The ANOVA for the design of experiment results revealed that the operating time is affected by the percent water recovery to be achieved and the flocculant dosage over the range studied. The average specific flux and rejection for the radioactive contaminants were not affected by the factors evaluated over the range studied. The 3 month long on-site pilot testing on the impacted leachate revealed that the CRA-LTD leachate treatment process was robust and produced an effluent water quality that met the surface water discharge criteria mandated by the Canadian Nuclear Safety Commission and the local municipality. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Outer shell symmetry for double shell capsules with aluminum ablators

Double shell targets are a promising potential avenue to obtain robust neutron yield at current laser facilities. Similar to single shell designs, double shells require the symmetric implosion of an ablator in order to uniformly compress and heat a fuel volume, with the goal of achieving thermonuclear burn. Significant differences between double and single shells include the usage of an aluminum ablator as well as a reverse ramp laser pulse. In addition, double shells require a different convergence than single shells for fuel ignition. Numerical implosion studies at various energies with comparisons to experimental outcomes are required to gain confidence that simulations can capture the ablator shape from subscale to full scale. The current work builds on previous implosion experiments conducted at 1-MJ laser energy to confirm achieved ablator symmetry at 1.25 and 1.5 MJ. Average ablator P2 and P4 shapes measured in these experiments are within 5% of the simulated shape, which merits the platforms for further experimental studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Firmware Tampering Detection in Heavy-Duty Vehicles through J1939 CAN Analysis

Modern heavy-duty vehicles rely on complex networks of Electronic Control Units (ECUs) that communicate using the J1939 protocol. While this system makes it easier to update and configure vehicle components, it also opens the door to serious cybersecurity risks if not properly secured. This work investigates the potential for firmware tampering through the J1939 communication protocol, which enables ECU configuration and reprogramming over the Controller Area Network (CAN) bus. By monitoring CAN traffic during legitimate maintenance operations and reverse-engineering OEM diagnostic software, we identified common and proprietary J1939 message identifiers, authentication patterns, and vulnerabilities within Unified Diagnostic Services (UDS). These findings demonstrate that inadequate authentication mechanisms can allow malicious actors to alter ECU firmware or disable safety functions, posing severe operational and safety risks. Our analysis contributes to the development of vehicle intrusion detection systems capable of recognizing abnormal reprogramming activity and future firmware fingerprinting methods to verify software integrity across ECUs. This work highlights the importance of standardizing secure firmware authentication across manufacturers to strengthen cyber resilience in heavy-duty vehicle systems.

33 ADVANCED PROPULSION SYSTEMS↗

"Electrostructural Phase Changes" In Charged Particulate Clouds: Planetary and Astrophysical Implications

There is empirical evidence that freely-suspended triboelectrostatically charged particulate clouds of dielectric materials undergo rapid conversion from (nominally) monodispersed "aerosols" to a system of well-defined grain aggregates after grain motion or fluid turbulence ceases within the cloud. In United States Microgravity Laboratory Space Shuttle experiments USML-1 and USML-2, it was found that ballistically-energized grain dispersions would rapidly convert into populations of filamentary aggregates after natural fluid (air) damping of grain motion. Unless continuously disrupted mechanically, it was impossible to maintain a non-aggregated state for the grain clouds of sand-size materials. Similarly, ground- based experiments with very fine dust-size material produced the same results: rapid, impulsive "collapse" of the dispersed grains into well-defined filamentary structures. In both ground-based and microgravity experiments, the chains or filaments were created by long-range dipole electrostatic forces and dipole-induced dielectric interactions, not by monopole interactions. Maintenance of the structures was assisted by short-range static boundary adhesion forces and van der Waals interactions. When the aggregate containers in the USML experiments were disturbed after aggregate formation, the quiescently disposed filaments would rearrange themselves into fractal bundles and tighter clusters as a result of enforced encounters with one another. The long-range dipole interactions that bring the grains together into aggregates are a product of randomly-distributed monopole charges on the grain surfaces. In computer simulations, it has been shown that when the force vectors of all the random charges (of both sign) on a grain are resolved mathematically by assuming Coulombic interaction between them, the net result is a dipole moment on individual grains, even though the grains are electrically neutral insofar as there is no predominance, on their surface, of one charge sign over another. The random charges of both sign derive from natural grain-to-grain interactions that produce triboelectrification via charge exchange every time grain surfaces make contact with one another. The conversion from a random distribution of grains (upon which there are randomly distributed charges) into an organization of electrostatically-ordered aggregates, can be regarded (within the framework of granular-material science) as an "electrical or Coulombic phase change" of the particulate cloud. It is not totally dissimilar from the more normal phase-change concept in which, for example, a gas with long free-path-molecules suddenly becomes a solid as a result of structural ordering of the molecules (notably, also the result of electronic forces, albeit at a different scale). In both the gas-to-solid case, and the aerosol-to-aggregate case, the same materials and charges are present before and after the phase change, but their arrangement now has a higher degree of order and a lower-energy configuration. An input of energy into the system is required to reverse the situation. The aggregates in the USML experiments were observed to undergo at least two phase changes as noted above. The point about phase changes, and by implication, the "electrostructural" reorganizations in particulate clouds, is the following: (a) they can occur very rapidly, almost spontaneously, above a critical cloud density, (b) in going from a higher energy state to a lower energy state, they convert to a denser system, (c) energy must be required to reverse the situation, implying that energy is released during the high-to-low energy phase change. In applying this information to natural particulate clouds, some inferences can be made (it is stressed that reference is still to dielectric materials attracted by dipole forces). There are several natural settings to which the USML observations apply, and to which the phase-change implications likewise apply. Dense clouds of triboelectrically-charged, kinetically-energized grains are to be found in volcanic eruptions (particularly on earth), aeolian dust storms (particularly on Mars), meteorite impact ejecta curtains (on all planets), in "immature" debris rings around planets (e.g., that from which our own Moon may have condensed), and in gravitationally collapsing protoplanetary dust/planetesimal debris disks where dielectric granules are being increasingly brought into collisional relationships with one another (increasing both charge exchange and physical proximity). It is noted that in many of these cases, the degree of electrical charging on the grains is likely to be much higher than that in the USML experiments where charging was not enhanced above the "normal", naturally encountered level for the particular materials (quartz, glass, and various silicate minerals). Application of the phase-change concept suggests that volcanic, aeolian, and impact debris clouds may, under certain circumstances, undergo rapid, impulsive, or even catastrophic collapse into a denser state that will lead to rapid precipitation or fall-out of suspended particulates. Although this idea has been suggested previously by the author , the phase change concept possibly permits some new insights into cloud-system behavior. For example, in a protoplanetary debris disk, the work of gravity may suddenly be enhanced by electrostatically driven collapse of the system when materials reach a critical intergranular spacing or grain number density. This might reduce the rate of planet formation by orders of magnitude (considering the ratio of g-forces to electrostatic ones for very small grains in close proximity), and indeed, it might drive the collapse into a system configuration that would perhaps not be created by gravity alone. Additional information is contained in the original.

Marshall, J. R.↗

Numerical Investigation of Butterfly Valve Performance in Variable Valve Sizes, Positions and Flow Regimes

Reliability and efficiency of valves are necessary for precise control and sufficient heat-flow to heat application plants for the integrated energy systems of nuclear power plants (NPPs). Strategic Management Analysis Requirement and Technology (SMART) valves’ ability to control flow and assess environmental parameters stands out for these requirements. Their ability to sustain the downstream flow rate, prevent reverse flow, and maintain pressure in the heat transport loop is much more efficient with the integration of sensors and intelligent algorithms. For assessing valve performance and monitoring, mechanical design and operating conditions are two important parameters. In this study, the butterfly valves of three different sizes are simulated with water and steam using STAR-CCM+ in various flow regimes and positions to analyze performance parameters to strategize an automated control system for efficiently balancing the heat–transport network. Also, flow behavior is studied using velocity and pressure fields for valve–body geometry optimization. It can be observed, through performance parameters, that the valves are suitable for operation between 30° and 90° positions with significantly low loss coefficients and high flow coefficients, and the performance parameters follow a certain pattern in both water and steam flow in each scenario.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sequential Kalman tuning of the t -preconditioned Crank-Nicolson algorithm: efficient, adaptive and gradient-free inference for Bayesian inverse problems

Ensemble Kalman Inversion (EKI) has been proposed as an efficient method for the approximate solution of Bayesian inverse problems with expensive forward models. However, when applied to the Bayesian inverse problem EKI is only exact in the regime of Gaussian target measures and linear forward models. Here, in this work we propose embedding EKI and Flow Annealed Kalman Inversion, its normalizing flow (NF) preconditioned variant, within a Bayesian annealing scheme as part of an adaptive implementation of the t-preconditioned Crank-Nicolson (tpCN) sampler. The tpCN sampler differs from standard pCN in that its proposal is reversible with respect to the multivariate t-distribution. The more flexible tail behaviour allows for better adaptation to sampling from non-Gaussian targets. Within our Sequential Kalman Tuning (SKT) adaptation scheme, EKI is used to initialize and precondition the tpCN sampler for each annealed target. The subsequent tpCN iterations ensure particles are correctly distributed according to each annealed target, avoiding the accumulation of errors that would otherwise impact EKI. We demonstrate the performance of SKT for tpCN on three challenging numerical benchmarks, showing significant improvements in the rate of convergence compared to adaptation within standard SMC with importance weighted resampling at each temperature level, and compared to similar adaptive implementations of standard pCN. The SKT scheme applied to tpCN offers an efficient, practical solution for solving the Bayesian inverse problem when gradients of the forward model are not available. Code implementing the SKT schemes for tpCN is available at https://github.com/RichardGrumitt/KalmanMC.

97 MATHEMATICS AND COMPUTING↗

Approximate singly excited states from a two-component Hartree-Fock reference

For many molecules, relaxing the spin symmetry constraint on the wave function results in the lowest energy mean-field solution. The two-component Hartree-Fock (2cHF) method relaxes all spin symmetry constraints, and the wave function is no longer an eigenfunction of the total spin, spin projection, or time-reversal symmetry operators. For ground state energies, 2cHF is a superior mean-field method for describing spin-frustrated molecules. For excited states, the utility of 2cHF is uncertain. Here, we implement the 2cHF extensions of two single-reference excited state methods, the two-component configuration interaction singles and time-dependent Hartree-Fock. We compare the results to the analogous methods based off of the unrestricted Hartree-Fock approximation, as well as the full configuration interaction for three small molecules with distinct 2cHF solutions, and discuss the nature of the 2cHF excited state solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic Capacity of Quantum Processes

Thermodynamics imposes restrictions on what state transformations are possible. In the macroscopic limit of asymptotically many independent copies of a state—as for instance in the case of an ideal gas—the possible transformations become reversible and are fully characterized by the free energy. In this Letter, we present a thermodynamic resource theory for quantum processes that also becomes reversible in the macroscopic limit, a property that is especially rare for a resource theory of quantum channels. We identify a unique single-letter and additive quantity, the thermodynamic capacity, that characterizes the “thermodynamic value” of a quantum channel, in the sense that the work required to simulate many repetitions of a quantum process employing many repetitions of another quantum process becomes equal to the difference of the respective thermodynamic capacities. On a technical level, we provide asymptotically optimal constructions of universal implementations of quantum processes. A challenging aspect of this construction is the apparent necessity to coherently combine thermal engines that would run in different thermodynamic regimes depending on the input state. Furthermore, our results have applications in quantum Shannon theory by providing a generalized notion of quantum typical subspaces and by giving an operational interpretation to the entropy difference of a channel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗