This work explores the influence of blend composition, network architecture, and hydrogen bonding on the material properties of crosslinked epoxy networks, focusing on the glass transition temperature (Tg) and Young’s modulus (Y). We used coarse-grained molecular dynamics simulations to simulate varying compositions of stiff and flexible components in epoxy monomer blends with varying excess of curative. We find that, without hydrogen bonding, networks of any composition show a monotonically increasing Tg with decreasing excess curative, consistent with theory. In contrast, we find that when hydrogen bonding is introduced, the binary blend networks show significant enhancement in Tg for lightly crosslinked systems. This result contributes to an explanation of the anomalous Tg behavior observed experimentally in these systems. We further find that Y is generally enhanced by hydrogen bonds, especially below Tg, demonstrating that hydrogen bonding has a significant influence on mechanical properties and can allow access to other desirable dynamic behavior, especially self-healing.
Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.
We report significantly enhanced device performance in long wavelength interband cascade lasers (ICLs) by employing a recently proposed innovative quantum well (QW) active region containing strained InAsP layers. These ICLs were able to operate at wavelengths near 14.4 μm, the longest ever demonstrated for III-V interband lasers, implying great potential of ICLs to cover an even wider wavelength range. Also, by applying the aforesaid QW active region configuration on ICLs at relatively short wavelengths, ICLs were demonstrated at a low threshold current density (e.g., 13 A/cm2 at 80 K) and at temperatures up to 212 K near 12.4 μm, more than 50 K higher than the previously reported ICLs with the standard W-shape QW active region at similar wavelengths. This suggests that the QW active region with InAsP layers can be used to improve device performance at the shorter wavelengths.
Researchers have the potential to be exposed to a wide variety of hazards inherent to the equipment they use and maintain. When equipment does not function as expected, researchers sometimes reach out to their vendors for assistance. Early diagnostic or troubleshooting interactions between researcher and vendor are often conducted over the telephone and can lead to researchers performing work outside of their area of expertise and exposure to unknown hazards. This type of interaction significantly contributed to an incident where during diagnostic activities a researcher accidentally contacted, and discharged, a capacitor in an X-ray diffraction instrument. While this incident did not produce a serious injury, if the capacitor discharge path had occurred hand-to-hand across the heart, a serious injury may have been possible.
Electronic structure calculations on small systems such as H2, H2O, LiH, and BeH2 with chemical accuracy are still a challenge for the current generation of noisy intermediate-scale quantum (NISQ) devices. One of the reasons is that due to the device limitations, only minimal basis sets are commonly applied in quantum chemical calculations, which allows one to keep the number of qubits employed in the calculations at a minimum. However, the use of minimal basis sets leads to very large errors in the computed molecular energies as well as potential energy surface shapes. One way to increase the accuracy of electronic structure calculations is through the development of small basis sets better suited for quantum computing. In this work, we show that the use of adaptive basis sets, in which exponents and contraction coefficients depend on molecular structure, provides an easy way to dramatically improve the accuracy of quantum chemical calculations without the need to increase the basis set size and thus the number of qubits utilized in quantum circuits. As a proof of principle, we optimize an adaptive minimal basis set for quantum computing calculations on an H2 molecule, in which exponents and contraction coefficients depend on the H-H distance, and apply it to the generation of H2 potential energy surface on IBM-Q quantum devices. The adaptive minimal basis set reaches the accuracy of the double-zeta basis sets, thus allowing one to perform double-zeta quality calculations on quantum devices without the need to utilize twice as many qubits in simulations. This approach can be extended to other molecular systems and larger basis sets in a straightforward manner.
Wang, Hanyu; Self, Ethan C.; Addamane, Sadhvikas J.; Rouleau, Christopher M.; Wixom, Ryan R.; Browning, Katie L.; Veith, Gabriel M.; Liang, Liyuan; Browning, James F.
Here, we report deposition of hematite/Pd thin films on silicon wafers via radio frequency (RF) magnetron sputtering and subsequent characterization for future in situ neutron reflectometry studies. Following deposition, the hematite/Pd thin films were characterized as prepared and after annealing in air for 2h at 400, 500, and 600 °C, respectively. Raman spectroscopy, grazing incidence x-ray diffraction, and neutron reflectometry (NR) were used to characterize the structure and chemical compositions of the thin films. The results indicate that pure α-Fe2O3 (hematite) films were produced, free from other iron oxide phases and impurities. NR data reveal that one intermediate layer between the Pd layer and the hematite layer was formed during sputtering deposition processes. The fitted scattering length density (SLD) of the as-deposited hematite layer is 70% of the theoretical SLD value, indicating that the grains are loosely packed in the RF-deposited hematite films. After annealing at elevated temperatures, the hematite films show increased SLD values but remain comparable to that of preannealed.
Earlier studies have proven how ducted fuel injection (DFI) substantially reduces soot for low- and mid-load conditions in heavy-duty engines, without significant adverse effects on other emissions. Nevertheless, no comprehensive DFI study exists showing soot reductions at high- and full-load conditions. This study investigated DFI in a single-cylinder, 1.7-L, optical engine from low- to full-load conditions with a low-net-carbon fuel consisting of 80% renewable diesel and 20% biodiesel. Over the tested load range, DFI reduced engine-out soot by 38.1-63.1% compared to conventional diesel combustion (CDC). This soot reduction occurred without significant detrimental effects on other emission types. Thus, DFI reduced the severity of the soot-NOx tradeoff at all tested conditions. While DFI delivered considerable soot reductions in the present study, previous DFI studies at low- and mid-load conditions delivered larger soot reductions (>90%) compared to CDC operation at the same conditions. Therefore, the DFI configuration used here has been deemed nonoptimal (in terms of parameters such as the injector-spray and piston geometries), and several improvements are recommended for future studies with high-load DFI. These improvements include employing better spray-duct alignment, a deeper piston bowl with a smaller injector umbrella angle, and a fuel injector that opens and closes faster. The study also suggests future research to make DFI ready for commercialization, such as metal-engine tests to ensure desirable DFI performance over an engine's complete speed/load map. Overall, this study supports the continued development and commercialization of DFI to meet upcoming emissions regulations for heavy-duty vehicles. Specifically, multicylinder engine experiments and CFD simulations should be utilized to optimize the performance and clarify the full potential of DFI.
Physical fatigue can have adverse effects on humans in extreme environments. Therefore, being able to predict fatigue using easy to measure metrics such as heart rate (HR) signatures has potential to have an impact in real-life scenarios. We apply a functional logistic regression model that uses HR signatures to predict physical fatigue, where physical fatigue is defined in a data-driven manner. Data were collected using commercially available wearable devices on 47 participants hiking the 20.7-mile Grand Canyon rim-to-rim trail in a single day. Fitted model provides good predictions and interpretable parameters for real-life application.
A predictive thermodynamic model is utilized for the calculation of fuel properties of oxymethylene dimethyl ethers (OME3–4), surrogates for gasoline, diesel and aviation fuel, as well as alcohol blends with gasoline and diesel. The alcohols used for these blends are methanol, ethanol, propanol, butanol and pentanol; their mixing ratio ranges from 10 to 50% by volume. The model is based on the Perturbed-Chain Statistical Association Fluid Theory (PC-SAFT) equation of state (EoS) and Vapor Liquid Equilibrium (VLE) calculations at constant temperature, density and composition. The model includes the association term, with the assumption of two association sites (2B scheme), to enable the modeling of alcohols. The pure-component parameters are estimated based on the Group Contribution (GC) method of various sources, as well as a parametrization model specifically designed for the case of OME3–4. The results of the computational model for the density, vapor pressure and distillation curves at various conditions, including high-pressure, high-temperature (HPHT), are compared to experimental and computational data available in the literature. In the cases where no measurements are available for the surrogates, experimental data for the corresponding target fuel are used, taking into consideration the inherent deviation in properties between real and surrogate fuel. Overall, the results are in good agreement with the data from the literature, with the average deviation not exceeding 12% for temperature (Kelvin) on the distillation curves, 10% for density and 46% for vapor pressure and the general trend being captured successfully. The use of different pure component parameter estimation techniques can further improve the prediction quality in the cases of OME3–4 and the aviation fuel surrogate, especially for the vapor pressure, leading to an average deviation lower than 18%. These results demonstrate the predictive capabilities of the model, which extend to a wide range of fuel types and pressure/temperature conditions. Through this investigation, the present work aims to establish the limits of applicability of this thermodynamic property prediction methodology.
Porous liquids (PLs) based on the zeolitic imidazole framework ZIF-8 are attractive systems for carbon capture since the hydrophobic ZIF framework can be solvated in aqueous solvent systems without porous host degradation. However, solid ZIF-8 is known to degrade when exposed to CO2 in wet environments, and therefore the long-term stability of ZIF-8-based PLs is unknown. Through aging experiments, the long-term stability of a ZIF-8 PL formed using the water, ethylene glycol, and 2-methylimidazole solvent system was systematically examined, and the mechanisms of degradation were elucidated. The PL was found to be stable for several weeks, with no ZIF framework degradation observed after aging in N2 or air. However, for PLs aged in a CO2 atmosphere, formation of a secondary phase occurred within 1 day from the degradation of the ZIF-8 framework. From the computational and structural evaluation of the effects of CO2 on the PL solvent mixture, it was identified that the basic environment of the PL caused ethylene glycol to react with CO2 forming carbonate species. These carbonate species further react within the PL to degrade ZIF-8. The mechanisms governing this process involves a multistep pathway for PL degradation and lays out a long-term evaluation strategy of PLs for carbon capture. Additionally, it clearly demonstrates the need to examine the reactivity and aging properties of all components in these complex PL systems in order to fully assess their stabilities and lifetimes.
Here, using a quantum computer to speed up one step in a textbook approach to generating random numbers proves to be a savvy strategy, and one that could make good use of quantum computers that will be available in the near future.
Astra, deployed in 2018, was the first petascale supercomputer to utilize processors based on the ARM instruction set. The system was also the first under Sandia's Vanguard program which seeks to provide an evaluation vehicle for novel technologies that with refinement could be utilized in demanding, large-scale HPC environments. In addition to ARM, several other important first-of-a-kind developments were used in the machine, including new approaches to cooling the datacenter and machine. This article documents our experiences building a power measurement and control infrastructure for Astra. While this is often beyond the control of users today, the accurate measurement, cataloging, and evaluation of power, as our experiences show, is critical to the successful deployment of a large-scale platform. While such systems exist in part for other architectures, Astra required new development to support the novel Marvell ThunderX2 processor used in compute nodes. In addition to documenting the measurement of power during system bring up and for subsequent on-going routine use, we present results associated with controlling the power usage of the processor, an area which is becoming of progressively greater interest as data centers and supercomputing sites look to improve compute/energy efficiency and find additional sources for full system optimization.
Here, the effects of gamma and proton irradiation, and of forward bias minority carrier injection, on photo-response were investigated for InAsSb/AlAsSb pBn mid-wave infrared (MWIR) detectors with an engineered majority-carrier barrier. Room-temperature gamma irradiation had an insignificant effect on 77 K photo-response. Gamma irradiation at 77 K detector temperature, however, decreased in situ photo-response by 19% after a cumulative dose of ~ 500 krad(Si). Subsequent forward bias minority carrier injection had no effect on photo-response. The 77 K detectors irradiated with 30 MeV protons up to 2 Mrad(Si) had photo-response degraded by up to 70%, but here forward bias minority carrier (hole) injection caused up to 12% recovery that persisted more than 30 min. These results suggest a mitigation strategy for maintaining the photo-response of similar detectors in radiation environments that cause displacement damage defects.
Vanadium dioxide (VO2) undergoes a metal-insulator phase transition at ∼70 °C and has attracted substantial interest for potential applications in electronics, including those in neuromorphic computing. The vanadium-oxygen system has a rather complicated phase diagram, and controlling the stoichiometry and the phase of thin films of vanadium oxides is a well-known challenge. We explore the novel combination of two methods of VO2 thin film deposition using off-axis RF magnetron sputtering on (100)- and (111)-oriented yttria-stabilized zirconia (YSZ) substrates: reactive sputtering of vanadium in an oxygen environment and sputtering of vanadium metal followed by oxidation to VO2. Interestingly, the reactive sputtering process on both substrate orientations yields the metastable semiconducting VO2 (B) phase, which is structurally stabilized by the YSZ surface. The metal sputtering and oxidation process on YSZ produces mainly the equilibrium monoclinic (or M1) phase of VO2 that exhibits a metal-insulator transition. Using this method, we obtained thin films of (010)-textured polycrystalline VO2 (M1) that show a metal-insulator transition with an on/off ratio larger than 1000.
The ability to rapidly screen material performance in the vast space of high entropy alloys is of critical importance to efficiently identify optimal hydride candidates for various use cases. Given the prohibitive complexity of first principles simulations and large-scale sampling required to rigorously predict hydrogen equilibrium in these systems, we turn to compositional machine learning models as the most feasible approach to screen on the order of tens of thousands of candidate equimolar high entropy alloys (HEAs). Critically, we show that machine learning models can predict hydride thermodynamics and capacities with reasonable accuracy (e.g. a mean absolute error in desorption enthalpy prediction of ∼5 kJ molH2−1) and that explainability analyses capture the competing trade-offs that arise from feature interdependence. We can therefore elucidate the multi-dimensional Pareto optimal set of materials, i.e., where two or more competing objective properties can't be simultaneously improved by another material. This provides rapid and efficient down-selection of the highest priority candidates for more time-consuming density functional theory investigations and experimental validation. Various targets were selected from the predicted Pareto front (with saturation capacities approaching two hydrogen per metal and desorption enthalpy less than 60 kJ molH2−1) and were experimentally synthesized, characterized, and tested amongst an international collaboration group to validate the proposed novel hydrides. Additional top-predicted candidates are suggested to the community for future synthesis efforts, and we conclude with an outlook on improving the current approach for the next generation of computational HEA hydride discovery efforts.
Phase formation and stability of five component compositionally complex rare earth zirconates (5RE 2 Zr 2 O 7 ) were investigated by X‐ray diffraction and electron microprobe analysis. Zirconates with different rare earth compositions (LaNdSmEuDy, LaNdSmEuYb, LaNdEuErYb, LaNdDyErYb, SmEuDyYHo, LaYHoErYb, and DyYHoErYb) were synthesized at 1700°C and 2000°C by the solid‐state method to investigate the effect of A‐site site disorder (δ A ) on phase stability. Increased site disorder results from mixed cation occupancy with localized crystallographic strain and bond disorder. Compositions LaNdSmEuDy (δ A = 4.6) and LaNdSmEuYb (δ A = 6.0) produced a single pyrochlore phase and compositions SmDyYHoErYb (δ A = 2.8), LaYHoErYb (δ A = 6.2), and DyYHoErYb (δ A = 1.7) produced a single fluorite phase. High δ A compositions LaNdEuErYb (δ A = 6.9) and LaNdDyErYb (δ A = 7.2) produced a pyrochlore and fluorite phase mixture at 1700°C. Single phase was obtained for the latter composition at 2000°C. Of the single phase compositions calcined at 1700°C, LaNdSmEuYb and LaYHoErYb (both with largest δ A ) showed decomposition to mixed fluorite and pyrochlore phases during lower temperature anneals, indicating entropic stabilization. Comparison with prior work shows a temperature dependence of the critical δ A for phase stability, and compositions near it are expected to be entropy stabilized.
The penetration of wind power generation into the power grid has been accelerated in recent times due to the aggressive emission targets set by governments and other regulatory authorities. Although wind power has the advantage of being environment-friendly, wind as a resource is intermittent in nature. In addition, wind power contributes little inertia to the system as most wind turbines are connected to the grid via power electronic converters. These negative aspects of wind power pose serious challenges to the frequency security of power systems as penetration increases. In this work, an approach is proposed where an energy storage system (ESS) is used to mitigate frequency security issues of wind-integrated systems. ESSs are well equipped to supply virtual inertia to the grid due to their fast-acting nature, thus replenishing some of the energy storage capability of displaced inertial generation. In this work, a probabilistic approach is proposed to estimate the amount of inertia required by a system to ensure frequency security. Reduction in total system inertia due to the displacement of conventional synchronous generation by wind power generation is considered in this approach, while also taking into account the loss of inertia due to forced outages of conventional units. Monte Carlo simulation is employed for implementing the probabilistic estimation of system inertia. An ESS is then sized appropriately, using the system swing equation, to compensate for the lost inertia. The uncertainty associated with wind energy is modeled into the framework using an autoregressive moving average technique. Effects of increasing the system peak load and changing the wind profile on the expected system inertia are studied to illustrate various factors that might affect system frequency security. The proposed method is validated using the IEEE 39-bus test system.
Jiang, Kunyao; Tang, Jingyu; House, Stephen D.; Xu, Chengchao; Xiao, Kelly; Porter, Lisa M.; Davis, Robert F.
Ga2O3 films were deposited on (100) MgAl2O4 spinel substrates at 550, 650, 750, and 850 °C using metal-organic chemical vapor deposition and investigated using x-ray diffraction and transmission electron microscopy. A phase-pure γ-Ga2O3-based material having an inverse spinel structure was formed at 850 °C; a mixture of the γ-phase and β-Ga2O3 was detected in films grown at 750 °C. Only β-Ga2O3 was determined in the films deposited at 650 and 550 °C. A β- to γ-phase transition occurred from the substrate/film interface during growth at 750 °C. The growth and stabilization of the γ-phase at the outset of film growth at 850 °C was affected by the substantial Mg and Al chemical interdiffusion from the MgAl2O4 substrate observed in the energy-dispersive x-ray spectrum. Further, atomic-scale investigations via scanning transmission electron microscopy of the films grown at 750 and 850 °C revealed a strong tetrahedral site preference for Ga and an octahedral site preference for Mg and Al. It is postulated that the occupation of these atoms in these particular sites drives the β-Ga2O3 to γ-phase transition and markedly enhances the thermal stability of the latter phase at elevated temperatures.
You, Huaiqian Q.; Xu, Xiao; Yu, Yue; D'Elia, Marta; Foster, John; Silling, Stewart A.
Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.
Modern density functional theory (DFT) is a powerful tool for accurately predicting self-consistent material properties such as equations of state, transport coefficients and opacities in high energy density plasmas, but it is generally restricted to conditions of local thermodynamic equilibrium (LTE) and produces only averaged electronic states instead of detailed configurations. We propose a simple modification to the bound-state occupation factor of a DFT-based average-atom model that captures essential non-LTE effects in plasmas—including autoionization and dielectronic recombination—thus extending DFT-based models to new regimes. Finally, we then expand the self-consistent electronic orbitals of the non-LTE DFT-AA model to generate multi-configuration electronic structure and detailed opacity spectra. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.
High-altitude balloons carrying infrasound sensor payloads can be leveraged toward monitoring efforts to provide some advantages over other sensing modalities. On 10 July 2020, three sets of controlled surface explosions generated infrasound waves detected by a high-altitude floating sensor. One of the signal arrivals, detected when the balloon was in the acoustic shadow zone, could not be predicted via propagation modeling using a model atmosphere. Considering that the balloon’s horizontal motion showed direct evidence of gravity waves, we examined their role in infrasound propagation. Implementation of gravity wave perturbations to the wind field explained the signal detection and aided in correctly predicting infrasound travel times. Our results show that the impact of gravity waves is negligible below 20 km altitude; however, their effect is important above that height. The results presented here demonstrate the utility of balloon-borne acoustic sensing toward constraining the source region of variability, as well as the relevance of complexities surrounding infrasound wave propagation at short ranges for elevated sensing platforms.
Understanding titanium particle combustion processes is critical not only for characterizing existing pyrotechnic systems but also for creating new igniter designs. In order to characterize titanium particle combustion processes, morphologies, and temperatures, simultaneous spatially-resolved electric field holography and imaging pyrometry techniques were used to capture post-ignition data at up to 7 kHz. Due to the phase and thermal distortions present in the combustion cloud, traditional digital in-line holography techniques fail to capture accurate data. In this work, electric field holography techniques are used in order to cancel distortions and capture the three-dimensional spatial locations and diameters of the particles. In order to estimate the projected surface temperatures of the titanium particles, an imaging pyrometry method that ratios emission at 750 and 850 nm is utilized. Using these diagnostics, joint statistics are collected for particle size, morphology, velocity, and temperature. Results show that, early in the combustion process, the titanium particles are primarily oxidized by potassium perchlorate inside the igniter cup, resulting in projected surface temperatures near 3000 K. Later in the process, the particles interact with ambient air, resulting in lower surface temperatures around 2400 K and the formation of flame zones. These results are consistent with adiabatic flame temperature predictions as well as particle morphology observations of a titanium core with a TiO2 surface. Late stage particle expansion, star fragmentation, and molten droplet breakup events are also observed using the time-resolved morphology and temperature diagnostics. These results illustrate the different stages of titanium particle combustion in pyrotechnic environments, which can be used to inform improvements in next-generation igniters.
It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. We also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.
We develop an adaptive method for quantum state preparation that utilizes randomness as an essential component and that does not require classical optimization. Instead, a cost function is minimized to prepare a desired quantum state through an adaptively constructed quantum circuit, where each adaptive step is informed by feedback from gradient measurements in which the associated tangent space directions are randomized. We provide theoretical arguments and numerical evidence that convergence to the target state can be achieved for almost all initial states. We investigate different randomization procedures and develop lower bounds on the expected cost function change, which allows for drawing connections to barren plateaus and for assessing the applicability of the algorithm to large-scale problems.
SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.
Avoiding stress concentrations is essential to achieve robust parts since failure tends to originate at such concentrations. With recent advances in multimaterial additive manufacturing, it is possible to alter the stress (or strain) distribution by adjusting the material properties in selected locations. Here, we investigate the use of grayscale digital light processing (g-DLP) 3D printing to create modulus gradients around areas of high stress. These gradients prevent failure by redistributing high stresses (or strains) to the neighboring material. The improved material distributions are calculated using finite element analysis. The much-enhanced properties are demonstrated experimentally for thin plates with circular, triangular, and elliptical holes. This work suggests that multimaterial additive manufacturing techniques like g-DLP printing provide a unique opportunity to create tougher engineering materials and parts.
Many highly pixelated organic scintillator detection systems would benefit from independent readout of each scintillator pixel. Recent advances in Silicon Photomultiplier (SiPM) technology makes this goal feasible, however the data acquisition from potentially hundreds or thousands of channels requires a low-cost and compact solution. For pixelated neutron detection with organic scintillators, the capability to distinguish between neutron and gamma interactions using Pulse Shape Discrimination (PSD) is required along with pulse charge and time of arrival. The TOFPET2 ASIC from PETsys Electronics is a 64-channel readout chip providing pulse time and charge integration measurements from SiPMs, and is specifically designed for time-of-flight positron-emission tomography. Using an 8 × 8 array of 6 mm × 6 mm J-series SiPMs from SensL/OnSemi (ArrayJ-60035-64P-PCB), we have studied the energy and PSD performance of the TOFPET2 ASIC using a 4 × 4 array of 6 mm × 6 mm × 30 mm trans-Stilbene crystals from Inrad Optics and a custom SiPM routing board from PETsys Electronics. Using a time-over-threshold method, we measure a maximum PSD figure-of-merit of approximately 1.2 at 478 keV (the Compton edge of 662 keV) for a J-series SiPM operating at an over-voltage of 3V.
Plastic scintillators are widely used as radiation detection media in homeland security and nuclear physics applications. Their attributes include low cost, scalability to large detector volumes, and additive compounding to enable additional material and detection features, such as pulse shape discrimination (PSD), gamma-ray spectroscopy, aging resistance, and coincidence timing. However, traditional chemically cured plastic scintillators (CCS) require long reaction times, and hazardous wet chemical procedures performed by specially trained personnel, and can leave residual monomer, resulting in deleterious optical and material properties. Here, we synthesize melt blended scintillators (MBSs) in 2.5 days using easily accessible solid-state compounding of commercially-available poly(styrene) with 30-60 wt% fluorene-based compound 'P2' to create monolithic detectors with < 100 ppm residual monomer, in several form factors. The best scintillation performance was recorded for 60 wt% P2 in Styron 665, including gamma-ray light yield 139% of EJ- 200 commercial scintillator and PSD figure of merit (FOM) value of 2.65 at 478 keVee, approaching P2 organic glass scintillator (OGS). The capability of MBS to generate fog-resistant scintillators and poly(methyl methacrylate) (PMMA)-based scintillators for use in challenging environments is also demonstrated.