Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Exact two component”

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 OSTI2020

Implementing Chiral Three-Body Forces in Terms of Medium-Dependent Two-Body Forces

Three-nucleon (3N) forces are an indispensable ingredient for accurate few-body and many- body nuclear structure and reaction theory calculations. While the direct implementation of chiral 3N forces can be technically very challenging, a simpler approach is given by employing instead a medium-dependent NN interaction V med that reflects the physics of three-body forces at the two-body normal-ordered approximation. We review the derivation and construction of Vmed from the chiral 3N interaction at next-to-next-to-leading order (N2LO), consisting of a long-range 2π-exchange term, a mid-range 1π-exchange component and a short-range contact-term. Several applications of V med to the equation of state of cold nuclear and neutron matter, the nucleon single-particle potential in nuclear matter, and the nuclear quasiparticle interaction are discussed. We also explore differences in using local vs. nonlocal regulating functions on 3N forces and make direct comparisons to exact results at low order in perturbation theory expansions for the equation of state and single-particle potential. We end with a discussion and numerical calculation of the in-medium NN potential V med from the next-to-next-to-next-to-leading order (N3LO) chiral 3N force, which consists of a series of long-range and short-range terms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The density distribution and physical origins of intermittency in supersonic, highly magnetized turbulence with diverse modes of driving

The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

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

Liquid lithium as divertor material to mitigate severe damage of nearby components during plasma transients

The successful operation of thermonuclear fusion reactors such as ITER, DEMO, and future commercial plants is mainly determined by the optimum choice of materials for various components. The objective of this work is to accurately and comprehensively simulate the entire device in 3D to predict pros and cons of various materials, e.g., liquid lithium in comparison to tungsten and carbon to predict future ITER-like and DEMO divertor performances. We used our comprehensive HEIGHTS simulation package to investigate ITER-like components response during transient events in exact 3D geometry. Starting from the lost hot core plasma particles through SOL, deposition on the divertor surface, and the generation of secondary plasma of divertor materials. Our simulations predicted significant reduction in the heat loading and damage to the divertor nearby and internal components in the case when lithium is used on the divertor plates. While if tungsten or carbon are used on the divertor plate, significant melting areas and vaporization spots can occur (less for carbon) on the reflector, dome, and stainless steel tubes, and even parts of the first walls can melt due to the high radiation power of the secondary divertor plasma. Lithium photon radiation deposition into the divertor and nearby surfaces was decreased by two orders of magnitude compared to tungsten and by one order of magnitude compared to carbon. This analysis showed that using liquid lithium for ITER-like surfaces and future DEMO can lead to significant enhancement in components lifetime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Light cone dynamics in excitonic states of two-component Bose and Fermi gases

Herein we consider the non-equilibrium dynamics of two-component one dimensional quantum gases in the limit of extreme population imbalance where the minority species has but a single particle. We consider the situation where the gas is prepared in a state with a single spatially localized exciton: the single particle of the minority species is spatially localized while the density of the majority species in the vicinity of the minority particle sees a depression. Remarkably, we are able to consider cases where the gas contains on the order of $N$ = 100 particles, comparable to that studied in experiments on cold atomic gases. We are able to do so by exploiting the integrability of the gas together with the observation that the excitonic state can be constructed through a simple superposition of exact eigenstates of the gas. The number of states in this superposition, rather than being exponentially large in the number of particles, scales linearly with $N$. We study the evolution of such spatially localized states in both strongly interacting Bose and Fermi gases. The behavior of the light cones when the interaction strength and density of the gas is varied can be understood from exact results for the spin excitation spectrum in these systems. We argue that the light cone in both cases exhibits scaling collapse. However unique to the Bose gas, we show that the presence of gapped finite-momentum roton-like excitations provide the Bose gas dynamics with secondary light cones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A semiempirical transparency model for dual energy cargo radiography applications

Cargo containers passing through ports are scanned by non-intrusive inspection systems to search for concealed illicit materials. By using two photon beams with different energy spectra, dual energy inspection systems are sensitive to both the area density and the atomic number of cargo contents. Most literature on the subject assumes a simple exponential attenuation model for photon intensity in which only free streaming photons are detected. However, this approximation neglects second order effects such as scattering, leading to a biased model and thus incorrect material predictions. This work studies the accuracy of the free streaming model by comparing it to simulation outputs, finding that the model shows poor atomic number reconstruction accuracy at high- Z and suffers significantly if the source energy spectra and detector response function are not known exactly. To address these challenges, this work introduces a semiempirical transparency model which modifies the free streaming model by rescaling different components of the mass attenuation coefficient, allowing the model to capture secondary effects ignored by the free streaming model. The semiempirical model displays improvement agreement with simulated results at high- Z and shows excellent extrapolation to materials and thicknesses which were not included during the calibration step. Furthermore, this work demonstrates that the semiempirical model yields accurate atomic number predictions even when the source spectra and detector response are not known exactly. In conclusion, using the semiempirical model, manufacturers can perform a simple calibration to enable more precise Z reconstruction capabilities, which has the potential to significantly improve the performance of existing dual energy radiographic systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Mechanistic Studies of a Primitive Homolog of Nitrogenase Involved in Coenzyme F430 Biosynthesis

Methyl-coenzyme M reductase (MCR) is the key enzyme in the biological formation and anaerobic oxidation of methane (AOM). Methane is a potent greenhouse gas and the major component of natural gas. Given the abundance of natural gas reserves in remote areas, there is great current interest in a scalable bio-based process for the conversion of methane to liquid fuel or other high-value commodity chemicals. MCR holds much promise for use in such a methane bioconversion strategy. However, MCR cannot currently be produced in an active form in a heterologous host, due in large part to the lack of genetic and biochemical information about the production of holo MCR. In an effort to overcome this deficiency, our laboratory elucidated the biosynthetic pathway of the unique nickel-containing coenzyme of MCR, F430. The key step in coenzyme F430 biosynthesis (Cfb) was found to involve an unprecedented reductive cyclization reaction that converts Ni-sirohydrochlorin a , c -diamide to 15,17 3 -seco-F430-17 3 -acid. This remarkable transformation, which involves a 6-electron reduction of the isobacteriochlorin ring system, cyclization of the c -acetamide side chain to form a γ-lactam ring, and the formation of 7 stereocenters, is catalyzed by a primitive homolog of nitrogenase (CfbCD). Nitrogenase is a two-component metalloenzyme that catalyzes the ATP-dependent reduction of dinitrogen to ammonia (nitrogen fixation). Homologs of nitrogenase are also involved in the biosynthesis of the photosynthetic pigments chlorophyll and bacteriochlorophyll. Phylogenetic analysis of the CfbCD complex suggests that it is representative of a more ancient lineage of the nitrogenase superfamily, and a thorough investigation of its structure and function is likely to shed light on the mechanisms and evolution of these important metalloenzymes that catalyze multi-electron redox reactions. Moreover, a detailed understanding of the mechanism of the CfbCD complex may aid in the development of specific inhibitors to help reduce natural greenhouse gas emissions and can be exploited for the heterologous production of MCR for methane bioconversion. Towards these goals, the following Specific Aims will be pursued to determine the: 1) Identity of the CfbCD reaction product. The exact reaction catalyzed by CfbCD, including the number of electrons transferred and whether it involves enzymatic cyclization, is unclear. Several approaches, including reaction stoichiometry measurements, spectroelectrochemistry, and magnetic resonance spectroscopy will be applied to elucidate the structure of the reaction product and establish whether CfbCD is a reductase or reductive cyclase. 2) Structure, conformational dynamics, and oligomerization state changes of CfbCD. Significant insight into the mechanism and allosteric regulation of CfbCD can be obtained by assessing changes in the structure and dynamics of the complex during the catalytic cycle. To accomplish this, a combination of size-exclusion chromatography, hydrogen-deuterium exchange mass spectrometry, molecular dynamics simulations, and high-resolution structural methods will be employed. 3) Source, order, and stereochemistry of proton additions during CfbCD catalysis. Details regarding the order and stereochemistry of proton additions during the CfbCD reaction will be uncovered using a combined spectroscopic and computational approach. Complementary mechanistic studies employing site-directed mutagenesis and substrate analogs will establish the identity of active site acid residues and the possible involvement of substrate-assisted catalysis during the CfbCD reaction.

09 BIOMASS FUELS↗
DOE OSTI2025

Evaluating an Earth system model from a water manager perspective

Abstract. The large spatial scale of global Earth system models (ESMs) is often cited as an obstacle to using the output by water resource managers in localized decisions. Recent advances in computing have improved the fidelity of hydrological responses in ESMs through increased connectivity between model components. However, the models are seldom evaluated for their ability to reproduce metrics that are important for and resonate with practitioners or that allow practitioners to situate higher-resolution model outputs within a cascade of uncertainty stemming from different models and scenarios. We draw on the combined experience of the author team and water manager workshop participants to identify salient water management metrics and evaluate whether they are credibly reproduced over the conterminous USA by the Community Earth System Model v2 (CESM2) Large Ensemble. We find that, while the exact values may not match the observations, aspects such as interannual variability can be reproduced by CESM2 for the mean wet day precipitation and length of dry spells. CESM2 also captures the proportion of total annual precipitation that derives from the heaviest rain days in watersheds that are not snow-dominated. Aggregating the 7 d mean daily runoff to two-digit Hydrological Unit Code (HUC2) watersheds also shows that rain-dominated regions capture the timing and interannual variability of annual maximum and minimum flows. We conclude that there is potential for far greater use of large-ensemble ESMs, such as CESM2, in long-range water management decisions to supplement high-resolution regional projections.

Tye, Mari R. (ORCID:0000000324911020)↗
DOE OSTI2023

Coupled Cluster Study of the Heats of Formation of UF 6 and the Uranium Oxyhalides, UO 2 X 2 (X = F, Cl, Br, I, and At)

The atomization enthalpies of the U(VI) species UF6 and the uranium oxyhalides UO 2 X 2 (X=F, Cl, Br, I, At) were calculated using a composite relativistic Feller-Peterson-Dixon (FPD) approach based on scalar relativistic DKH3-CCSD(T) with extrapolations to the CBS limit. The inherent multideterminantal nature of the U atom was mitigated by utilizing the singly charged atomic cation in all calculations with correction back to the neutral asymptote via the accurate ionization energy of the U atom. The effects of SO coupling were recovered using full 4-component CCSD(T) with contributions due to the Gaunt Hamiltonian calculated using Dirac-Hartree-Fock. The final atomization enthalpy for UF 6 (752.2 kcal/mol) was within 2.5 kcal/mol of the experimental value, but unfortunately the latter carries a ±2.4 kcal/mol uncertainty that is predominantly due to the experimental uncertainty in the formation enthalpy of U atom. The analogous value for UO 2 F 2 (607.6 kcal/mol) was in nearly exact agreement with experiment, but the latter has a stated experimental uncertainty of ±4.3 kcal/mol. The FPD atomization enthalpy for UO 2 Cl 2 (540.4 kcal/mol) was within the experimental error limits of ±5.5 kcal/mol. FPD atomization energies for the non-U-containing molecules (used for reaction enthalpies) H 2 O and HX (X=F, Cl, Br, I, At) were within at most 0.3 kcal/mol of their experimental values where available. The FPD atomization enthalpies, together with FPD reaction enthalpies for two different reactions, were used to determine heats of formation for all species of this work with estimated uncertainties of ±4 kcal/mol. The calculated heat of formation for UF 6 (-511.0 kcal/mol) is within 2.5 kcal/mol of the accurately-known (±0.45 kcal/mol) experimental value.

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

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Visualizing Comparisons of Bill of Materials

Protecting critical infrastructure from cyber attacks, natural disasters, and other disruptions is a priority of the U.S. Government. Critical infrastructure includes providing electricity to homes and businesses, supplying natural gas for heating, and producing renewable energy sources. A loss of these services, as seen in the Solarwinds supply chain attack in 2020 , Texas snowstorm of 2021, the Colonial Pipeline cyber incident of 2021, and the Washington power substation attacks in 2022 result in high costs to consumers, disruption of everyday life, and even death. To protect the infrastructure, we first have to know what equipment we are protecting. The complexity of distributed manufacturing and development coupled with the increasing prevalence of cyber and supply chain attacks necessitates a greater understanding of the hardware and software components that comprise equipment in critical infrastructure. When a vulnerability in a single software library can have disastrous consequences, it is vital to understand critical equipment and systems at a granular level. This need has led to increased energy around the development and incorporation of bill-of-materials (BOM) into existing asset management practices to aid in mitigating, and responding to future attacks \cite{noauthor_software_nodate}. While much of the current research is devoted to creating BOMs, it is equally important to develop methodologies for leveraging BOMs to answer questions, such as: How has my software changed? Are two pieces of equipment equivalent? Does this piece of equipment that just arrived match my historical information? In this work, we demonstrate how BOMs can be represented by graph structures. We then describe how these structures can be fed into a graph comparison algorithm to produce a novel interactive visualization that allows us to not only identify differences in BOMs, but show exactly where they are in the product.

Jones, Rebecca D.↗
DOE OSTI2021

Stability of the particle-hole Pfaffian state and the 5 2 -fractional quantum Hall effect

Here we present a method for the exact construction of the fully particle-hole symmetric Pfaffian (PH-Pfaffian) ground state and its charged excitations on a sphere. We adopt the Moore-Read state, but with a nonholomorphic pairing component as in previous studies, and project it to the lowest Landau level. We study the energetics as well as other properties of these states and find that in a pure system interacting with the Coulomb forces the PH-Pfaffian cannot compete with either the Moore-Read state or its particle-hole conjugate, the anti-Pfaffian state, as an explanation for the 5/2 effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Nontrivial fusion of Majorana zero modes in interacting quantum-dot arrays

Motivated by recent experimental reports of Majorana zero modes (MZMs) in quantum-dot systems at the “sweet spot,” where the electronic hopping t ℎ is equal to the superconducting coupling Δ, we study the time-dependent spectroscopy corresponding to the nontrivial fusion of MZMs. The term “nontrivial” refers to the fusion of Majoranas from different original pairs of MZMs, each with well-defined parities. We employ an experimentally accessible time-dependent real-space local density-of-states (LDOS) method to investigate the nontrivial MZM fusion outcomes in canonical chains and in a Y-shaped array of interacting electrons. In the case of quantum-dot chains where two pairs of MZMs are initially disconnected, after fusion we find equal-height peaks in the electron and hole components of the LDOS, signaling nontrivial fusion into both the vacuum I and fermion Ψ channels with equal weight. For π-junction quantum-dot chains, where the superconducting phase has opposite signs on the left and right portions of the chain, after the nontrivial fusion we observed the formation of an exotic two-site MZM near the center of the chain, coexisting with another single-site MZM. Furthermore, we also studied the fusion of three MZMs in the Y-shaped geometry. In this case, after the fusion we observed the novel formation of another exotic multisite MZM, with properties depending on the connection and geometry of the central region of the Y-shaped quantum-dot array.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Reionization history constraints from neural network based predictions of high-redshift quasar continua

ABSTRACT Observations of the early Universe suggest that reionization was complete by z ∼ 6, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Lyα damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Lyα, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening intergalactic medium (IGM). In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Lyα. Our Quasar Spectra from Artificial Neural Network based predictive Regression Algorithm(QSANNdRA) improves the error in this reconstruction compared to the state-of-the-art principal component analysis (PCA) based model in the literature by 14.2 per cent on average, and provides an improvement of 6.1 per cent on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1 per cent and 16.8 per cent when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at z = 7.0851 and ULAS J1342+0928 at z = 7.5413, respectively. Using our more accurate reconstructions of these two z > 7 quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{H\, \small{I}} = 0.25^{+0.05}_{-0.05}$ at z = 7.0851 and $\bar{x}_\mathrm{H\, \small{I}} = 0.60^{+0.11}_{-0.11}$ at z = 7.5413. Our results are consistent with the literature and favour a rapid end to reionization.

Ďurovčíková, Dominika↗
DOE OSTI2020

Tilted outer and inner structures in edge-on galaxies?

ABSTRACT Tilted and warped discs inside tilted dark matter haloes are predicted from numerical and semi-analytical studies. In this paper, we use deep imaging to demonstrate the likely existence of tilted outer structures in real galaxies. We consider two SB0 edge-on galaxies, NGC 4469 and NGC 4452, which exhibit apparent tilted outer discs with respect to the inner structure. In NGC 4469, this structure has a boxy shape, inclined by ΔPA ≈ 3° with respect to the inner disc, whereas NGC 4452 harbours a discy outer structure with ΔPA ≈ 6°. In spite of the different shapes, both structures have surface brightness profiles close to exponential and make a large contribution (∼30 per cent) to the total galaxy luminosity. In the case of NGC 4452, we propose that its tilted disc likely originates from a former fast tidal encounter (probably with IC 3381). For NGC 4469, a plausible explanation may also be galaxy harassment, which resulted in a tilted or even a tumbling dark matter halo. A less likely possibility is accretion of gas-rich satellites several Gyr ago. New deep observations may potentially reveal more such galaxies with tilted outer structures, especially in clusters. We also consider galaxies, mentioned in the literature, where a central component (a bar or a bulge) is tilted with respect to the stellar disc. According to our numerical simulations, one of the plausible explanations of such observed ‘tilts’ of the bulge/bar is a projection effect due to a not exactly edge-on orientation of the galaxy coupled with a skew angle of the triaxial bulge/bar.

Mosenkov, Aleksandr V.↗
DOE OSTI2025

Hetero-Orbital Two-Component Fractional Quantum Hall States in Bilayer Graphene

A two-dimensional electron system exposed to a strong magnetic field produces a plethora of strongly interacting fractional quantum Hall (FQH) states, the complex topological orders of which are revealed through exotic emergent particles, such as composite fermions, and fractionally charged Abelian and non-Abelian anyons. Much insight has been gained by the study of multicomponent FQH states, where spin and pseudospin indices of the electron contribute additional correlation. Traditional multicomponent FQH states develop in situations where the components share the same orbital states and the resulting interactions are pseudospin independent; this homo-orbital nature is also crucial to their theoretical understanding. Here, we study “hetero-orbital” two-component FQH states, in which the orbital index is part of the pseudospin, rendering the multicomponent interactions strongly SU(2) anisotropic in the pseudospin space. Such states, obtained in bilayer graphene at the isospin transition between 𝑁 = 0 and 𝑁 = 1 electron Landau levels, are markedly different from previous homo-orbital two-component FQH states. In particular, we observe strikingly different behaviors for the parallel-vortex and reverse-vortex attachment composite fermion states, and an anomalously strong two-component 2/5 state over a wide range of magnetic field before it abruptly disappears at a high field. Our findings, combined with detailed theoretical calculations, reveal the surprising robustness of the hetero-orbital FQH effects, significantly enriching our understanding of FQH physics in this novel regime.

Bilayer graphene↗
DOE OSTI2021

Explicit physics-informed neural networks for nonlinear closure: The case of transport in tissues

In upscaling methods, closures for nonlinear problems present a well-known challenge. While a number of theoretical methods have been proposed for handling such closures, nonlinearities still remain a significant obstacle for many problems. In this work, we use a combination of formal upscaling and data-driven machine learning for explicitly closing a nonlinear transport and reaction process in multiscale tissues. The classical effectiveness factor model is used to formulate the macroscale reaction kinetics. We train a multilayer perceptron network using training data generated by direct numerical simulations over microscale examples. Once trained, the network is used in an algorithm for numerically solving the upscaled (coarse-grained) differential equation describing mass transport and reaction in two example tissues. The network is described as being explicit in the sense that the network is trained using macroscale concentrations and gradients of concentration as components of the feature space rather than incorporating them as part of a constraint in the optimization process. Network training and solutions to the macroscale transport equations were computed for two different tissues. The two tissue types (brain and liver) exhibit markedly different geometrical complexity and spatial scale (cell size and sample size). The upscaled solutions for the average concentration are compared with numerical solutions derived from the microscale concentration fields by a posteriori averaging. There are three outcomes of this work of particular note. 1) Our overall approach results in an upscaled nonlinear PDE. The PDE is closed using a neural network, and our approach results in the definition of the classical effectiveness factor for effecting closure. 2) We identify particular source terms for the closure problem that are important for representing the structure of the closure. These source terms involve macroscale concentrations and their gradients. We adopt these source terms to use as explicit features in the learning algorithm. We find the trained networks that include the macroscale source terms generate models that are able to predict the correction factor with increased fidelity over those that do not. 3) We find that the trained network exhibits good generalizability, and it is able to predict the effectiveness factor with high fidelity for realistically-structured tissues despite the significantly different scale and geometrical complexity of the two example tissue types. This latter result emphasizes our purposeful connection between conventional averaging methods with the use of machine learning for closure; this contrasts with some machine learning methods for upscaling where the exact form of the macroscale equation remains unknown.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Understanding the effect of density functional choice and van der Waals treatment on predicting the binding configuration, loading, and stability of amine-grafted metal organic frameworks

Metal organic frameworks (MOFs) are crystalline, three-dimensional structures with high surface areas and tunable porosities. Made from metal nodes connected by organic linkers, the exact properties of a given MOF are determined by node and linker choice. MOFs hold promise for numerous applications, including gas capture and storage. M2(4,4′-dioxidobiphenyl-3,3′-dicarboxylate)—henceforth simply M2(dobpdc), with M = Mg, Mn, Fe, Co, Ni, Cu, or Zn—is regarded as one of the most promising structures for CO2 capture applications. Further modification of the MOF with diamines or tetramines can significantly boost gas species selectivity, a necessity for the ultra-dilute CO2 concentrations in the direct-air capture of CO2. There are countless potential diamines and tetramines, paving the way for a vast number of potential sorbents to be probed for CO2 adsorption properties. The number of amines and their configuration in the MOF pore are key drivers of CO2 adsorption capacity and kinetics, and so a validation of computational prediction of these quantities is required to suitably use computational methods in the discovery and screening of amine-functionalized sorbents. In this work, we study the predictive accuracy of density functional theory and related calculations on amine loading and configuration for one diamine and two tetramines. In particular, we explore the Perdew–Burke–Ernzerhof (PBE) functional and its formulation for solids (PBEsol) with and without the Grimme-D2 and Grimme-D3 pairwise corrections (PBE+D2/3 and PBEsol+D2/3), two revised PBE functionals with the Grimme-D2 and Grimme-D3 pairwise corrections (RPBE+D2/3 and revPBE+D2/3), and the nonlocal van der Waals correlation (vdW-DF2) functional. We also investigate a universal graph deep learning interatomic potential’s (M3GNet) predictive accuracy for loading and configuration. These results allow us to identify a useful screening procedure for configuration prediction that has a coarse component for quick evaluation and a higher accuracy component for detailed analysis. Our general observation is that the neural network-based potential can be used as a high-level and rapid screening tool, whereas PBEsol+D3 gives a completely qualitatively predictive picture across all systems studied, and can thus be used for high accuracy motif predictions. We close by briefly exploring the predictions of relative thermal stability for the different functionals and dispersion corrections.

Chemistry↗