Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster computing”

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 829 records · Page 46

Projected coupled cluster theory: Optimization of cluster amplitudes in the presence of symmetry projection

Methods which aim at universal applicability must be able to describe both weak and strong electronic correlation with equal facility. Such methods are in short supply. The combination of symmetry projection for strong correlation and coupled cluster theory for weak correlation offers tantalizing promise to account for both on an equal footing. In order to do so, however, the coupled cluster portion of the wave function must be optimized in the presence of the symmetry projection. Finally, this paper discusses how this may be accomplished, and shows the importance of doing so for both the Hubbard model Hamiltonian and the molecular Hamiltonian, all with a computational scaling comparable to that of traditional coupled cluster theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sci-Thur AM: YIS – 06: A Monte Carlo study of macro- and microscopic dose descriptors and the microdosimetric spread using detailed cellular models

Purpose: To develop Monte Carlo models of cell clusters to investigate the relationships between macro- and microscopic dose descriptors, quantify the microdosimetric spread in energy deposition for subcellular targets, and determine how these results depend on the computational model. Methods: Microscopic tissue structure is modelled as clusters of 13 to 150 cells, with cell (nuclear) radii between 5 and 10 microns (2 and 9 microns). Energy imparted per unit mass (specific energy or dose) is scored in the nucleus (D{sub nuc}) and cytoplasm (D{sub cyt}) for incident photon energies from 20 to 370 keV. Dose-to-water (D{sub w,m}) and dose-to-medium (D{sub m,m}) are compared to D{sub nuc} and D{sub cyt}. Single cells and single nuclear cavities are also simulated. Results: D{sub nuc} and D{sub cyt} are sensitive to the surrounding environment with deviations of up to 13% for a single nucleus/cell compared with a multicellular cluster. These dose descriptors vary with cell and nucleus size by up to 10%. D{sub nuc} and D{sub cyt} differ from D{sub w,m} and D{sub m,m} by up to 32%. The microdosimetric spread is sensitive to whether cells are arranged randomly or in a hexagonal lattice, and whether subcellular compartment sizes are sampled from a normal distribution or are constant throughout the cluster. Conclusions: D{sub nuc} and D{sub cyt} are sensitive to cell morphology, elemental composition and the presence of surrounding cells. The microdosimetric spread was investigated using realistic elemental compositions for the nucleus and cytoplasm, and depends strongly on subcellular compartment size, source energy and dose.

60 APPLIED LIFE SCIENCES↗

ARENA: Asynchronous Reconfigurable Accelerator Ring to Enable Data-Centric Parallel Computing

The next generation HPC and data centers are likely to be reconfigurable and data-centric due to the trend of hardware specialization and the emergence of data-driven applications. In this work, we propose ARENA – an asynchronous reconfigurable accelerator ring architecture as a potential scenario on how the future HPC and data centers will be like. Despite using the coarse-grained reconfigurable arrays (CGRAs) as the substrate platform, our key contribution is not only the CGRA-cluster design itself, but also the ensemble of a new architecture and programming model that enables asynchronous tasking across a cluster of reconfigurable nodes, so as to bring specialized computation to the data rather than the reverse. We presume distributed data storage without asserting any prior knowledge on the data distribution. Hardware specialization occurs at runtime when a task finds the majority of data it requires are available at the present node. In other words, we dynamically generate specialized CGRA accelerators where the data reside. The asynchronous tasking for bringing computation to data is achieved by circulating the task token, which describes the dataflow graphs to be executed for a task, among the CGRA cluster connected by a fast ring network. Evaluations on a set of HPC and data-driven applications across different domains show that ARENA can provide better parallel scalability with reduced data movement (53.9 percent). Compared with contemporary compute-centric parallel models, ARENA can bring on average 4.37× speedup. The synthesized CGRAs and their task-dispatchers only occupy 2.93mm 2 chip area under 45nm process technology and can run at 800MHz with on average 759.8mW power consumption. ARENA also supports the concurrent execution of multi-applications, offering ideal architectural support for future high-performance parallel computing and data analytics systems.

97 MATHEMATICS AND COMPUTING↗

Modeling mesoscale fission gas behavior in UO2 by directly coupling the phase field method to spatially resolved cluster dynamics

Abstract Fission gas release within uranium dioxide nuclear fuel occurs as gas atoms diffuse through grains and arrive at grain boundary (GB) bubbles; these GB bubbles grow and interconnect with grain edge bubbles; and grain edge tunnels grow and connect to free surfaces. In this study, a hybrid multi-scale/multi-physics simulation approach is presented to investigate these mechanisms of fission gas release at the mesoscale. In this approach, fission gas production, diffusion, clustering to form intragranular bubbles, and re-solution within grains are included using spatially resolved cluster dynamics in the Xolotl code. GB migration and intergranular bubble growth and coalescence are included using the phase field method in the MARMOT code. This hybrid model couples Xolotl to MARMOT using the MultiApp and Transfer systems in the MOOSE framework, with Xolotl passing the arrival rate of gas atoms at GBs and intergranular bubble surfaces to MARMOT and MARMOT passing evolved GBs and bubble surface positions to Xolotl. The coupled approach performs well on the two-dimensional simulations performed in this work, producing similar results to the standard phase field model when Xolotl does not include fission gas clustering or re-solution. The hybrid model performs well computationally, with a negligible cost of coupling Xolotl and MARMOT and good parallel scalability. The hybrid model predicts that intragranular fission gas clustering and bubble formation results in up to 70% of the fission gas being trapped within grains, causing the increase in the intergranular bubble fraction to slow by a factor of six. Re-solution has a small impact on the fission gas behavior at 1800 K but it has a much larger impact at 1000 K, resulting in a twenty-times increase in the concentration of single gas atoms within grains. Due to the low diffusion rate, this increase in mobile gas atoms only results in a small acceleration in the growth of the intergranular bubble fraction. Finally, the hybrid model accounts for migrating GBs sweeping up gas atoms. This results in faster intergranular bubble growth with smaller initial grain sizes, since the additional GB migration results in more immobile gas clusters reaching GBs.

Kim, Dong-Uk↗

Stability of {sup 248–254}Cf isotopes against alpha and cluster radioactivity

Stability of {sup 248–254}Cf nuclei against alpha and cluster emissions is studied within our Coulomb and proximity potential model (CPPM). It is found that these nuclei are stable against light clusters (except alpha particles) and unstable against heavy cluster (A{sub 2}≥40) emissions. For heavy cluster emissions the daughter nuclei lead to doubly magic {sup 208}Pb or the neighbouring one. The effects of quadrupole and hexadecapole deformations of parent nuclei, daughter nuclei and emitted cluster on half lives are also studied. The computed alpha decay half life values (including quadrupole deformation β{sub 2}) are in close agreement with experimental data. Inclusion of quadrupole deformation reduces the height and width of the barrier (increases the barrier penetrability) and hence half life decreases. -- Highlights: •{sup 248–254}Cf parents are stable against light clusters (except alpha particles) and are unstable against heavy clusters ({sup 46}Ar, {sup 48,50}Ca etc.). •For the case of heavy cluster emissions the daughter nuclei are doubly magic {sup 208}Pb or neighbouring one. •The alpha decay half lives are in agreement with experimental data. •The cluster decay half lives decrease with the inclusion of quadrupole deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Grover-QAOA for 3-SAT: quadratic speedup, fair-sampling, and parameter clustering

Abstract The SAT problem is a prototypical NP-complete problem of fundamental importance in computational complexity theory with many applications in science and engineering; as such, it has long served as an essential benchmark for classical and quantum algorithms. This study shows numerical evidence for a quadratic speedup of the Grover Quantum Approximate Optimization Algorithm (G-QAOA) over random sampling for finding all solutions to 3-SAT (All-SAT) and Max-SAT problems. G-QAOA is less resource-intensive and more adaptable for these problems than Grover’s algorithm, and it surpasses conventional QAOA in its ability to sample all solutions. We show these benefits by classical simulations of many-round G-QAOA on thousands of random 3-SAT instances. We also observe G-QAOA advantages on the IonQ Aria quantum computer for small instances, finding that current hardware suffices to determine and sample all solutions. Interestingly, a single-angle-pair constraint that uses the same pair of angles at each G-QAOA round greatly reduces the classical computational overhead of optimizing the G-QAOA angles while preserving its quadratic speedup. We also find parameter clustering of the angles. The single-angle-pair protocol and parameter clustering significantly reduce obstacles to classical optimization of the G-QAOA angles.

Zhang, Zewen (ORCID:000000032258613X)↗

Laser Desorption Ionization Quadrupole Ion Trap Time-of-Flight Mass Spectrometry of Au{sub m}Fe{sub n}{sup +/–} Clusters Generated from Gold-Iron Nanoparticles and their Giant Nanoflowers. Electrochemical and/or Plasma Assisted Synthesis

Gold nanoparticles (NP) with average diameter ~100 nm synthesized from tetrachloroauric acid solution using stainless steel as a reducing agent were found to contain iron. Applying simultaneously high frequency (HF) plasma discharge in solution during the electrochemical reduction, giant gold-iron nanoflowers with average size ~1000–5000 nm were formed. Scanning electron microscopy (SEM) shows the morphology of the nanopowders produced as polygonal yet nearly spherical, whereas iron content in both products determined by energy dispersive X-ray analysis (EDX) was found to be at ~2.5 at. %. Laser desorption ionization (LDI) of both nanomaterials and mass spectrometric analysis show the formation of Au{sub m}Fe{sub n}{sup +/–} (m = 1–35; n = 1–3) clusters. Structure of few selected clusters in neutral or monocharged forms were computed by density functional theory (DFT) calculations and it was found that typical distances of an iron nucleus from adjacent gold nuclei lie in the interval 2.5 to 2.7 Å. Synthetized Au-Fe nanoparticles were found stable for at least 2 mo at room temperature (even in aqueous solution) without any stabilizing agent. Produced Au-Fe nanoparticles in combination with standard MALDI matrices enhance ionization of peptides and might find use in nanomedicine. .

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Stress-induced microcrack density evolution in β-eucryptite ceramics: Experimental observations and possible route to strain hardening

In order to investigate their microcracking behaviour, the microstructures of several β-eucryptite ceramics, obtained from glass precursor and cerammed to yield different grain sizes and microcrack densities, were characterized by laboratory and synchrotron x-ray refraction and tomography. Here, results were compared with those obtained from scanning electron microscopy (SEM). In SEM images, the characterized materials appeared fully dense but computed tomography showed the presence of pore clusters. Uniaxial tensile testing was performed on specimens while strain maps were recorded and analyzed by Digital Image Correlation (DIC). X-ray refraction techniques were applied on specimens before and after tensile testing to measure the amount of the internal specific surface (i.e., area per unit volume). X-ray refraction revealed that (a) the small grain size (SGS) material contained a large specific surface, originating from the grain boundaries and the interfaces of TiO 2 precipitates; (b) the medium (MGS) and large grain size (LGS) materials possessed higher amounts of specific surface compared to SGS material due to microcracks, which decreased after tensile loading; (c) the precursor glass had negligible internal surface. The unexpected decrease in the internal surface of MGS and LGS after tensile testing is explained by the presence of compressive regions in the DIC strain maps and further by theoretical arguments. It is suggested that while some microcracks merge via propagation, more close mechanically, thereby explaining the observed X-ray refraction results. Lastly, the mechanisms proposed would allow the development of a strain hardening route in ceramics.

36 MATERIALS SCIENCE↗

A machine learning approach to determine the elastic properties of printed fiber-reinforced polymers

This work focuses on the simultaneous determination of the elastic constants and the fiber orientation state for a short fiber-reinforced polymer composite by performing a minimum of experimental tests. Here we introduce a methodology that enables the inverse determination of fiber orientation state and the in-situ polymer properties by performing tensile tests at the composite coupon level. We demonstrate the approach for the extrusion deposition additive manufacturing (EDAM) process to illustrate one application of the methodology, but the development is such that it can be applied to short fiber-reinforced polymer (SFRP) systems processed via other methods. Currently, developing composites additive manufacturing digital twins require extensive material characterization. In particular, the mechanical characterization of the orthotropic elastic properties of a composite involves extensive sample preparation and testing, therefore the elasticity tensor is generally populated using a micromechanics model. This, however, requires measuring the fiber orientation state in addition to knowing the constituent material properties. Experimentally measuring the fiber orientation state can be tedious and time consuming. Further, optical methods are limited to resolving the orientation of cylindrical fibers or cluster of non-cylindrical fibers, and computed tomography (CT) methods scan regions of volume that are much smaller than a full printed bead. Therefore, we propose a methodology, accelerated by machine learning, to identify the anisotropic mechanical properties and fiber orientation state at the same time. Early results show that inference of the fiber orientation and composite properties is possible with as few as three tensile tests. Our results show that a combination of the choice of the micromechanics model and reliable set of experiments can yield the nine elastic constants, as well as, the fiber orientation state.

36 MATERIALS SCIENCE↗

The high level trigger and express data production at STAR

To meet the demands of the Beam Energy Scan phase-II (BES-II) program, the STAR experiment at the Relativistic Heavy Ion Collider (RHIC) developed a dual real-time framework consisting of a High Level Trigger (HLT) and an Express Data Production system (xProduction). The HLT operates online within the Data Acquisition (DAQ) chain on a dedicated multi-core CPU cluster with the option to offload compute-intensive kernels to Xeon Phi coprocessors. It uses parallelized algorithms, such as the Cellular Automaton (CA) Track Finder, to perform rapid tracking, vertexing, and event filtering. This allows it to select events of interest in real time and provide immediate feedback on detector and beam conditions. In contrast, the xProduction workflow runs concurrently and independently of the DAQ loop. It applies near offline-quality calibration and reconstruction within hours of data collection. The xProduction input is the express data stream, whose content can be enriched by HLT trigger/priority selections under DAQ/HLT resource constraints, and it uses the STAR calibration/conditions framework, incorporating online calibration/QA information when available. This enables early preliminary physics analysis, including the reconstruction of rare signals, such as hyperons and hypernuclei. It also provides collaboration-wide access to analysis-ready datasets. Together, the HLT and xProduction systems form a complementary architecture: the HLT performs online event selection while the xProduction chain delivers high-quality results within a short amount of time. This integrated framework has enabled the prompt reconstruction of the $^5_Λ$ He hypernucleus with high statistical significance and the efficient processing of hundreds of millions of heavy-ion collision events. In conclusion, its demonstrated scalability and robustness establish a model for future high-luminosity experiments requiring both online event filtering and rapid access to analysis-quality data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗

Full-Dimensional Quantum Dynamics of SiO in Collision with H 2

Here, we report the first full-dimensional potential energy surface (PES) and quantum mechanical close-coupling calculations for scattering of SiO due to H 2 . The full-dimensional interaction potential surface was computed using the explicitly correlated coupled-cluster (CCSD(T)-F12b) approach and fitted using an invariant polynomial approach. Pure rotational quenching cross sections from initial states v 1 = 0, j 1 = 1–5 of SiO in collision with H 2 are calculated for collision energies between 1.0 and 5000 cm –1 . State-to-state rotational rate coefficients are calculated at temperatures between 5 and 1000 K. The rotational rate coefficients of SiO with para-H 2 (p-H 2 ) are compared with previous approximate findings which were obtained using SiO-He PESs or scaled from SiO-He rate coefficients. Rovibrational state-to-state and total quenching cross sections and rate coefficients for initially excited SiO (v 1 = 1, j 1 = 0 and 1) in collisions with p-H 2 (v 2 = 0, j 2 = 0) and ortho-H 2 (o-H 2 ) (v 2 = 0, j 2 = 1) are also obtained. The application of the current collisional rate coefficients to astrophysics is briefly discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Preorganized Hydrogen Bond Network and Its Effect on Anion Stability

Rigid tricyclic locked in all axial 1,3,5-cyclohexanetriol derivatives with 0–3 trifluoromethyl groups were synthesized and photoelectron spectra of their conjugate bases and chloride anion clusters are reported along with density functional computations. The resulting vertical and adiabatic detachment energies provide measures of the anion stabilization due to the hydrogen bond network and inductive effects. The latter mechanism is found to be transmitted through space via hydrogen bonds

Environmental Molecular Sciences Laboratory↗

A Scalable Online Monitoring System Based on Elasticsearch for Distributed Data Acquisition in Cms

The part of the CMS Data Acquisition (DAQ) system responsible for data readout and event building is a complex network of interdependent distributed applications. To ensure successful data taking, these programs have to be constantly monitored in order to facilitate the timeliness of necessary corrections in case of any deviation from specified behaviour. A large number of diverse monitoring data samples are periodically collected from multiple sources across the network. Monitoring data are kept in memory for online operations and optionally stored on disk for post-mortem analysis. We present a generic, reusable solution based on an open source NoSQL database, Elasticsearch, which is fully compatible and non-intrusive with respect to the existing system. The motivation is to benefit from an off the-shelf software to facilitate the development, maintenance and support efforts. Elasticsearch provides failover and data redundancy capabilities as well as a programming language independent JSON-over-HTTP interface. The possibility of horizontal scaling matches the requirements of a DAQ monitoring system. The data load from all sources is balanced by redistribution over an Elasticsearch cluster that can be hosted on a computer cloud. In order to achieve the necessary robustness and to validate the scalability of the approach the above monitoring solution currently runs in parallel with an existing in-house developed DAQ monitoring system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Rationalization of the Hubbard U parameter in CeO{sub x} from first principles: Unveiling the role of local structure in screening

The density functional theory (DFT)+U method has been widely employed in theoretical studies on various ceria systems to correct the delocalization bias in local and semi-local DFT functionals with moderate computational cost. We present a systematic and quantitative study, aiming to gain better understanding of the dependence of Hubbard U on the local atomic arrangement. To rationalize the Hubbard U of Ce 4f, we employed the first principles linear response method to compute Hubbard U for Ce in ceria clusters, bulks, and surfaces. We found that the Hubbard U varies in a wide range from 4.3 eV to 6.7 eV, and exhibits a strong correlation with the Ce coordination number and Ce–O bond lengths, rather than the Ce 4f valence state. The variation of the Hubbard U can be explained by the changes in the strength of local screening due to O → Ce intersite transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CC2 oscillator strengths within the local framework for calculating excitation energies (LoFEx)

In a recent work [P. Baudin and K. Kristensen, J. Chem. Phys. 144, 224106 (2016)], we introduced a local framework for calculating excitation energies (LoFEx), based on second-order approximated coupled cluster (CC2) linear-response theory. LoFEx is a black-box method in which a reduced excitation orbital space (XOS) is optimized to provide coupled cluster (CC) excitation energies at a reduced computational cost. In this article, we present an extension of the LoFEx algorithm to the calculation of CC2 oscillator strengths. Two different strategies are suggested, in which the size of the XOS is determined based on the excitation energy or the oscillator strength of the targeted transitions. The two strategies are applied to a set of medium-sized organic molecules in order to assess both the accuracy and the computational cost of the methods. The results show that CC2 excitation energies and oscillator strengths can be calculated at a reduced computational cost, provided that the targeted transitions are local compared to the size of the molecule. Finally, to illustrate the potential of LoFEx for large molecules, both strategies have been successfully applied to the lowest transition of the bivalirudin molecule (4255 basis functions) and compared with time-dependent density functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗

Solving the nuclear pairing model with neural network quantum states

In this work, we present a variational Monte Carlo method that solves the nuclear many-body problem in the occupation number formalism exploiting an artificial neural network representation of the groundstate wave function. A memory-efficient version of the stochastic reconfiguration algorithm is developed to train the network by minimizing the expectation value of the Hamiltonian. We benchmark this approach against widely used nuclear many-body methods by solving a model used to describe pairing in nuclei for different types of interaction and different values of the interaction strength. Despite its polynomial computational cost, our method outperforms coupled-cluster and provides energies that are in excellent agreement with the numerically-exact full configuration interaction values.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗