SEARCH · Engineering Papers
Results for “Computational Complexity”
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.
Fundamental Insights into the Reactivity and Utilization of Open Metal Sites in Cu(I)-MFU-4 l
Metal organic frameworks (MOFs) having open metal sites have the potential to approach the activity of homogeneous organometallic complexes, thus combining the advantages of homogeneous and heterogeneous catalysis. We pre-sent a fundamental study of the effectiveness of incorporating open metal sites into MOFs. We have modeled the binding of a series of adsorbates in a Cu(I) substituted MOF, Cu(I)-MFU-4l using density functional theory, and compared the activity of the Cu(I) sites in Cu(I)-MFU-4l with two different Cu(I) scorpionate complexes. The computational results confirm the single-site nature of the Cu(I) active site. This is further supported by complementary experiments to measure the chemisorption uptake inside our synthesized samples in order to estimate the amount of active Cu sites present. We observed a level of chemisorption that is roughly half the theoretical maximum, which implies that only half of the Cu atoms incorporated into MFU-4l via metal-ion exchange are able to act as binding sites. We speculate that the inactive Cu atoms are coordinately saturated Cu(II) sites. Furthermore, our work suggests that the performance of Cu(I)-MFU-4l could be significantly increased by optimizing the metal exchange and activation processes.
MXene Nanosheets Functionalized with Cu Atoms for Urea Adsorption in Aqueous Media
Ti 3 C 2 T x MXene is an emerging family of two-dimensional materials, and because of its large specific surface area, it has potential for many applications. Herein, a new application using Cu-doped MILD (minimally intensive layer delamination) synthesized Ti 3 C 2 T x MXene for urea removal is demonstrated. The doping of Cu on MXene results in an increase in its affinity for urea adsorption as compared to the pristine MILD synthesized MXene due to the formation of the Cu–urea complex. Previous computational studies have shown that the adsorption energies of urea on the MXene surface can be improved in the presence of Cu. The valence state of Cu in the doped MILD synthesized MXene, which binds on to the surface via Ti–O–Cu linkage, is between 0 and +1 as verified by XAS and XPS. As the optimal urea adsorption occurs on Cu as a single atom site, an increase in Cu doping on MXene does not increase urea removal due to Cu agglomeration. Furthermore, looking at the adsorption behaviour, it seems that Cu-doped MXene follows the monolayer adsorption on homogenous surface model.
Nuclear quantum effects enter the mainstream
Atomistic simulations of chemical, biological and materials systems have become increasingly precise and predictive owing to the development of accurate and efficient techniques that describe the quantum mechanical behaviour of electrons. Nevertheless, the overwhelming majority of such simulations still assumes that the nuclei behave as classical particles. Historically, this approximation could sometimes be justified owing to the complexity and computational overhead. However, neglecting nuclear quantum effects has become one of the largest sources of error, especially when systems containing light atoms are treated using current state-of-the-art descriptions of chemical interactions. Over the past decade, this realization has spurred a series of methodological advances that have dramatically reduced the cost of including these important physical effects in the structure and dynamics of chemical systems. We discuss how these developments are now allowing nuclear quantum effects to become a mainstream feature of molecular simulations. These advances have led to new insights into phenomena that are relevant to different areas of science — from biochemistry to condensed matter — and open the door to many exciting future opportunities.
Let’s get our hands dirty: a comprehensive evaluation of DAQDB, key-value store for petascale hot storage
Data acquisition systems are a key component for successful data taking in any experiment. The DAQ is a complex distributed computing system and coordinates all operations, from the data selection stage of interesting events to storage elements. For the High Luminosity upgrade of the Large Hadron Collider, the experiments at CERN need to meet challenging requirements to record data with a much higher occupancy in the detectors. The DAQ system will receive and deliver data with a significantly increased trigger rate, one million events per second, and capacity, terabytes of data per second. An effective way to meet these requirements is to decouple real-time data acquisition from event selection. Data fragments can be temporarily stored in a large distributed key-value store. Fragments belonging to the same event can be then queried on demand, by the data selection processes. Implementing such a model relies on a proper combination of emerging technologies, such as persistent memory, NVMe SSDs, scalable networking, and data structures, as well as high performance, scalable software. In this paper, we present DAQDB (Data Acquisition Database) — an open source implementation of this design that was presented earlier, with an extensive evaluation of this approach, from the single node to the distributed performance. Furthermore, we complement our study with a description of the challenges faced and the lessons learned while integrating DAQDB with the existing software framework of the ATLAS experiment.
A moment-convergence method for stochastic analysis of biochemical reaction networks
Traditional moment-closure methods need to assume that high-order cumulants of a probability distribution approximate to zero. However, this strong assumption is not satisfied for many biochemical reaction networks. Here, we introduce convergent moments (defined in mathematics as the coefficients in the Taylor expansion of the probability-generating function at some point) to overcome this drawback of the moment-closure methods. As such, we develop a new analysis method for stochastic chemical kinetics. This method provides an accurate approximation for the master probability equation (MPE). In particular, the connection between low-order convergent moments and rate constants can be more easily derived in terms of explicit and analytical forms, allowing insights that would be difficult to obtain through direct simulation or manipulation of the MPE. In addition, it provides an accurate and efficient way to compute steady-state or transient probability distribution, avoiding the algorithmic difficulty associated with stiffness of the MPE due to large differences in sizes of rate constants. Applications of the method to several systems reveal nontrivial stochastic mechanisms of gene expression dynamics, e.g., intrinsic fluctuations can induce transient bimodality and amplify transient signals, and slow switching between promoter states can increase fluctuations in spatially heterogeneous signals. The overall approach has broad applications in modeling, analysis, and computation of complex biochemical networks with intrinsic noise.
Plasma transport in an Eulerian AMR code
A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions to flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unit of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.
Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow
A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.
Opinion dynamics in financial markets via random networks
We investigate financial market dynamics by introducing a heterogeneous agent-based opinion formation model. In this work, we organize individuals in a financial market according to their trading strategy, namely, whether they are noise traders or fundamentalists. The opinion of a local majority compels the market exchanging behavior of noise traders, whereas the global behavior of the market influences the decisions of fundamentalist agents. We introduce a noise parameter, q , to represent the level of anxiety and perceived uncertainty regarding market behavior, enabling the possibility of adrift financial action. We place individuals as nodes in an Erdös-Rényi random graph, where the links represent their social interactions. At any given time, individuals assume one of two possible opinion states ±1 regarding buying or selling an asset. The model exhibits fundamental qualitative and quantitative real-world market features such as the distribution of logarithmic returns with fat tails, clustered volatility, and the long-term correlation of returns. We use Student’s t distributions to fit the histograms of logarithmic returns, showing a gradual shift from a leptokurtic to a mesokurtic regime depending on the fraction of fundamentalist agents. Furthermore, we compare our results with those concerning the distribution of the logarithmic returns of several real-world financial indices.
Tritium Breeding Ratio Evaluation of Solid Breeder Concepts for the FESS-FNSF
This paper presents a parametric study of the Fusion Energy System Studies-Fusion Nuclear Science Facility’s (FNSF’s) tritium breeding performance for several solid breeder concepts, neutron multiplying materials, and blanket materials, assuming volume fractions based on the most recent FNSF design as a realistically representative fusion facility. In this study, we initially surveyed the tritium breeding ratio (TBR) of several solid breeder concepts by employing a simplified but efficient one-dimensional (1-D) infinite cylinder reduced-order model (ROM). Parametric studies were performed with the ROMs for the full range of breeder-to-multiplier ratios to identify the optimum mixture compositions for each breeder type that would lead to a maximum TBR. These optimized breeder-multiplier combinations were then homogenized with FNSF blanket component materials to estimate their impacts on the TBR. Subsequently, as a validation step for the optimal designs, TBR calculations were performed using a more realistic modified 1-D ROM with inner and outer breeding regions, as well as with a fully detailed 22.5-deg three-dimensional (3-D) sector of the FNSF to assess the impact of geometry details on the TBR. The differences between the two 1-D models were negligible, while the ROMs were able to correctly predict trends and identify the maximum and minimum TBR cases, as well as show consistent biases relative to the results produced by the full 3-D, 22.5-deg sector for specific breeder/multiplier combinations. Solid breeder concepts such as Li 2 O, Li 4 SiO 4 , and Li 8 ZrO 6 outperformed all others in this study in terms of TBR performance when combined with all the neutron multiplier materials selected. Here, an underlying goal of this study was to develop and improve rapid and reliable ROMs to aid designers during parametric optimizations of highly complex and computationally expensive fusion models.
SALSAA: a statistical approach to line shapes from an average atom
Ion-Stark line broadening is a key density diagnostic for hot dense plasmas relevant to inertial fusion and astrophysics. It is caused by interactions of a radiating ion with nearby perturbing ions, whose electric microfields lead to changes in bound-bound transition energies. Ion-Stark broadening becomes increasingly difficult to compute for complex, many-electron ions with myriad transitions. In this paper, we propose a simplified approach to ion-Stark broadening based on self-consistent ion distributions and electronic structure from an average-atom model. We find that this approach reproduces the line shape predictions of one traditional method for high-n K-shell emission lines in aluminum ions with accuracy sufficient for density diagnostics in thermal plasmas with equal ion and electron temperatures. We expect that this approach can be extended to provide a reasonable picture of ion-Stark broadening in many-electron ions, enabling rapid calculations of line profiles in complex spectra.
State-to-state rate coefficients for HCS+ in rotationally inelastic collisions with H2 at low temperatures
ABSTRACT HCS+ ions have been detected in several regions of the interstellar medium (ISM), but an accurate determination of the chemical-physical conditions in the molecular clouds where this molecule is observed requires detailed knowledge of the collisional rate coefficients with the most common colliders in those environments. In this work, we study the dynamics of rotationally inelastic collisions of HCS+ + H2 at low temperature, and report, for the first time, a set of rate coefficients for this system. We used a recently developed potential energy surface for the HCS+–H2 van der Waals complex and computed state-to-state rotational rate coefficients for the lower rotational states of HCS+ in collision with both para- and ortho-H2, analysing the influence of the computed rate coefficients on the determination of critical densities. Additionally, the computed rate coefficients are compared with those obtained by scaling the ones from HCS+ in collision with He (an approximation that is sometimes used when data is lacking), and large differences are found. Furthermore, the approximation of using the rates for the HCO+ + H2 collision as a rough approximation for those of the HCS+ + H2 system is also evaluated. Finally, the complete set of de-excitation rate coefficients for the lowest 30 rotational states of HCS+ by collision with H2 is reported from 5 to 100 K.
Exploratory calculation of 𝐾 L → 𝜇 + 𝜇 − decay from lattice QCD at physical pion mass
We compute the complex, long-distance two-photon-exchange amplitude which contributes to the rare 𝐾 L → 𝜇 + 𝜇 − decay from lattice QCD. We use a 24 3 × 64 physical-pion-mass gauge field ensemble at an inverse lattice spacing of 1.023 GeV and a QED ∞ -based formalism. Our implementation strategies for all five non-SU(3)-flavor-suppressed diagram topologies are given in detail. We achieve a 25% statistical precision on the dispersive part of this long-distance amplitude. This calculation is carried out with 2+1 quark flavors and therefore requires the addition of counterterms to compensate for the absence of the Glashow-Iliopoulos-Maiani mechanism. These counterterms are not included in the current calculation and will be the subject of a second paper. Although a direct comparison to experiment cannot yet be made because of those omitted counterterms, the present exploratory calculation allows one to identify principal sources of statistical uncertainty in this calculation. The precision of our results is limited by the reconstruction of the physical contribution of the 𝜂 intermediate state, for which various strategies are tested and compared.
Lieb-Robinson Bounds with Exponential-in-Volume Tails
Lieb-Robinson bounds demonstrate the emergence of locality in many-body quantum systems. Intuitively, Lieb-Robinson bounds state that, with local or exponentially decaying interactions, the correlation that can be built up between two sites separated by distance 𝑟 after a time 𝑡 decays as exp (𝑣𝑡 −𝑟), where 𝑣 is the emergent Lieb-Robinson velocity. In many problems, it is important to also capture how much of an operator grows to act on 𝑟 𝑑 sites in 𝑑 spatial dimensions. Perturbation theory and cluster expansion methods suggest that, at short times, these volume-filling operators are suppressed as exp (−𝑟 𝑑 ). We confirm this intuition, showing that, for 𝑟 >𝑣𝑡, the volume-filling operator is suppressed by exp (−(𝑟−𝑣𝑡) 𝑑 /(𝑣𝑡) 𝑑−1 ). This closes a conceptual and practical gap between the cluster expansion and the Lieb-Robinson bound. We then present two very different applications of this new bound. Firstly, we obtain improved bounds on the classical computational resources necessary to simulate many-body dynamics with error tolerance 𝜀 for any finite time 𝑡: as 𝜀 becomes sufficiently small, only 𝜀 −O(𝑡 𝑑−1 ) resources are needed. A protocol that likely saturates this bound is given. Secondly, we prove that disorder operators have volume-law suppression near the “solvable (Ising) point” in quantum phases with spontaneous symmetry breaking, which implies a new diagnostic for distinguishing many-body phases of quantum matter.
Reinforcement learning based schemes to manage client activities in large distributed control systems
Not Available
Simulation based inference for efficient theory space sampling: An application to supersymmetric explanations of the anomalous muon g − 2
Not Available
Intractability of Electronic Structure in a Fixed Basis
Not Available
Connecting relativistic density functional theory to microscopic calculations
The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.