Results for “complex computing” — Page 49

DOE OSTI2022

Modeling deformation and failure in AlSi-polyester abradable sealcoating material using microstructure-based finite element simulation

A plasma-sprayed aluminum-silicon (AlSi)/polyester coating is applied in modern gas turbine engines as an abradable sealcoating to maintain tight clearances between the rotating blades and the static casing. While running the engine, the rotating blades “rub” with the abradable coating, which results in extreme strain rate (up to 10 6 ) dynamics and a high-temperature environment. Due to the difficulty of collecting direct measurements, predictive computational models are important for analyzing the deformation and failure of the abradable material, to help meet the design target of avoiding damage to the blade tip and maintaining high fuel efficiency. In this research, a microstructure-based finite element (FE) computational model was developed to capture the complex mechanical behavior of the AlSi/polyester microstructure. The model is based on a virtual representative-volume-element (RVE) of a metal-polymer microstructure, reconstructed from x-ray computed tomography. It models the plastic deformation of, and damage to, each AlSi and polyester constituents, as well as the failure at their interface. The model was calibrated and validated with uniaxial tension and compression experiments, conducted at two temperatures (298 K and 533 K) at an applied strain rate of 10 3 - 10 4 s -1 . The material exhibited strongly asymmetric tension-compression behavior and a sensitivity to temperature, which was well captured by the model. The model was further applied to investigate changes in mechanical behavior due to variations in constituents’ volume fractions, which provides guidance to the microstructural design of AlSi/polyester abradable materials. The model is expected to facilitate the development of improved abradable materials by bypassing the conventional trial-and-error approach and extensive testing requirements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

M + (M=Ca, Ba) Cations Bound to Molecular Cavities: A New Strategy for Incorporating Molecular Quantum States into Quantum Information

This project pursued a novel strategy for incorporating multiple qubits per ion into ion-trap based quantum computing (ITQC) involving Ca + and Ba + . By forming molecular complexes of these cations with molecular-scale cages, we hypothesized that molecular energy levels could be incorporated into quantum computing while retaining key properties of the atomic ions intact. We experimented with a variety of molecular cages and found that Na + , K + , Rb + , Ca 2+ , Sr 2+ , and Ba 2+ could be captured and brought into the gas phase efficiently by imbedding them inside [2.2.2]-benzocryptand. IR and UV spectra of these cage complexes are sensitive to the size and charge state of the ion, reporting on the structures and binding properties of the cage complexes. UV photofragmentation of the Ba 2+ -Acetate -1 -BzCrypt complex produces Ba + -BzCrypt, the complex targeted for exploration in the original hypothesis. Follow-on funding is needed to pursue the spectroscopy of this complex as a target for ITQC.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗

DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗

DOE OSTI2021

Mesoscopic Modeling and Rapid Simulation of Incremental Changes in Epidemic Scenarios on GPUs

In simulation-based studies and analyses of epidemics, a major challenge lies in resolving the conflict between fidelity of models and the speed of their simulation. Another related challenge arises in dealing with the large number of what–if scenarios that need to be explored. Here, we describe new computational methods that together provide an approach to dealing with both challenges. A mesoscopic modeling approach is described that strikes a middle ground between macroscopic models based on coupled differential equations and microscopic models built on fine-grained behaviors at the individual entity level. The mesoscopic approach offers the ability to incorporate complex compositions of multiple layers of dynamics even while retaining the potential for aggregate behaviors at varying levels. It also is an excellent match to the accelerator-based architectures of modern computing platforms in which graphical processing units (GPUs) can be exploited for fast simulation via the parallel execution mode of single instruction multiple thread (SIMT). The challenge of simulating a large number of scenarios is addressed via a method of sharing model state and computation across a tree of what–if scenarios that are localized, incremental changes to a large base simulation. A combination of the mesoscopic modeling approach and the incremental what–if scenario tree evaluation has been implemented in the software on modern GPUs. Synthetic simulation scenarios are presented to demonstrate the computational characteristics of our approach. Results from the experiments with large population data, including USA, UK, and India, illustrate the modeling methodology and computational performance on thousands of synthetically generated what–if scenarios. Execution of our implementation scaled to 8192 GPUs of supercomputing platforms demonstrates the ability to rapidly evaluate what–if scenarios several orders of magnitude faster than the conventional methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

HDBind: encoding of molecular structure with hyperdimensional binary representations

Traditional methods for identifying “hit” molecules from a large collection of potential drug-like candidates rely on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have a significant limitation in that they require exceptional computing capabilities for even relatively small collections of molecules. Increasingly large and complex state-of-the-art deep learning approaches have gained popularity with the promise to improve the productivity of drug design, notorious for its numerous failures. However, as deep learning models increase in their size and complexity, their acceleration at the hardware level becomes more challenging. Hyperdimensional Computing (HDC) has recently gained attention in the computer hardware community due to its algorithmic simplicity relative to deep learning approaches. The HDC learning paradigm, which represents data with high-dimension binary vectors, allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas (computer vision, bioinformatics, mass spectrometery, remote sensing, edge devices, etc.). To the best of our knowledge, our work is the first to consider HDC for the task of fast and efficient screening of modern drug-like compound libraries. We also propose the first HDC graph-based encoding methods for molecular data, demonstrating consistent and substantial improvement over previous work. We compare our approaches to alternative approaches on the well-studied MoleculeNet dataset and the recently proposed LIT-PCBA dataset derived from high quality PubChem assays. We demonstrate our methods on multiple target hardware platforms, including Graphics Processing Units (GPUs) and Field Programmable Gate Arrays (FPGAs), showing at least an order of magnitude improvement in energy efficiency versus even our smallest neural network baseline model with a single hidden layer. Our work thus motivates further investigation into molecular representation learning to develop ultra-efficient pre-screening tools. We make our code publicly available at https://github.com/LLNL/hdbind.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Adiabatic quantum support vector machines

Adiabatic quantum computers can solve difficult optimization problems (e.g., the quadratic unconstrained binary optimization problem), and they seem well suited to train machine learning models. In this paper, we describe an adiabatic quantum approach for training support vector machines. We show that the time complexity of our quantum approach is an order of magnitude better than the classical approach. Next, we compare the test accuracy of our quantum approach against a classical approach that uses the Scikit-learn library in Python across five benchmark datasets (Iris, Wisconsin Breast Cancer (WBC), Wine, Digits, and Lambeq). We show that our quantum approach obtains accuracies on par with the classical approach. Finally, we perform a scalability study in which we compute the total training times of the quantum approach and the classical approach with an increasing number of features and an increasing number of data points in the training dataset. In conclusion, our scalability results show that the quantum approach obtains a 3.5–4.5x speedup over the classical approach on datasets with many (millions of) features.

Computational Complexity↗

DOE OSTI2023

Simultaneous estimation of multiple eigenvalues with short-depth quantum circuit on early fault-tolerant quantum computers

We introduce a multi-modal, multi-level quantum complex exponential least squares (MM-QCELS) method to simultaneously estimate multiple eigenvalues of a quantum Hamiltonian on early fault-tolerant quantum computers. Our theoretical analysis demonstrates that the algorithm exhibits Heisenberg-limited scaling in terms of circuit depth and total cost. Notably, the proposed quantum circuit utilizes just one ancilla qubit, and with appropriate initial state conditions, it achieves significantly shorter circuit depths compared to circuits based on quantum phase estimation (QPE). Numerical results suggest that compared to QPE, the circuit depth can be reduced by around two orders of magnitude under several settings for estimating ground-state and excited-state energies of certain quantum systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Electronic structure with direct diagonalization on a D-wave quantum annealer

Quantum chemistry is regarded to be one of the first disciplines that will be revolutionized by quantum computing. Although universal quantum computers of practical scale may be years away, various approaches are currently being pursued to solve quantum chemistry problems on near-term gate-based quantum computers and quantum annealers by developing the appropriate algorithm and software base. This work implements the general Quantum Annealer Eigensolver (QAE) algorithm to solve the molecular electronic Hamiltonian eigenvalue-eigenvector problem on a D-Wave 2000Q quantum annealer. The approach is based on the matrix formulation, efficiently uses qubit resources based on a power-of-two encoding scheme and is hardware-dominant relying on only one classically optimized parameter. We demonstrate the use of D-Wave hardware for obtaining ground and excited electronic states across a variety of small molecular systems. The approach can be adapted for use by a vast majority of electronic structure methods currently implemented in conventional quantum-chemical packages. The results of this work will encourage further development of software such as qbsolv which has promising applications in emerging quantum information processing hardware and has expectation to address large and complex optimization problems intractable for classical computers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Fast All-sky Radiation Model for Solar applications (FARMS) [SWR-16-18]

The Fast All-sky Radiation Model for Solar applications (FARMS) is used to compute cloudy irradiance. Radiative transfer (RT) models simulating broadband solar radiation have been widely used by atmospheric scientists to model solar resources for various energy applications such as operational forecasting. Due to the complexity of solving the RT equation, the computation under cloudy conditions can be extremely time consuming though many approximations (e.g. two-stream approach and delta-M truncation scheme) have been utilized. Thus, a more efficient RT model is crucial for model developers as a new option for approximating solar radiation at the land surface with minimal loss of accuracy. We have developed a fast all-sky radiation model for solar applications (FARMS) using the simplified clear-sky RT model, REST2, and simulated cloud transmittances and reflectances from the Rapid Radiation Transfer Model (RRTM) with a sixteen-stream Discrete Ordinates Radiative Transfer (DISORT). Simulated lookup tables (LUTs) of cloud transmittances and reflectances were created by varying cloud optical thicknesses, cloud particle sizes, and solar zenith angles. Equations with optimized parameters were fitted to the cloud transmittances and reflectances to develop the model. Using this model the all-sky solar irradiance at the land surface can be computed rapidly by combining REST2 with the cloud transmittances and reflectances. This new RT model is more than 1000 times faster than those currently utilized in solar resource assessment and forecasting since it does not explicitly solve the RT equation for each individual cloud condition. Our results indicate the accuracy of the fast radiative transfer model is comparable to or better than two-stream approximation in term of computing cloud transmittance and solar radiation.

Xie, Yu↗

DOE OSTI2023

Fast All-sky Radiation Model for Solar applications (FARMS) 2023 [SWR-16-18]

The Fast All-sky Radiation Model for Solar applications (FARMS) is used to compute cloudy irradiance. Radiative transfer (RT) models simulating broadband solar radiation have been widely used by atmospheric scientists to model solar resources for various energy applications such as operational forecasting. Due to the complexity of solving the RT equation, the computation under cloudy conditions can be extremely time consuming though many approximations (e.g. two-stream approach and delta-M truncation scheme) have been utilized. Thus, a more efficient RT model is crucial for model developers as a new option for approximating solar radiation at the land surface with minimal loss of accuracy. We have developed a fast all-sky radiation model for solar applications (FARMS) using the simplified clear-sky RT model, REST2, and simulated cloud transmittances and reflectances from the Rapid Radiation Transfer Model (RRTM) with a sixteen-stream Discrete Ordinates Radiative Transfer (DISORT). Simulated lookup tables (LUTs) of cloud transmittances and reflectances were created by varying cloud optical thicknesses, cloud particle sizes, and solar zenith angles. Equations with optimized parameters were fitted to the cloud transmittances and reflectances to develop the model. Using this model the all-sky solar irradiance at the land surface can be computed rapidly by combining REST2 with the cloud transmittances and reflectances. This new RT model is more than 1000 times faster than those currently utilized in solar resource assessment and forecasting since it does not explicitly solve the RT equation for each individual cloud condition. Our results indicate the accuracy of the fast radiative transfer model is comparable to or better than two-stream approximation in term of computing cloud transmittance and solar radiation.

Xie, Yu↗

DOE OSTI2024

Modeling Clustered DNA Damage by Ionizing Radiation Using Multinomial Damage Probabilities and Energy Imparted Spectra

Simple and complex clustered DNA damage represent the critical initial damage caused by radiation. In this paper, a multinomial probability model of clustered damage is developed with probabilities dependent on the energy imparted to DNA and surrounding water molecules. The model consists of four probabilities: (A) direct damage of sugar-phosphate moieties leading to SSB, (B) OH− radical formation with subsequent SSB and BD formation, (C) direct damage to DNA bases, and (D) energy imparted to histone proteins and other molecules in a volume not leading to SSB or BD. These probabilities are augmented by introducing probabilities for the relative location of SSB using a ≤10 bp criteria for a double-strand break (DSB) and for the possible success of a radical attack that leads to SSB or BD. Model predictions for electrons, 4He, and 12C ions are compared to the experimental data and show good agreement. Thus, the developed model allows an accurate and rapid computational method to predict simple and complex clustered DNA damage as a function of radiation quality and to explore the resulting challenges to DNA repair.

Biochemistry & Molecular Biology↗

DOE OSTI2021

Immobilization of “Capping Arene” Cobalt(II) Complexes on Ordered Mesoporous Carbon for Electrocatalytic Water Oxidation

We report the synthesis, characterization, and electrocatalytic water oxidation activity of two cobalt complexes, (6-FP)Co(NO 3 ) 2 (1) (6-FP = 8,8'-(1,2-phenylene)diquinoline) and (5-FP)Co(NO 3 ) 2 (2) (5-FP = 1,2-bis(N-7-azaindolyl)benzene), containing "capping arene" bidentate ligands with nitrogen atom donors. The cobalt complexes 1 and 2 were supported on ordered mesoporous carbon (OMC) by π-π stacking, resulting in heterogenized cobalt materials 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, respectively, and studied for electrocatalytic water oxidation. We find that 6-FP-Co-OMC-1 exhibits an overpotential of 355 mV for a current density of 10 mA cm -2 and a turnover frequency (TOF) of ~0.53 s -1 at an overpotential of 400 mV at pH 14. 6-FP-Co-OMC-1 exhibits activity that is ~1.6 times that of 5-FP-Co-OMC-2, which gives a TOF of 0.32 s -1 at 400 mV overpotential. The structural stability of the single-atom Co site was demonstrated for 6-FP-Co-OMC-1 using X-ray absorption spectroscopy for the molecular complex supported on OMC, but slow degradation in catalyst activity can be attributed to eventual formation of Co oxide clusters. DFT computations of electrocatalytic water oxidation using the molecular complexes as models provide a description of the catalytic mechanism. These studies reveal that the mechanism for O-O bond formation involves an intermediate Co IV oxo complex that undergoes an intramolecular reductive O-O coupling to form a Co II -OOH species. Further, the calculations predict that the molecular 6-FP-Co structure is more active for electrocatalytic water oxidation than 5-FP-Co, which is consistent with experimental studies of 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, highlighting the possibility that the ligand structure influences the catalytic activity of the supported molecular catalysts.

25 ENERGY STORAGE↗

DOE OSTI2022

Virtual Infrastructure Twins: Software Testing Platforms for Computing-Instrument Ecosystems

Science ecosystems are being built by federating computing systems and instruments located at geographically distributed sites over wide-area networks. These computing-instrument ecosystems are expected to support complex workflows that incorporate remote, automated AI-driven science experiments. Their realization, however, requires various designs to be explored and software components to be developed, in order to support the orchestration of distributed computations and experiments. It is often too expensive, infeasible, or disruptive for the entire ecosystem to be available during the typically long software development and testing periods. We propose a Virtual Infrastructure Twin (VIT) of the ecosystem that emulates its network and computing components, and incorporates its instrument software simulators. It provides a software environment nearly identical to the ecosystem to support early development and testing, and design space exploration. We present a brief overview of previous digital infrastructure twins that culminated in the VIT concept, including (i) the virtual science network environment for developing software-defined networking solutions, and (ii) the virtual federated science instrument environment for testing the federation software stack and remote instrument control software. We briefly describe VITs for Nion microscope steering and access to GPU systems.

Rao, Nageswara↗

DOE OSTI2023

Virtual Infrastructure Twins: Software Testing Platforms for Computing-Instrument Ecosystems

Science ecosystems are being built by federating computing systems and instruments located at geographically distributed sites over wide-area networks. These computing-instrument ecosystems are expected to support complex workflows that incorporate remote, automated AI-driven science experiments. Their realization, however, requires various designs to be explored and software components to be developed, in order to support the orchestration of distributed computations and experiments. It is often too expensive, infeasible, or disruptive for the entire ecosystem to be available during the typically long software development and testing periods. We propose a Virtual Infrastructure Twin (VIT) of the ecosystem that emulates its network and computing components, and incorporates its instrument software simulators. It provides a software environment nearly identical to the ecosystem to support early development and testing, and design space exploration. We present a brief overview of previous digital infrastructure twins that culminated in the VIT concept, including (i) the virtual science network environment for developing software-defined networking solutions, and (ii) the virtual federated science instrument environment for testing the federation software stack and remote instrument control software. We briefly describe VITs for Nion microscope steering and access to GPU systems.

Rao, Nageswara↗

DOE OSTI2025

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗

DOE OSTI2021

CFD Evaluation of Pressure Change Along Coolant Passages in Sodium-Cooled Fast Reactor with Nek5000

To support the design efforts of advanced sodium-cooled fast reactors (SFRs), a series of computational fluid dynamics (CFD) simulations are performed to investigate the pressure change along various flow passages in the proposed SFR system. The simulations are carried out with the state-of-the-art spectral element flow solver, Nek5000. Two specific case studies are presented in this paper: the flow exiting the axial neutron reflector channels and the flow entering the fuel pin bundle. Due to the high Reynolds numbers expected, a Reynolds-averaged Navier-Stokes (RANS) approach is necessary to model the turbulence. A newly developed regularized k-omega RANS model is adopted in the related CFD calculations. The first case study explores the effect of Reynolds number on the pressure change when flow exits the reflector channels. The pressure change in this case has two major contributors: the change due to wall friction and the Bernoulli effect. It is noted that the nondimensional pressure loss follows a log-linear trend up to Re = 10 5 , and then the trend is flattened. In the second case study, the advanced NekNek coupling capability is tested where an integral domain can be divided into multiple subdomains with coupling interfaces, which would greatly ease the meshing process of complex engineering geometries and potentially save computational resources. The preliminary results obtained so far confirm the consistency between the NekNek results and those produced by regular Nek5000 simulation. The presented work demonstrates the readiness and flexibility of the related CFD techniques, which is part of the broader effort to leverage cutting-edge CFD to inform the advanced nuclear reactor designs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗