Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic temperature”

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 1,063 records · Page 59

Disordered interfaces of alkaline aluminate salt hydrates provide glimpses of Al 3+ coordination changes

Hypothesis: The precipitation and dissolution of aluminum-bearing mineral phases in aqueous systems often proceed via changes in both aluminum coordination number and connectivity, complicating molecular-scale interpretation of the transformation mechanism. Here, the thermally induced transformation of crystalline sodium aluminum salt hydrate, a phase comprised of monomeric octahedrally coordinated aluminate which is of relevance to industrial aluminum processing, has been studied. Because intermediate aluminum coordination states during melting have not previously been detected, it is hypothesized that the transition to lower coordinated aluminum ions occurs within a highly disordered quasi-two-dimensional phase at the solid-solution interface. Experiments and simulations: In this work, in situ X-ray diffraction (XRD), Raman and 27 Al nuclear magnetic resonance (NMR) spectroscopy were used to monitor the melting transition of nonasodium aluminate hydrate (NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6H 2 O). A mechanistic interpretation was developed based on complementary classical molecular dynamics (CMD) simulations including enhanced sampling. A reactive forcefield was developed to bridge speciation in the solution and in the solid phase. Findings: In contrast to classical dissolution, aluminum coordination change proceeds through a dynamically stabilized ensemble of intermediate states in a disordered layer at the solid-solution interface. In both melting and dissolution of NSA, octahedral, monomeric aluminum transition through an intermediate of pentahedral coordination. The intermediate dehydroxylates to form tetrahedral aluminate (Al(OH)$_4^–$) in the liquid phase. This coordination change is concomitant with a breaking of the ionic aluminate-sodium ion linkages. The solution phase Al(OH)$_4^–$ ions subsequently polymerize into polynuclear aluminate ions. However, there are some differences between bulk melting and interfacial dissolution, with the onset of the surface-controlled process occurring at a lower temperature (~30 °C) and the coordination change taking place more gradually as a function of temperature. This work to determine the local structure and dynamics of aluminum in the disordered layer provides a new basis to understand mechanisms controlling aluminum phase transformations in highly alkaline solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating structure and dynamics of unentangled poly(dimethyl- co -diphenyl)siloxane via molecular dynamics simulation

Polysiloxane is one of the most important polymeric materials in technological use. Polydimethylsiloxane displays glass-like mechanical properties at low temperatures. Incorporation of phenyl siloxane, via copolymerization for example, improves not only the low-temperature elasticity but also enhances its performance over a wide range of temperatures. Copolymerization with the phenyl component can significantly change the microscopic properties of polysiloxanes, such as chain dynamics and relaxation. However, despite much work in the literature, the influence of such changes is still not clearly understood. In this work, we systematically study the structure and dynamics of random poly(dimethyl- co -diphenyl)siloxane via atomistic molecular dynamics simulations. As the molar ratio Φ of the diphenyl component increases, we find that the size of the linear copolymer chain expands. At the same time, the chain-diffusivity slows down by over an order of magnitudes. In conclusion, the reduced diffusivity appears to be a result of a complex interplay between the structural and dynamic changes induced by phenyl substitution.

36 MATERIALS SCIENCE↗

Data from: "Moisture rivals temperature in limiting photosynthesis by trees establishing beyond their cold-edge range limit under ambient and warmed conditions"

This archive contains data files that were used to draw conclusions in “Moisture rivals temperature in limiting photosynthesis by trees establishing beyond their cold-edge range limit under ambient and warmed conditions”, by Moyes et al., 2015. All field research was completed in common garden plots set up as part of the Alpine Treeline Warming Experiment (ATWE) on Niwot Ridge, Colorado, USA.There are two main data file formats in this archive: comma-separated values (.csv), and Microsoft Excel (.xls and .xlsx). .xlsx files can be read using Microsoft Excel and Google Sheets, and .csv files can be read using any simple text editor program, such as TextEdit (Mac) and Notepad (Windows). This .pdf data user’s guide can be read using Adobe Acrobat Reader, or any other compatible software. Seedling photographs and their corresponding leaf area-processed images are available in .jpg/.JPG image format, and can be opened using Preview (Mac) and Photos (Windows). To provide additional spatial context, two types of geospatial files are also published in this data package: ESRI shapefiles (.shp) and .kml files. Shapefiles are compatible with any GIS software able to read the file type (such as QGIS or ESRI’s ArcGIS suite), and .kml files can be opened with Google Earth or Google Maps. Figures 3 and 4 in the publication contain data from Moyes et al. 2013. This publication is cited in the References section in this archive, and data files can be accessed via the Alpine Treeline Warming Experiment project portal on ESS-DIVE. ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- Climate change is altering plant species distributions globally, and warming is expected to promote uphill shifts in mountain trees. However, at many cold-edge range limits, such as alpine treelines in the western United States, tree establishment may be colimited by low temperature and low moisture, making recruitment patterns with warming difficult to predict.- We measured response functions linking carbon (C) assimilation and temperature- and moisture-related microclimatic factors for limber pine (Pinus flexilis) seedlings growing in a heating × watering experiment within and above the alpine treeline. We then extrapolated these response functions using observed microclimate conditions to estimate the net effects of warming and associated soil drying on C assimilation across an entire growing season.- Moisture and temperature limitations were each estimated to reduce potential growing season C gain from a theoretical upper limit by 15–30% (c. 50% combined). Warming above current treeline conditions provided relatively little benefit to modeled net assimilation, whereas assimilation was sensitive to either wetter or drier conditions.- Summer precipitation may be at least as important as temperature in constraining C gain by establishing subalpine trees at and above current alpine treelines as seasonally dry subalpine and alpine ecosystems continue to warm.

54 ENVIRONMENTAL SCIENCES↗

Metal-insulator and magnetic phase diagram of Ca 2 RuO 4 from auxiliary field quantum Monte Carlo and dynamical mean field theory

Layered perovskite ruthenium oxides exhibit a striking series of metal-insulator and magnetic-nonmagnetic phase transitions easily tuned by temperature, pressure, epitaxy, and nonlinear drive. In this work, we combine results from two complementary state-of-the-art many-body methods, auxiliary field quantum Monte Carlo and dynamical mean field theory, to determine the low-temperature phase diagram of Ca 2 RuO 4 . Both methods predict a low-temperature, pressure-driven metal-insulator transition accompanied by a ferromagnetic-antiferromagnetic transition. Further, the properties of the ferromagnetic state vary nonmonotonically with pressure and are dominated by the ruthenium $d_{xy}$ orbital, while the properties of the antiferromagnetic state are dominated by the $d_{xz}$ and $d_{yz}$ orbitals. Differences in the details of the predictions of the two methods are analyzed. This work is theoretically important as it presents the first application of the auxiliary field quantum Monte Carlo method to an orbitally degenerate system with both Mott and Hunds physics and provides an important comparison of the dynamical mean field and auxiliary field quantum Monte Carlo methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic Heterogeneity and Kovacs’ Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. Here, we used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water’s structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wafer-Scale Synthesis of 2D Materials by an Amorphous Phase-Mediated Crystallization Approach

The interest in the wafer-scale growth of two-dimensional (2D) materials, including transition metal dichalcogenides (TMDCs), has been rising for transitioning from lab-scale devices to commercial-scale systems. Among various synthesis techniques, physical vapor deposition, such as pulsed laser deposition (PLD), has shown promise for the wafer-scale growth of 2D materials. However, due to the high volatility of chalcogen atoms (e.g., S and Se), films deposited by PLD usually suffer from a lack of stoichiometry and chalcogen deficiency. To mitigate this issue, excess chalcogen is necessary during the deposition, which results in problems like uniformity or not being repeatable. This study demonstrates a condensed-phase or amorphous phase-mediated crystallization (APMC) approach for the wafer-scale synthesis of 2D materials. This method uses a room-temperature PLD process for the deposition and formation of amorphous precursors with controlled thicknesses, followed by a post-deposition crystallization process to convert the amorphous materials to crystalline structures. This approach maintains the stoichiometry of the deposited materials throughout the deposition and crystallization process and enables the large-scale synthesis of crystalline 2D materials (e.g., MoS 2 and WSe 2 ) on Si/SiO 2 substrates, which is critical for future wafer-scale electronics. We show that the thickness of the layers can be digitally controlled by the number of laser pulses during the PLD phase. Optical spectroscopy is used to monitor the crystallization dynamics of amorphous layers as a function of annealing temperature. The crystalline quality, domain sizes, and the number of layers were explored using nanoscale and atomistic characterization (e.g., AFM, STEM, and EDS) along with electrical characterization to explore process–structure–performance relationships. Finally, this growth technique is a promising method that could potentially be adopted in conventional semiconductor industries for wafer-scale manufacturing of next-generation electronic and optoelectronic devices.

2D electronics↗

Influence of geochemical reactions on the creep behavior of Mt. Simon sandstone

The capture and subsurface storage of carbon dioxide is a sustainable option that is currently pursued worldwide to mitigate greenhouse gas effect. However, predicting the long-term mechanical integrity of CO 2 underground storage systems remains a challenge. To address that question, it is essential to understand the influence of fluid-rock chemo-mechanical interactions on the long-term and on the time-dependent mechanical properties. In turn, the long-term mechanical response and the time-dependent mechanical behavior can be represented by the creep response. We investigate the impact of CO 2 -induced geochemical reaction on the creep response of Mt. Simon sandstone with a 50–400 micron grain size. We perform static and dynamic flow experiments on Mt. Simon sandstone specimens under geological conditions, at a temperature of T = 50 - 53 °C and for CO 2 pressures of P = 8:62; 17:2 MPa under both static flow and dynamic flow-through conditions. After aging, we employ creep indentation testing, high-resolution SEM-EDS, computer vision, machine learning, and micromechanics modeling to probe changes on the microstructure and mechanical properties. Following both static and dynamic flow-through experiments, we observe a 10–22% decrease in quartz volume fraction and an increase in both the microporosity (7–28%) and nanoporosity (60–65%). Additional CO 2 -induced microstructural changes include an enlargement of pore throats and the formation of channels. These observations point to the presence of CO 2 -induced K-feldspar dissolution and clay dissolution reactions. Here, the macroscopic creep behavior is logarithmic and the macroscopic creep modulus varies depending on the microporosity and the relative quartz and feldspar content. As a result of these geochemical reactions and of the related microstructural changes, a 55-60% decrease in the macroscopic logarithmic creep modulus is predicted.

58 GEOSCIENCES↗

A pathway to synthesizing single-crystal Fe and FeCr films

Nuclear reactor environments provide a unique scientific and engineering challenge wherein materials must tolerate prolonged exposure to concurrent irradiation, elevated temperatures, and corrosive media. However, uncontrolled variability in material composition and structure often prohibits truly single-variable experiments that can reveal basic aspects of environmental damage. Here, Magnetron sputtering is used to provide a more controlled model system for these fundamental studies, yielding reproducible single-crystal Fe and FeCr thin films containing 8 and 18 at.% Cr. Electron microscopy is used to determine the systematic correlations between growth conditions and the resulting film microstructure and surface morphology. It is found that the substrate temperature and applied radio frequency (RF) bias can be tuned to obtain consistent homogeneous and single crystal films with a minimal amount of Ar impurities from the RF bias process. Epitaxial, single-crystal Fe films are obtained on MgO substrates at 500 °C with 10 Watt (W) RF bias deposition. However, when Cr is alloyed with Fe, higher substrate temperatures (600 °C) and applied RF biases (15 W) are required to achieve a similar epitaxial single-crystal FeCr film. Accelerated molecular dynamics simulations reveal that Cr impedes surface transport, explaining the need for higher temperature and bias during the growth of the Cr-bearing films.

36 MATERIALS SCIENCE↗

Kuramoto synchronization of quantum tunneling polarons for describing the dynamic structure in cuprate superconductors

A major open topic in cuprates is the interplay between the lattice and electronic dynamics and the importance of their coupling to the mechanism of high-temperature superconductivity (HTSC). As evidenced by extended x-ray absorption fine structure (EXAFS) experiments, anharmonic structural effects are correlated with the charge dynamics and the transition to a superconducting phase in different HTSC compounds. Here, in this paper, we describe how structural anharmonic effects can be coupled to electronic and lattice dynamics in cuprate systems by performing the exact diagonalization of a prototype anharmonic many-body Hamiltonian on a relevant 6-atom cluster and show that the EXAFS results can be understood as a Kuramoto synchronization transition between coupled internal quantum tunneling of polarons associated with the two-site distribution of the copper–apical oxygen (Cu-O ap ) pair in the dynamic structure. The transition is driven by the anharmonicity of the lattice vibrations and promotes the pumping of charge, initially stored at the apical oxygen reservoirs, into the copper-oxide plane. Simultaneously, a finite projection of the internal quantum tunneling polaron extends to the copper–planar oxygen (Cu-O pl ) pair. All these findings allow an interpretation based on an effective quantum-mechanical triple-well potential associated with the oxygen sites of the 6-atom cluster, which accurately represents the phase synchronization of apical oxygens and lattice-assisted charge transfer to the CuO 2 plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Theory of finite-temperature two-dimensional structural transformations in group-IV monochalcogenide monolayers

An account of two-dimensional (2D) structural transformations in 2D ferroelectrics predicts an evolution from a structure with Pnm 2 1 symmetry into a structure with square P4/nmm symmetry and is consistent with experimental evidence, while another argues for a transformation into a structure with rectangular Pnmm symmetry. An analysis of the assumptions made in these models is provided here, and six fundamental results concerning these transformations are contributed as follows: (i) Softened phonon modes produce rotational modes in these materials. (ii) The transformation to a structure with P4/nmm symmetry occurs at the lowest critical temperature T c . (iii) The hypothesis that one unidirectional optical vibrational mode underpins the 2D transformation is unwarranted. (iv) Being successively more constrained, a succession of critical temperatures ( T c < T c ' < T c " ) occurs in going from molecular dynamics calculations with the NPT and NVT ensembles onto models with unidirectional oscillations. (v) The choice of exchange-correlation functional impacts the estimate of the critical temperature. (vi) Importantly, the correct physical picture of these transformations is one in which rotational modes confer a topological character to the 2D transformation via the proliferation of vortices.

2d ferroelectrics↗

Energy migration and scintillation kinetics in compositionally complex (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce single crystal scintillator

It is well-established that compositional tuning through binary admixture can improve scintillation performance in several materials systems, including Ce-activated garnets. Although recent work on ternary or quaternary cation admixture shows promise, the impact of this increased compositional complexity on thermal stability and carrier-defect dynamics has not been addressed. Here, we investigate a compositionally complex garnet, (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce (GYTLAG), grown by the Czochralski method using temperature-dependent photoluminescence (PL), PL decay, and thermoluminescence (TL). PL and PL decay measurements support a thermally activated Tb 3+ -Ce 3+ energy transfer, where Tb 3+ emission dominates below 60 K, but Ce 3+ emission increases from 20-300 K. Thermal quenching of Ce 3+ emission occurs around T 50 = 508 K, with an activation energy of 0.6 eV. TL and wavelength-resolved TL spectra from 20–500 K show that GYTLAG contains similar trap groups to LuAG but with a broader distribution of glow peaks below room temperature, possibly caused by quaternary cation mixing. A combination of dose dependence, partial cleaning and initial rise, and glow curve fitting to a first order continuous Gaussian distribution model are used to understand the contribution of electronic point defects to scintillation decay and afterglow at room temperature. Furthermore, these results inform how increased compositional complexity influences recombination dynamics in garnet scintillators.

Compositionally complex↗

Analysis of Frost Prevention for Heat Recovery Ventilation in an Arctic Facility

Sub-Arctic environments present a unique challenge for the continuous operation of mechanical ventilation systems in occupied buildings where cold temperatures create complications for maintaining adequate air quality, efficiency, and system operability. Heat recovery ventilation is most difficult at the coldest temperatures, when there is the greatest availability to transfer heat between airstreams, due to the potential for frost accumulation in the exhaust airstream. Frost limits heat transfer and creates the potential for subsequent blockage of airflow. A successful frost prevention strategy is necessary to ensure that recovery ventilation systems function in locations with a large amount of annual heating degree days. Commonly, this is achieved by preheating outdoor air above freezing temperatures via airflow through a glycol-water mix hydronic coil before heat exchange occurs. Preheating the outdoor air proactively inhibits frost formation but is associated with a penalty to recovery efficiency. A commercial-scale heat recovery ventilator in a facility in Fairbanks, Alaska was analyzed for sensible recovery effectiveness and exhaust temperature at different winter conditions. The system utilized the common defrost strategy of preheating the air above-freezing with a static temperature control scheme. The heat transfer efficiency of the recovery ventilator was observed to be poor at various winter conditions with elevated exhaust temperatures indicating lost heat recovery effectiveness. Using prior findings on strategies for sub-freezing temperatures along with observed field conditions it was determined that a dynamic relationship between the coil and the exterior conditions would reduce the energy required to preheat the air and allow the heat exchanger to maximize effectiveness due to temperate difference. This reduction in preheated temperature of the air is achieved by modulating the temperature of the coil to keep the discharge air stream above a certain dewpoint temperature threshold, likely leading to energy savings for other sub-Arctic facilities using this frost control strategy.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗

Dynamic Simulation Modeling and Control of a Desiccant Assisted Direct-expansion Air Handling Unit

Desirable built environments demand simultaneous regulation of thermal comfort and indoor air quality (IAQ) with energy-efficient operation of heating, ventilation and air conditioning (HV AC) systems, which involves controls of temperature, humidity and airborne contaminants simultaneously. This paper presents the efforts of dynamic modeling and initial development control strategy for a desiccant-assisted multi-functional air handling unit (AHU) coupled with a direct-expansion rooftop unit (RTU) system, which aims to achieve multiple functions for indoor environment conditioning with energy efficient control. The RTU-AHU system includes a desiccant wheel for dehumidification and a conceptual direct air capture (DAC) filtering device for CO2 regulation. A Modelica-based dynamic model is developed for this conceptual system, and a simple decentralized control strategy is designed, which combines a differential-enthalpy based AHU return-air ratio control, a demand-controlled ventilation, and supply-air temperature humidity control via the RTU and DW controls. The proposed control method is evaluated with the Modelica simulation model for a selected set of scenarios

Pan, Chao↗

Expansion dynamics of femtosecond laser-induced plasmas: Influence of thermophysical plasma properties

This study investigates the expansion dynamics of femtosecond laser-induced plasmas, emphasizing the impact of plasma thermophysical properties and ambient gas composition. Through shadowgraphy experiments and multiphase computational fluid dynamics (CFD) simulations, the influence of parameters such as heat capacity, molecular weight, and thermal conductivity on plume morphology, shockwave evolution, and energy dissipation mechanisms is examined. A mixture multiphase model is implemented to capture the interaction between the plasma and the surrounding gas. Simulation results reveal that plasma expansion is strongly inertia-driven. Results show that differences in plasma properties and ambient conditions affect the shape and temperature distribution of the expanding plume. The early-stage dynamics are primarily dictated by pressure forces, whereas thermal and viscous effects play a growing role in the plume's behavior during later stages of expansion. The CFD findings show the necessity of accurate initial condition characterization, including crater geometry and plasma pressure and temperature, for reliable modeling of plasma evolution in laser ablation processes.

CFD modeling↗

Tailoring the Local Environment of Platinum in Single‐Atom Pt 1 /CeO 2 Catalysts for Robust Low‐Temperature CO Oxidation

Abstract A single‐atom Pt 1 /CeO 2 catalyst formed by atom trapping (AT, 800 °C in air) shows excellent thermal stability but is inactive for CO oxidation at low temperatures owing to over‐stabilization of Pt 2+ in a highly symmetric square‐planar Pt 1 O 4 coordination environment. Reductive activation to form Pt nanoparticles (NPs) results in enhanced activity; however, the NPs are easily oxidized, leading to drastic activity loss. Herein we show that tailoring the local environment of isolated Pt 2+ by thermal‐shock (TS) synthesis leads to a highly active and thermally stable Pt 1 /CeO 2 catalyst. Ultrafast shockwaves (>1200 °C) in an inert atmosphere induced surface reconstruction of CeO 2 to generate Pt single atoms in an asymmetric Pt 1 O 4 configuration. Owing to this unique coordination, Pt 1 δ+ in a partially reduced state dynamically evolves during CO oxidation, resulting in exceptional low‐temperature performance. CO oxidation reactivity on the Pt 1 /CeO 2 _TS catalyst was retained under oxidizing conditions.

Jiang, Dong↗

Computational fluid dynamics investigations of flow, heat transfer, and oxidation in heat recovery steam generator

Modern heat recovery steam generators (HRSGs) operate at elevated temperatures, leading to the formation of oxides inside the tubes of heat exchangers (HXs). This oxide growth reduces the heat recovery efficiency. Moreover, after reaching a certain critical thickness, some oxide scales detach from the tube surface (exfoliation), causing erosion damage to the components downstream. Predicting the metal temperature distribution and associated oxide thickness in the HXs of an HRSG can aid in mitigating these problems. Here, a computational fluid dynamics (CFD) model was developed within the commercial code STAR-CCM+ for the prediction of fluid flow, conjugate heat transfer, and associated oxidation in HRSGs. Moreover, a new Porous Media Model (PMM) method was developed to model the fin effect on the heat transfer in HX, which can substantially reduce the prohibitive computational costs of fin meshes. The developed CFD model was used to conduct a high-fidelity simulation of a real-scale HRSG to investigate flow, heat transfer, and oxide growth. The calculated oxide thickness on different tubes can be used to identify HX regions that require oxide-resistant coatings to prevent exfoliation and ensuing damages. Furthermore, this CFD framework can serve as a reference for future studies that intend to model and investigate high-temperature oxidation in HXs used for any applications.

42 ENGINEERING↗

Collapse of the waiting time effect in a spin glass

Waiting time measurements are reported here for a bulk polycrystalline sample of the spin glass Cu 0.95 Mn 0.05 . We observe a suppression of aging as the glass temperature is approached from below. These results are interpreted in terms of a spin-glass correlation length approaching and reaching crystallite sizes. They provide a direct experimental connection between the correlation length scale and relaxation times in the spin-glass state. In particular, these measurements track the value of the dynamic critical exponent as the system transitions from the critical temperature regime to its low-temperature, quasistatic nonequilibrium phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗