Publications

Ge-Sb-Te (GST) alloys are leading phase-change materials for data storage due to the fast phase transition between amorphous and crystalline states. Ongoing research aims at improving the stability of the amorphous phase to improve retention. This can be accomplished by the introduction of carbon as a dopant to Ge2Sb2Te5, which is known to alter the short- and mid-range structure of the amorphous phase and form covalently bonded C clusters, both of which hinder crystallization. The relative importance of these processes as a function of C concentration is not known. We used molecular dynamics simulation based on density functional theory to study how carbon doping affects the atomic structure of GST-C. Carbon doping results in an increase in tetrahedral coordination, especially of Ge atoms, and this is known to stabilize the amorphous phase. We observe an unexpected, non-monotonous trend in the number of tetrahedral bonded Ge with the amount of carbon doping. Our simulations show an increase in the number of tetrahedral bonded Ge up to 5 at.% C, after which the number saturates and begins to decrease above 14 at.% C. The carbon atoms aggregate into clusters, mostly in the form of chains and graphene flakes, leaving less carbon to disrupt the GST matrix at higher carbon concentrations. Different degrees of carbon clustering can explain divergent experimental results for recrystallization temperature for carbon doped GST.

High-nickel-content layered oxides are among the most promising electric vehicle battery cathode materials. However, their interfacial reactivity with electrolytes and tendency toward oxygen release (possibly yielding reactive 1O2) remain degradation concerns. Elucidating the most relevant (i.e., fastest) interfacial degradation mechanism will facilitate future mitigation strategies. We apply screened hybrid density functional (HSE06) calculations to compare the reaction kinetics of LixNiO2 surfaces with ethylene carbonate (EC) with those of O2 release. On both the (001) and (104) facets, EC oxidative decomposition exhibits lower activation energies than O2 release. Our calculations, coupled with previously computed liquid-phase reaction rates of 1O2 with EC, strongly question the role of “reactive 1O2” species in electrolyte oxidative degradation. The possible role of other oxygen species is discussed. To deal with the challenges of modeling LixNiO2 surface reactivity, we emphasize a “local structure” approach instead of pursuing the global energy minimum.

During February 2023, a total of 32 individual DAS systems acted jointly as a global seismic monitoring network. The aim of this Global DAS Month campaign was to coordinate a diverse network of organizations, instruments, and file formats in order to gain knowledge and move toward the next generation of earthquake monitoring networks. During this campaign, 156 earthquakes of magnitude 5 or larger were reported by the USGS and contributors shared data for 60 min after each event’s origin time. Participating systems represent a variety of manufacturers, a range of recording parameters, and varying cable emplacement settings (e.g., shallow burial, borehole, subaqueous, dark fiber). Monitored cable lengths vary between 152 and 120129 m, with channel spacing between 1 and 49 m. The data has a total size of 6.8 TB, and is available for free download. Organizing and executing the Global DAS Month has produced a unique dataset for further exploration and highlighted areas of further development for the seismological community to address.

The cis- form of diaminodibenzocyclooctane (DADBCO, C16H18N2) is of interest as a negative coefficient of thermal expansion (CTE) material. The crystal structure was determined through single-crystal X-ray diffraction at 100 K and is presented herein.

Abstract not provided.

Monopulse is a technique for determining the Direction of Arrival (DOA) of a radar echo by comparing the simultaneous signal responses from two or more antenna beams or apertures. Two principal architectures are employed: 1) amplitude-comparison monopulse, and 2) phase-comparison monopulse. For a constrained-size fully and uniformly illuminated aperture, there is no meaningful difference between the DOA angle precision achievable by an amplitude monopulse architecture versus a phase monopulse

​Protocols play an essential role in Advance Reactor systems. A diverse set of protocols are available to these reactors. Advanced Reactors benefit from technologies that can minimize their resource utilization and costs. Evaluation frameworks are often used when assessing protocols and processes related to cryptographic security systems. The following report discusses the various characteristics associated with these protocol evaluation frameworks, and derives a novel evaluative framework.

The ShakeAlert Earthquake Early Warning (EEW) system aims to issue an advance warning to residents on the West Coast of the United States seconds before the ground shaking arrives, if the expected ground shaking exceeds a certain threshold. However, residents in tall buildings may experience much greater motion due to the dynamic response of the buildings. Therefore, there is an ongoing effort to extend ShakeAlert to include the contribution of building response to provide a more accurate estimation of the expected shaking intensity for tall buildings. Currently, the supposedly ideal solution of analyzing detailed finite element models of buildings under predicted ground-motion time histories is not theoretically or practically feasible. The authors have recently investigated existing simple methods to estimate peak floor acceleration (PFA) and determined these simple formulas are not practically suitable. Instead, this article explores another approach by extending the Pacific Earthquake Engineering Research Center (PEER) performance-based earthquake engineering (PBEE) to EEW, considering that every component involved in building response prediction is uncertain in the EEW scenario. While this idea is not new and has been proposed by other researchers, it has two shortcomings: (1) the simple beam model used for response prediction is prone to modeling uncertainty, which has not been quantified, and (2) the ground motions used for probabilistic demand models are not suitable for EEW applications. In this article, we address these two issues by incorporating modeling errors into the parameters of the beam model and using a new set of ground motions, respectively. We demonstrate how this approach could practically work using data from a 52-story building in downtown Los Angeles. Using the criteria and thresholds employed by previous researchers, we show that if peak ground acceleration (PGA) is accurately estimated, this approach can predict the expected level of human comfort in tall buildings.

Quantitative risk assessment (QRA) is highly dependent on data, leading to more robust models as new and updated data is acquired. The Hydrogen Plus Other Alternative Fuels Risk Assessment (HyRAM+) QRA capabilities include calculations of individual risk from leaks in a gaseous hydrogen facility due to the potential effects of jet fires and explosions. Leak frequencies are acquired through statistical analysis of published data from a variety of sources and industries. The filter leak frequencies in previous versions of the HyRAM+ software are substantially greater than the leak frequencies of other components, leading to QRA results for gaseous hydrogen in which filters consistently dominate the overall risk. Data that were previously used to derive the filter leak frequencies were reevaluated for applicability and additional data points were added to update the filter leak frequencies. The new frequencies are more comparable to leak frequencies for other components.

Abstract not provided.

This paper presents a new approach for autonomous motion planning for aircraft suffering from a loss-of-thrust emergency. Specifically, we show how modifications to the Closed-Loop Rapidly exploring Random Trees (CL-RRT) framework combined with controlled energy dissipation can enable rapid and effective kinodynamic motion planning. This CL-RRT Glide algorithm uses closed-loop prediction not only for node connections but also to estimate the remaining energy and prune infeasible paths. This greatly speeds up the search process, which is essential for emergency situations. In addition, we improve the ability of the gliding aircraft to reach a goal position and energy state. We do so by creating a Dissipative Total Energy Control Scheme (TECS). Dissipative TECS enables the glider to lose excess altitude in order to reach a desired energy level. Simulation results illustrate how the proposed methods enable faster motion planning. We also integrate the system into a small unmanned aerial vehicle system and experimentally demonstrate autonomous glide planning and execution during a motor-failure event. This type of algorithm can primarily benefit unmanned aircraft but can also serve to assist pilots in stressful emergency situations.

Photonic crystal nanobeam cavities are valued for their small mode volume, CMOS compatibility, and high coupling efficiency-crucial features for various low-power photonic applications and quantum information processing. However, despite their potential, nanobeam cavities often suffer from low quality factors due to fabrication imperfections that create surface states and optical absorption. In this work, we demonstrate InP nanobeam cavities with up to 140% higher quality factors by applying a coating of Al₂O₃ via atomic layer deposition to terminate dangling bonds and reduce surface absorption. Additionally, changing the deposition thickness allows precise tuning of the cavity mode wavelength without compromising the quality factor. This Al₂O₃ atomic layer deposition approach holds great promise for optimizing nanobeam cavities that are well-suited for integration with a wide range of photonic applications.

Abstract not provided.

All freely available plane-of-array (POA) transposition models and photovoltaic (PV) temperature and performance models in pvlib-python and pvpltools-python were examined against multiyear field data from Albuquerque, New Mexico. The data include different PV systems composed of crystalline silicon modules that vary in cell type, module construction, and materials. These systems have been characterized via IEC 61853-1 and 61853-2 testing, and the input data for each model were sourced from these system-specific test results, rather than considering any generic input data (e.g., manufacturer's specification [spec] sheets or generic Panneau Solaire [PAN] files). Six POA transposition models, 7 temperature models, and 12 performance models are included in this comparative analysis. These freely available models were proven effective across many different types of technologies. The POA transposition models exhibited average normalized mean bias errors (NMBEs) within ±3%. Most PV temperature models underestimated temperature exhibiting mean and median residuals ranging from −6.5°C to 2.7°C; all temperature models saw a reduction in root mean square error when using transient assumptions over steady state. The performance models demonstrated similar behavior with a first and third interquartile NMBEs within ±4.2% and an overall average NMBE within ±2.3%. Although differences among models were observed at different times of the day/year, this study shows that the availability of system-specific input data is more important than model selection. For example, using spec sheet or generic PAN file data with a complex PV performance model does not guarantee a better accuracy than a simpler PV performance model that uses system-specific data.

It is commonly assumed that cleaning photovoltaic (PV) modules is unnecessary when the inverter is undersized because clipping will sufficiently mask the soiling losses. Clipping occurs when the inverter's AC size is smaller than the overall modules' DC capacity and leads to the conversion of only part of the PV-generated DC energy into AC. This study evaluates the validity of this assumption, theoretically investigating the current magnitude of clipping and its effect on soiling over the contiguous United States. This is done by modelling energy yield, clipping and soiling across a grid of locations. The results show that in reality, under the current deployment trends, inverter undersizing minimally affects soiling, as it reduces these losses by no more than 1%absolute. Indeed, clipping masks soiling in areas where losses are already low, whereas it has a negligible effect where soiling is most significant. However, the mitigation effects might increase under conditions of lower performance losses or more pronounced inverter undersizing. In any case, one should take into account that degradation makes clipping less frequent as systems age, also decreasing its masking effect on soiling. Therefore, even if soiling was initially mitigated by the inverter undersizing, its effect would become more visible with time.

The Example Problems Manual supplements the User's Manual and the Theory Manual. The goal of the Example Problems Manual is to reduce learning time for complex end to end analyses. These documents are intended to be used together. See the User's Manual for a complete list of the options for a solution case. All the examples are part of the salinas test suite. Each runs as is.

Recent experimental findings have shown that tantalum single crystals display strong anisotropy during Taylor impact testing in stark contrast to isotropic deformation in polycrystalline counterparts. In this study, a coupled dislocation dynamics and finite element model was developed to simulate the complex stress field under dynamic loading of a Taylor impact test and track the intricate evolution of the dislocation microstructure. Our model allowed us to investigate detailed motion of dislocations and their mutual interactions and the effect of varying simulation parameters, such as sample size, initial dislocation density, crystallographic orientation, and temperature. Simulation results show good agreement with experimental observations and shed light on the mechanical response at small-scale under extreme loading conditions. In addition, resolved shear stress analysis incorporating the effect of shear stress from impact was performed to quantitatively support and provide a means to understand the model predictions of the impact foot shape.

We report the results of the second charged-particle transport coefficient code comparison workshop, which was held in Livermore, California on 24-27 July 2023. This workshop gathered theoretical, computational, and experimental scientists to assess the state of computational and experimental techniques for understanding charged-particle transport coefficients relevant to high-energy-density plasma science. Data for electronic and ionic transport coefficients, namely, the direct current electrical conductivity, electron thermal conductivity, ion shear viscosity, and ion thermal conductivity were computed and compared for multiple plasma conditions. Additional comparisons were carried out for electron-ion properties such as the electron-ion equilibration time and alpha particle stopping power. Overall, 39 participants submitted calculated results from 18 independent approaches, spanning methods from parameterized semi-empirical models to time-dependent density functional theory. In the cases studied here, we find significant differences—several orders of magnitude—between approaches, particularly at lower temperatures, and smaller differences—roughly a factor of five—among first-principles models. We investigate the origins of these differences through comparisons of underlying predictions of ionic and electronic structure. The results of this workshop help to identify plasma conditions where computationally inexpensive approaches are accurate, where computationally expensive models are required, and where experimental measurements will have high impact.

This report documents analysis to determine whether a hydrogen jet flame impinging on a tunnel ceiling structure could result in permanent damage to the Callahan tunnel in Boston, Massachusetts. This tunnel ceiling structure consists of a passive fire protective board supported by stainless steel hangers anchored to the tunnel ceiling with epoxy. Three types of fire protective boards were considered to determine whether heat from the flame could reach the stainless-steel hangers and the epoxy and cause the ceiling structure to collapse. Heat transfer analyses performed showed that the temperature remains constant where the steel hangers are attached to the passive fire protective board. According to these results, the passive fire protective board should provide adequate protection to the tunnel structure in this release scenario. Tunnel structures with similar suspended fire-resistant liner board materials should protect the integrity of the structure against the extremely low probability of an impinging hydrogen jet flame.

Abstract not provided.

High-fidelity simulations are performed to characterize the turbulence-induced wall pressure fluctuations on a sharp cone at a 5.5-degree angle-of-attack in a Mach 8 flow. Wall-resolved large-eddy simulation (LES) and wall-modeled large-eddy simulation (WMLES) results are compared to measurements at several locations on the cone body. Simulations are also compared to each other, and WMLES show good comparison in the autospectra, but modest comparison in the coherence.

This document is intended to help users program the new mid-circuit measurement (MCM) and classical branching capabilities of the Quantum Scientific Computing Open User Testbed (QSCOUT). Here, we present and explain an exemplar “ping-pong teleportation” program that makes repeated MCM and branching calls. The program is written in Jaqal, the quantum assembly language used by QSCOUT. This document is intended to accompany a companion Jupyter notebook Exemplar_one_bit_teleportation_pingpong.ipynb.

This report provides technical guidance for the calibration of laboratory glassware to help the practitioner achieve traceability to the International System of Units and meet customer quality requirements. The discussion of traceability uses the National Institute of Standards and Technology’s seven essential elements of traceability as a framework. The guidance also includes how to determine when calibration is necessary, practical tips, and helpful references.

Current biogeochemical models produce carbon–climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first-order or Michaelis–Menten kinetics at the global scale. Nevertheless, a wider range of data with high-quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics-function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

NasGen provides a path for migration of structural models from NASTRAN bulk data format (BDF) into both an Exodus mesh file and an ASCII input file for Sierra Structural Dynamics (Salinas) and Solid Mechanics (Presto).