SEARCH · Search NASA

Results for “weak solutions”

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 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2022

VWC-BERT: Scaling Vulnerability–Weakness–Exploit Mapping on Modern AI Accelerators

Defending cybersystems needs accurate mapping of software and hardware vulnerabilities to generalized descriptions of weaknesses, and weaknesses to exploits. These mappings enable cyber defenders to build plans for effective defense and assessment of potential risks to a cybersystem. With close to 170k vulnerabilities, manual mapping is not a feasible option. However, automated mapping is challenging due to limited training data, computational intractability, and limitations in computational natural language processing. Tools based on breakthroughs in Transformer-based language models have been demonstrated to classify vulnerabilities with high accuracy. We make three key contributions in this paper: (1) We present a new framework, \VWCBERT, that augments the Transformer-based hierarchical multi-class classification framework of Das et al. (\textsc{V2W-BERT}) with the ability to map weaknesses to exploits. (2) We implement \VWCBERT~ on modern AI accelerator platforms using two data parallel techniques for the pre-training phase and demonstrate nearly linear speedups across NVIDIA and Graphcore accelerator platforms. We observe nearly linear speedups for up to 16 V100 and 8 A100 GPUs, and about 3.4$\times$ speedup for A100 relative to V100 GPUs. We also observe excellent speedups on Graphcore, with $5.7\times$ speedup on 128 IPUs relative to 16 IPUs. Enabled by scaling, we also demonstrate higher accuracy using a larger language model, RoBERTa-Large. We show up to 87\% accuracy for strict and up to 98\% accuracy for relaxed classification. (3) We develop a novel parallel link manager for the link prediction phase and demonstrate up to 21$\times$ speedup with 16 V100 GPUs relative to one V100 GPU, and thus reducing the runtime from 2.5 hours to 10 minutes. We believe that generalizability and scalability of \VWCBERT~ will benefit both the theoretical development and practical deployment of novel cyberdefense solutions and vulnerability classification.

Das, Siddhartha Shankar↗

DOE OSTI2020

A more natural composite Higgs model

Composite Higgs models provide an attractive solution to the hierarchy problem. However, many realistic models suffer from tuning problems in the Higgs potential. There are often large contributions from the UV dynamics of the composite resonances to the Higgs potential, and tuning between the quadratic term and the quartic term is required to separate the electroweak breaking scale and the compositeness scale. We consider a composite Higgs model based on the SU(6)/Sp(6) coset, where an enhanced symmetry on the fermion resonances can minimize the Higgs quadratic term. Moreover, a Higgs quartic term from the collective symmetry breaking of the little Higgs mechanism can be realized by the partial compositeness couplings between elementary Standard Model fermions and the composite operators, without introducing new elementary fields beyond the Standard Model and the composite sector. The model contains two Higgs doublets, as well as several additional pseudo-Nambu-Goldstone bosons. To avoid tuning, the extra Higgs bosons are expected to be relatively light and may be probed in the future LHC runs. The deviations of the Higgs couplings and the weak gauge boson couplings also provide important tests as they are expected to be close to the current limits in this model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

An OpenMP GPU-offload implementation of a non-equilibrium solidification cellular automata model for additive manufacturing

Here, in this paper, performance strategies on GPU-based HPC platforms of a cellular automata (CA) simulation code for non-equilibrium solidification, including nucleation, grain growth, solute partitioning and transport for the metal additive manufacturing (AM) process are investigated using OpenMP 4.5. To accurately report the speed-up for multicore CPUs and GPUs, a rigorous performance analysis employed optimizations appropriate for both CPU-only code (baseline) and GPU offload codes for an isothermal test problem. The performance results on Summit at the Oak Ridge Leadership Computing Facility indicate that using a precomputed list of interface cells significantly decreased the wall-clock time on GPUs. The speedup due to GPU acceleration was evaluated for a full Summit node and measured to be 1.8X when comparing a 6 MPI tasks run with 6 GPUs versus 36 MPI tasks on the CPU only. That speed-up was found to be 7.9X when comparing 6 MPI tasks with 6 GPUs versus the 6 MPI tasks running on the CPU only. Performance measurements showed that system total time is almost constant for runs with more than 96 MPI tasks (or GPUs), indicating that the GPU-accelerated code showed an excellent weak scaling performance. Finally, a rapid directional solidification problem was considered to demonstrate the CA code capability on Summit. It was found that a mesh size of at least 0.05 μm is recommended for the AM-like simulations in order to obtain accurate elongated grain microstructure and elongated subgrain features, which are in qualitative good agreement with experimental data. The results presented in this study indicate that the performance strategies on GPU-based HPC platforms for the CA code are appropriate for novel HPC exascale platforms.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae

The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than 50 years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi–Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Rare and Experimentally Challenging Supersymmetry Signatures

Supersymmetry has long played a central role in the search for physics beyond the Standard Model at colliders, providing a comprehensive and internally consistent framework for generating well-motivated experimental signatures. For more than 15 years of Large Hadron Collider (LHC) operation, the CMS and ATLAS Collaborations have achieved remarkable sensitivity to a wide range of supersymmetric signatures. Despite this unprecedented reach, no conclusive evidence for supersymmetry has emerged. If supersymmetry is nature's solution to outstanding questions in particle physics, it is necessarily challenging to find. In this article, we review supersymmetric signatures that are particularly rare or otherwise challenging, with a focus on searches at the LHC. We highlight experimental challenges related to detector constraints and analysis difficulties, in addition to model challenges in the interpretation and optimization of searches. Here, we also identify regions of signature space that remain comparatively unconstrained and therefore represent promising targets for future exploration.

Large Hadron Collider↗

DOE OSTI2020

Helicity at small x: oscillations generated by bringing back the quarks

We construct a numerical solution of the recently-derived large-N c &N f small-x helicity evolution equations with the aim to establish the small-x asymptotics of the quark helicity distribution beyond the large-N c limit explored previously in the same framework. (Here N c and N f are the numbers of quark colors and flavors.) While the large-N c helicity evolution involves gluons only, the large-N c &N f evolution includes contributions from quarks as well. We find that adding quarks to the evolution makes quark helicity distribution oscillate as a function of x. Our numerical results in the large-N c &N f limit lead to the x-dependence of the flavor-singlet quark helicity distribution which is well-approximated by $$ {\left.\Delta \Sigma \left(x,{Q}^2\right)\right|}_{\mathrm{large}\hbox{-} {N}_c\&{N}_f}\sim {\left(\frac{1}{x}\right)}^{\alpha_h^q}\cos \left[{\omega}_q\ln \left(\frac{1}{x}\right)+{\varphi}_q\right]. $$. The power $α^{q}_{h}$ exhibits a weak N f -dependence, and, for all Nf values considered, remains very close to $α^{q}_{h}$(N f =0)=$\left(4/\sqrt{3}\right)\sqrt{\alpha_s{N}_c/\left(2\pi \right)}$ obtained earlier in the large-N c limit. The novel oscillation frequency ω q and phase shift φ q depend more strongly on the number of flavors N f (with ω q = 0 in the pure-glue large-N c limit). The typical period of oscillations for ΔΣ is rather long, spanning many units of rapidity. We speculate whether the oscillations we find are related to the sign variation with x seen in the strange quark helicity distribution extracted from the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Biomass-Derived Polymeric Binders in Silicon Anodes for Battery Energy Storage Applications

The demand for portable electronic devices has increased rapidly during the past decade, and has driven a concordant growth in battery production. Since their development as a commercial energy storage solution in the 1990s, lithium-ion batteries (LIBs) have attracted significant attention in both science and industry due to their long cycle life, high energy density, low self-discharge rate, and high working voltage. Production of LIBs requires large amounts of a polymeric binder – commonly polyvinylidene difluoride (PVDF) – for processing and performance purposes. However, since this material is petrochemically-derived, it is far from “green” or sustainable. On the other hand, polymers and their building blocks are found widely throughout nature and can be renewably sourced from biomass at low cost; therefore, replacing PVDF with biomass-derived binders is a promising approach to reduce the environmental footprint of LIBs. Additionally, polymer binders play a critical role in next-generation battery performance. For instance, silicon (Si) is a promising high-capacity anode material for LIBs because of its high theoretical capacity (4200 mA h g –1 ), low working potential, and high abundance in Earth's crust. However, its huge volume change during charge/discharge tends to result in a shortened cycle life, since conventional binders interact only weakly with silicon's native surface and cannot maintain long-term integrity of the electrode. Naturally derived polymers have found better success in this role due to their high structural advantages. In this review, we summarize recent developments in silicon anode binders derived from various biomass sources, with a focus on polymer properties and their effect on battery performance. Further, we propose various perspectives based on our own assessment of these works, and provide a brief commentary on the future outlook of the field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Observation of Extraordinary Vibration Scatterings Induced by Strong Anharmonicity in Lead‐Free Halide Double Perovskites

Abstract Lead‐free halide double perovskites provide a promising solution for the long‐standing issues of lead‐containing halide perovskites, i.e., the toxicity of Pb and the low stability under ambient conditions and high‐intensity illumination. Their light‐to‐electricity or thermal‐to‐electricity conversion is strongly determined by the dynamics of the corresponding lattice vibrations. Here, the measurement of lattice dynamics is presented in a prototypical lead‐free halide double perovskite(Cs 2 NaInCl 6 ). The quantitative measurements and first‐principles calculations show that the scatterings among lattice vibrations at room temperature are at the timescale of ≈1 ps, which stems from the extraordinarily strong anharmonicity in Cs 2 NaInCl 6 . Further the degree of anharmonicity of each type of atom is quantitatively characterized in the Cs 2 NaInCl 6 single crystal, which stems from the interatomic forces, and demonstrate that this strong anharmonicity is synergistically contributed by the bond hierarchy, the tilting of the NaCl 6 and InCl 6 octahedral units, and the rattling of Cs + ions. Consequently, the crystalline Cs 2 NaInCl 6 possesses an ultralow thermal conductivity of ≈0.43 W mK −1 at room temperature, and a weak temperature dependence ofT −0.41 . These findings uncovered the underlying mechanisms behind the dynamics of lattice vibrations in double perovskites, which can largely benefit the design of optoelectronics and thermoelectrics based on halide double perovskites.

Chemistry↗

DOE OSTI2026

Toward Hydrogen Isotope Separations through Strong Hydrogen Adsorption at Open Copper(I) Sites in an Ultramicroporous Metal-Organic Framework

Metal-organic frameworks with coordinatively unsaturated metal sites (open metal sites) capable of engaging in orbital interactions with pi-acidic gases are of interest for enabling ambient-temperature gas separations, such as hydrogen isotope separations. In view of the weakly pi-acidic nature of H2, we sought to strengthen pi-backbonding-mediated H2 adsorption through pore confinement effects. Toward that end, we synthesized and characterized the ultramicroporous metal-organic framework CuxZn5-xCl4-yHz(bbta)3 (CuIZn-MFU-4; H2bbta = 1H,5H-benzo(1,2-d:4,5-d')bistriazole), featuring pi-basic trigonal pyramidal CuI sites that reside within 7 A of one another at their closest. Gas adsorption measurements reveal an H2 adsorption enthalpy of -38 kJ/mol, exceeding that of the larger-pore analog (CuIZn-MFU-4l; -33 kJ/mol) and representing the strongest H2 adsorption yet achieved in a metal-organic framework. The stronger H2 adsorption in CuIZn-MFU-4 is attributed to a combination of pore confinement effects and the increased ..sigma..-accepting nature of the CuI sites caused by a more electron-withdrawing bbta2- linker, as supported by structural, spectroscopic, and computational evidence. With the strongest H2 adsorption, equilibrium isotope effects in CuIZn-MFU-4 lead to a D2/H2 selectivity (as estimated by ideal adsorbed solution theory) of 1.35 even at 298 K, approaching the values reported below 200 K for conventional porous materials.

08 HYDROGEN↗

DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗

DOE OSTI2023

High-temperature ferromagnetism in Cr 1+x Pt 5-x P

In this work, we present the growth and basic magnetic and transport properties of Cr 1+x Pt 5-x P. We show that single crystals can readily be grown from a high-temperature solution created by adding dilute quantities of Cr to Pt-P based melts. Like other 1-5-1 compounds, Cr 1+x Pt 5-x P adopts a tetragonal P4/mmm, structure composed facesharing CrPts like slabs that are broken up along the c-axis by sheets of P atoms. EDS and X-ray diffraction measurements both suggest Cr 1+x Pt 5-x P has mixed occupancy between Cr and Pt atoms, similar to what is found in the closely related compound CrPtg, giving real compositions of Cr 1.5 Pt 4.5 P (x = 0.5). We report that Cr 1.5 Pt 4.5 P orders ferromagnetically at T c = 464.5 K with a saturated moment of ≈ 2.1 μs/Cr at 1.8 K. Likely owing to the strong spin-orbit coupling associated with the large quantity of high Z, Pt atoms, Cr 1.5 Pt 4.5 P has exceptionally strong planar anisotropy with estimated anisotropy fields of 345 kOe and 220 kOe at 1.8 K and 300 K respectively. The resistance of Cr 1.5 Pt 4.5 P has a metallic temperature dependence with relatively weak magnetoresistance. Electronic band structure calculations show that CrPt 5 P has a large peak in the density of states near the Fermi level which is split into spin majority and minority bands in the ferromagnetic state. Furthermore, the calculations suggest substantial hybridization between Cr-3d and Pt-5d states near the Fermi level, in agreement with the experimentally measured anisotropy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Data Centers and Digital Assurance Workshop 3 – Mitigations for Digital Assurance Risks

The third session of the TADA (Technical Assistance for Digital Assurance) Data Centers Cohort, held on November 18, 2025, focused on developing mitigation strategies for digital assurance risks identified in previous workshops. Hosted by Idaho National Laboratory (INL) and ScottMadden, the session emphasized the application of Cyber-Informed Engineering (CIE) to data center infrastructure, particularly at the utility–data center interface. Participants revisited and ranked key digital assurance risks, including architecture and interface weaknesses, governance gaps, and AI-enabled threats. The workshop introduced the 12 principles of CIE, advocating for consequence-focused design, engineered controls, and secure information architecture to proactively reduce cyber-physical vulnerabilities. These principles were applied to critical data center systems such as power distribution, UPS, cooling, SCADA/BMS, and grid-forming batteries. The session also addressed governance challenges at the interconnection boundary, highlighting the need for clear roles in telemetry sharing, firmware management, and trip settings. Special attention was given to emerging risks from behind-the-meter (BTM) generation, including reverse-power flow and the integration of small modular reactors (SMRs), which shift data centers from large loads to complex generation nodes. Participants explored how interconnection agreements can serve as enforceable instruments for digital assurance, and reviewed gaps in current standards such as NERC CIP, IEC 62443, and IEEE 1547. The workshop concluded with pathways to standardization, including model agreement language, state-level programs, and expanded NERC guidance. INL also presented tools and frameworks for secure procurement and supplier risk management, reinforcing the need for integrated engineering and policy solutions to secure the evolving data center–grid ecosystem. Session 3 of 3.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Unexpected Conformational Behavior of Poly(poly(ethylene glycol) methacrylate)-Poly(propylene carbonate)-Poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA) Amphiphilic Block Copolymers in Micellar Solution and at the Air-Water Interface

Self-assemblies (bulk micelles and Langmuir monolayers) of a new amphiphilic block copolymer, poly(poly(ethylene glycol) methacrylate)-poly(propylene carbonate)-poly(poly(ethylene glycol) methacrylate) (PPEGMA-PPC-PPEGMA), were investigated in dilute aqueous solution and at the air-water interface. An important feature that distinguishes PPEGMA-PPC-PPEGMA from conventional linear PEG-based polymer surfactants is that in the hydrophilic PPEGMA block the PEG moieties exist as side chains attached to the poly(methacrylate) (PMA) backbone. This study was performed using three PPEGMA-PPC-PPEGMA samples having an identical PPC molecular weight (5.6 kDa) and different PPEGMA molecular weights (7.3, 2.8 and 2.1 kDa on either side) (named as “G7C6G7”, “G3C5G3”, and “G2C6G2”, respectively). Cryo-TEM images revealed that G2C6G2 micelles form large clusters, whereas G7C6G7 micelles exist as unimers; the behavior of G3C6G3 is intermediate. The average core surface distance between two adjacent G2C6G2 micelles in a cluster was estimated to be only about 3.7 nm, which suggests that in the micelle corona layer only the PEG side chains (but not the PMA backbone) exist as hydrated brush chains. A calculation showed that at the small micelle separation distance observed, the van der Waals attractive forces become significant, and are thus likely responsible for the cluster formation. Steady shear rheometry measurements confirmed that G2C6G2 micelles are indeed only weakly bound in clusters, and these micelle clusters are thus prone to breakup even under mild steady shear. The micelle breakup force estimated from the shear viscosity data is in reasonable agreement with an estimate obtained from the theoretical van der Waals interaction potential. Langmuir monolayers formed at the air-water interface by PPEGMA-PPC-PPEGMA were also characterized; surface pressure-area measurements were performed. Detailed features of the surface pressure-area isotherms further supported that only the PEG side chains are hydrated in the subphase, while the PMA backbone remains unhydrated and situated on the water surface.

Lee, Jaewon↗

DOE OSTI2022

Charge Delocalization and Vibronic Couplings in Quadrupolar Squaraine Dyes

Squaraines are prototypical quadrupolar charge-transfer chromophores that have recently attracted much attention as building blocks for solution-processed photovoltaics, fluorescent probes with large two-photon absorption cross sections, and aggregates with large circular dichroism. Their optical properties are often rationalized in terms of phenomenological essential state models, considering the coupling of two zwitterionic excited states to a neutral ground state. As a result, optical transitions to the lowest S 1 excited state are one-photon allowed, whereas the next higher S 2 state can only be accessed by two-photon transitions. A further implication of these models is a substantial reduction of vibronic coupling to the ubiquitous high-frequency vinyl-stretching modes of organic materials. Here, in this study, we combine time-resolved vibrational spectroscopy, two-dimensional electronic spectroscopy, and quantum-chemical simulations to test and rationalize these predictions for nonaggregated molecules. We find small Huang–Rhys factors below 0.01 for the high-frequency, 1500 cm –1 modes in particular, as well as a noticeable reduction for those of lower frequency modes in general for the electronic S 0 → S 1 transition. The two-photon allowed state S 2 is well separated energetically from S 1 and has weak vibronic signatures as well. Thus, the resulting pronounced concentration of the oscillator strength in a narrow region relevant to the lowest electronic transition makes squaraines and their aggregates exceptionally interesting for strong and ultrastrong coupling of excitons to localized light modes in external resonators with chiral properties that can largely be controlled by the molecular architecture.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Experimental and computational kinetics study of the liquid-phase hydrogenation of C=C and C=O bonds

Solvent effects on adsorption equilibrium and reaction kinetics are evaluated for hydrogenation reactions catalyzed by Pd/alumina in a series of different solvents. Three reactants – cyclohexene, benzene and benzaldehyde – and three solvents – n-heptane (HEP), methylcyclohexane (MCHA), decalin (DL) – have been investigated. Kinetic analysis of hydrogenation of cyclohexene reveals that hydrogen adsorbs on different sites from those where cyclohexene and the solvents adsorb; however, the presence of hydrogen on these separate sites affects the heats of adsorption of the hydrocarbons. When the solvent is a weakly interacting linear alkane (HEP), the rate determining step of the reaction is the first hydrogenation of adsorbed cyclohexene. Furthermore, this conclusion is supported by DFT calculations that show a higher enthalpy barrier for the first hydrogenation than for the second, while statistical thermodynamics analysis validates the physical significance of the entropy of adsorption parameters derived from the kinetic fitting of experimental data. By contrast, with solvents such as MCHA and DL, which interact more strongly with the metal surface and compete for active sites with the reactant and the surface intermediate, the rate limiting step seems to shift to the second hydrogenation step.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗