Search NASA⌕ Search

SEARCH · Search NASA

Results for “ARS”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Future Warming and Intensification of Precipitation Extremes: A “Double Whammy” Leading to Increasing Flood Risk in California

This study focuses on quantifying future anthropogenic changes in surface runoff associated with extreme precipitation in California's Sierra Nevada. The method involves driving a land surface model with output from a high resolution regional atmospheric simulation of the most extreme atmospheric rivers (ARs). AR events were selected from an ensemble of global climate model simulations of historical and late 21st century climate under the “high-emission” RCP8.5 scenario. Average precipitation during the future ARs increases by ~25% but a much lower proportion falls as snow. The resulting future runoff increase is dramatic—nearly 50%, reflecting both the precipitation increase and simultaneous conversion of snow to rain. The “double whammy” impact on runoff is largest in the 2,000–2,500 m elevation band, where the snowfall loss and precipitation increase are both especially large. Overall, this huge increase in runoff during the most extreme AR events could present major flood control challenges for the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Atmospheric Rivers Impacting Northern California Exhibit a Quasi-Decadal Frequency

Periods of water surplus and deficit in Northern California follow a pronounced quasi-decadal cycle. This cycle is largely driven by the frequency of atmospheric rivers (ARs), affecting the region’s wet and dry periods. Our analyses demonstrate that the quasi-decadal cycle of AR frequency relies on moisture transport associated with the position and intensity of the Aleutian Low. In observations, the Aleutian Low is shown to covary with tropical Pacific sea surface temperature anomalies. A modeling experiment, which incorporates ocean observations from the equatorial Pacific into the fully coupled climate model, provides support that the quasi-decadal cycle of the Aleutian Low is forced by the tropical Pacific. Subsequently, the tropical Pacific modulates the wet season moisture transport toward California on decadal time scales, affecting AR frequency. These results provide metrics for improving interannual-to-decadal prediction of AR activity, which drives hydrological cycles in Northern California.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Atmospheric River Response to Arctic Sea Ice Loss in the Polar Amplification Model Intercomparison Project

The atmospheric river (AR) response to Arctic sea ice loss in the Northern hemisphere winter is investigated using simulations from the Polar Amplification Model Intercomparison Project. Results have shown that the midlatitude responses are dominated by dynamic effects. Poleward of around 60°N, the dynamic and thermodynamic effects cancel each other, resulting in relatively small responses. The response uncertainty can be characterized by leading uncertainty modes, with the responses over the Pacific and Atlantic projecting onto the northeastward extension and equatorward shift mode, respectively. In addition, the responses seem to be mean state-dependent: under the same forcing, models with more poleward-located climatological ARs tend to show stronger equatorward shifts over the Atlantic; over the Pacific, models with more westward-located climatological AR core tend to show stronger northeastward extensions. Finally, these relationships highlight the importance of improving the AR climatology representation on reducing the response uncertainty to Arctic sea ice loss.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Radiation exposure induces cross-species temporal metabolic changes that are mitigated in mice by amifostine

Exposure to acute, damaging radiation may occur through a variety of events from cancer therapy and industrial accidents to terrorist attacks and military actions. Our understanding of how to protect individuals and mitigate the effects of radiation injury or Acute Radiation Syndrome (ARS) is still limited. There are only a few Food and Drug Administration-approved therapies for ARS; whereas, amifostine is limited to treating low dose (0.7–6 Gy) radiation poisoning arising from cancer radiotherapy. An early intervention is critical to treat ARS, which necessitates identifying diagnostic biomarkers to quickly characterize radiation exposure. Towards this end, a multiplatform metabolomics study was performed to comprehensively characterize the temporal changes in metabolite levels from mice and non-human primate serum samples following γ-irradiation. The metabolomic signature of amifostine was also evaluated in mice as a model for radioprotection. The NMR and mass spectrometry metabolomics analysis identified 23 dysregulated pathways resulting from the radiation exposure. These metabolomic alterations exhibited distinct trajectories within glucose metabolism, phospholipid biosynthesis, and nucleotide metabolism. A return to baseline levels with amifostine treatment occurred for these pathways within a week of radiation exposure. Together, our data suggests a unique physiological change that is independent of radiation dose or species. Furthermore, a metabolic signature of radioprotection was observed through the use of amifostine prophylaxis of ARS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Formation of tungsten ethylidene complexes from diethyl complexes through a proton-catalyzed rearrangement of ethylene

Here, W(NAr) 2 Et 2 (Ar = 2,6-diisopropylphenyl) reacts with two equivalents of R F9 OH (OR F9 = OC(CF 3 ) 3 ) to yield W(NAr)(ArNH 2 )(OR F9 )2(C 2 H 4 ) complexes and ethane. In solution W(NAr)(ArNH 2 )(OR F9 ) 2 (C 2 H 4 ) decomposes to give R F9 OH, ethane, W(NAr)(OR F9 ) 2 (C 2 H 4 ), and W(NAr)(NHAr′)(ArNH 2 )(OR F9 ), in which Ar′ contains a dehydrogenated isopropyl group (Ar′ = (2-i-Pr)(6-CMe=CH 2 )C 6 H 3 ) coordinated to the metal. On a similar time scale W(NAr)(OR F9 ) 2 (CHCH 3 ) complexes are formed from W(NAr)(OR F9 ) 2 (C 2 H 4 ) through an ArNH 2 -catalyzed rearrangement of the ethylene ligand. W(NAr)(NHAr′)(ArNH 2 )(OR F9 ) reacts with cyclohexene to form methylenecyclohexene and complexes that contain an NHAr″ ligand where Ar″ is a disubstituted (methyl/aryl) alkylidene, (2-i-Pr)(6-CMe)C 6 H 3 ) that is tethered to the metal through the amido nitrogen. In contrast to W(NAr)(ArNH 2 )(OR F9 ) 2 (C 2 H 4 ), analogous OR F6 (OCMe(CF 3 ) 2 ) and OR F3 (OCMe 2 (CF 3 )) complexes are relatively stable at 22 °C.

Maji, Milan [University of California, Riverside, ↗
DOE OSTI2020

Observation of shock-front separation in multi-ion-species collisional plasma shocks

We observe shock-front separation and species-dependent shock widths in multi-ion-species collisional plasma shocks, which are produced by obliquely merging plasma jets of a He/Ar mixture (97% He and 3% Ar by initial number density) on the plasma liner experiment. Visible plasma emission near the He-I 587.6-nm and Ar-II 476.5–514.5-nm lines is simultaneously recorded by splitting a single visible image of the shock into two different fast-framing cameras with different narrow bandpass filters (589 ± 5 nm for observing the He-I line and 500 ± 25 nm for the Ar-II lines). For conditions in these experiments (pre-shock ion and electron densities ≈5×10 14 cm –3 , ion and electron temperatures of ≈2.2 eV, and relative plasma-merging speed of 22 km/s), the observationally inferred magnitude of He/Ar shock-front separation and the shock widths themselves are <1 cm, which correspond to ~50 post-shock thermal ion–ion mean free paths. These experimental length scales are in reasonable qualitative and quantitative agreement with results from 1D multi-fluid simulations using the CHICAGO code. However, we note there are differences between the experimentally inferred and simulation-predicted ionization states and line emission intensities, particularly in the post-shock region. Overall, the experimental and simulation results are consistent with theoretical predictions that the lighter He ions diffuse farther ahead within the overall shock front than the heavier Ar ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Decomposition of methane diluted with inert gas in an RF discharge cell

Decomposition of methane using non-thermal plasmas is an attractive route for producing hydrogen-rich gases and valuable carbon nanomaterials. Understanding how plasma discharge modes influence methane decomposition in optimizing plasma-assisted chemical conversion remains unexplored. This study explores the coupling between the discharge structure and product selectivity in RF capacitively coupled discharges operating in methane/inert gas mixtures in the pressure range of 2–3 torr. Under our experimental conditions the discharge exhibits mode transitions from uniform to striated in Ar and Kr and from diffuse to contracted in Ar and Kr with <5% CH 4 . The discharges in He and Ne remained uniform under our operating conditions, and their mixtures with CH 4 remained diffuse. A 0-d model for Ar/CH 4 discharge established a threshold for contraction while also asserting the importance of Ar$^{*}_{\textrm{m}}$ in the dissociation and ionization processes. The highest degree of methane decomposition, >99.7% with the main products of acetylene and graphitized solid carbon was achieved in the contracted discharge mode for both Kr or Ar with ⩽5% CH 4 . We demonstrate that contraction can play a crucial role in the effective decomposition of methane with value-added products and that both the electronic and thermal properties of plasma gas are responsible for this effect.

RF plasma↗
DOE OSTI2023

Trends in runaway electron plateau partial recombination by massive H 2 or D 2 injection in DIII-D and JET and first extrapolations to ITER and SPARC

Abstract Experimental trends in thermal plasma partial recombination resulting from massive D 2 injection into high-Z (Ar) containing runaway electron (RE) plateaus in DIII-D and JET are studied for the purpose of achieving sufficiently low electron density ( n e ≈ 10 18 m − 3 ) to increase RE final loss MHD levels. In both DIII-D and JET, thermal electron density n e is found to drop by ∼100 × when the thermal plasma partially recombines, with a minimum at a vacuum vessel-averaged D 2 density in the range 10 20 − 10 21 m − 3 . RE effective resistivity also drops after partial recombination, indicating expulsion of the Ar content. The n e level after partial recombination is found to increase as RE current is increased. The amount of initial Ar in the RE plateau is not observed to have a strong effect on partial recombination. Partial recombination timescales of order 5 ms in DIII-D and 15 ms in JET are observed. These basic trends and timescales are matched with a 1D diffusion model, which is then used to extrapolate to ITER and SPARC tokamaks. Within the approximations of this model, it is predicted that ITER will be able to achieve sufficiently low n e values on time scales faster than expected RE plateau vertical drift timescales (of order 100 ms), provided sufficient D 2 or H 2 is injected. In SPARC, it is predicted that achieving significant n e recombination will be challenging, due to the very high RE current density. In both ITER and SPARC, it is predicted that achieving low n e will be easier with Ar as a background impurity (rather than Ne).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

CUT&RUN identifies centromeric DNA regions of Rhodotorula toruloides IFO0880

ABSTRACT Rhodotorula toruloides has been increasingly explored as a host for bioproduction of lipids, fatty acid derivatives and terpenoids. Various genetic tools have been developed, but neither a centromere nor an autonomously replicating sequence (ARS), both necessary elements for stable episomal plasmid maintenance, has yet been reported. In this study, cleavage under targets and release using nuclease (CUT&RUN), a method used for genome-wide mapping of DNA–protein interactions, was used to identify R. toruloides IFO0880 genomic regions associated with the centromeric histone H3 protein Cse4, a marker of centromeric DNA. Fifteen putative centromeres ranging from 8 to 19 kb in length were identified and analyzed, and four were tested for, but did not show, ARS activity. These centromeric sequences contained below average GC content, corresponded to transcriptional cold spots, were primarily nonrepetitive and shared some vestigial transposon-related sequences but otherwise did not show significant sequence conservation. Future efforts to identify an ARS in this yeast can utilize these centromeric DNA sequences to improve the stability of episomal plasmids derived from putative ARS elements.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Evolving language of pediatric anxiety in electronic health records

Objectives This study aimed to identify and quantify semantic drift (ie, the change in semantic meaning over time) within expert-defined anxiety-related (AR) terminology and compare it to common electronic health record (EHR) vocabulary across longitudinal pediatric clinical notes. Materials and Methods A corpus of pediatric clinical notes from 2009 to 2022 was analyzed using computational methods. Semantic drift for each term was quantified using cosine similarity between annual temporal word embeddings. Contextual meaning was examined through changes in nearest neighbors across years. The Laws of Semantic Change were applied to assess the influence of word frequency and polysemy. Vocabulary terms were categorized as AR or common EHR. Results 98% of AR terminology maintained a cosine similarity between 0.00 and 0.50, indicating moderate semantic stability, whereas 90% of common EHR terms remained between 0.00 and 0.25, showing greater contextual stability overall. Frequent terms exhibited minimal change (Frequency Coefficient = 0.04), whereas highly polysemous or abbreviated terms showed less stability (Polysemy Coefficient = 0.630). AR terminology drifted more slowly than general EHR vocabulary (Type Coefficient = −0.179), further supported by significant year–type interactions (Coef = −0.09 to −0.523). Discussion Although anxiety-related terminology demonstrates slower semantic drift than general EHR vocabulary, subtle contextual shifts still occur that may affect downstream interpretability and retrieval in automated systems. Conclusion Continuous linguistic monitoring and adaptive modeling are essential to maintain semantic fidelity and ensure the long-term reliability of clinical decision support systems as healthcare documentation evolves.

Pediatric anxiety disorders↗
DOE OSTI2021

Nanoscale dynamics during self-organized ion beam patterning of Si. II. Kr + bombardment

Understanding the self-organized ion beam nanopatterning of elemental semiconductors, particularly silicon, is of intrinsic scientific and technological interest. This is the second component of a two-part coherent x-ray scattering and x-ray photon correlation spectroscopy (XPCS) investigation of the kinetics and fluctuation dynamics of nanoscale ripple development on silicon during 1 keV Ar + (part I) and Kr + bombardment at 65° polar angle. Here it is found that the ion-enhanced viscous flow relaxation is essentially equal for Kr + -induced patterning as previously found for Ar + patterning despite the difference in ion masses. However, the magnitude of the surface curvature-dependent roughening rate in the early-stage kinetics is larger for Kr + than for Ar + , consistent with expectations that the heavier ion gives an increased mass redistributive contribution to the initial surface instability. As with the Ar + case, fluctuation dynamics in the late stage show a peak in correlation times at the length scale corresponding to the dominant structural feature on the surface—the ripples. Finally, it is shown that speckle motion during the surface evolution can be analyzed to determine spatial inhomogeneities in erosion rate and ripple velocity. This allows the direction and speed of ripple motion to be measured in a real time experiment. In the present case, ripple motion is found to be into the projected direction of the ion source, in contrast to expectations from an existing sputter erosion driven model with parameters derived from binary collision approximation simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

New precise half-life measurement for the superallowed β+ emitter Ar34

We have measured the half-life of 34 Ar, the parent of a superallowed 0 + → 0 + β transition, to be 0.84646(35) s. With a precision of 0.04%, this result is an essential ingredient needed to complete a third pair of mirror superallowed transitions, 34 Ar → 34 Cl and 34 Cl → 34 S, the ratio of whose ft values would provide the most sensitive test yet of the isospin-symmetry-breaking corrections used to obtain V ud . In a sequence repeated thousands of times, we implanted a very pure 34 Ar beam into tape for 0.7 s, and then rapidly transported the collected source to a shielded 4π proportional gas counter, which detected the decay positrons from the combined decays of 34 Ar and its daughter 34 Cl. The data were analyzed as a linked parent-daughter pair. Our new result replaces one we published 13 years ago.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Calibration of the liquid argon ionization response to low energy electronic and nuclear recoils with DarkSide-50

DarkSide-50 has demonstrated the high potential of dual-phase liquid argon time projection chambers in exploring interactions of WIMPs in the GeV=c 2 mass range. The technique, based on the detection of the ionization signal amplified via electroluminescence in the gas phase, allows us to explore recoil energies down to the sub-keV range. We report here on the DarkSide-50 measurement of the ionization yield of electronic recoils down to ~180 eV er , exploiting 37 Ar and 39 Ar decays, and extrapolated to a few ionization electrons with the Thomas-Imel box model. Further, we present a model-dependent determination of the ionization response to nuclear recoils down to ~500 eV nr , the lowest ever achieved in liquid argon, using in situ neutron calibration sources and external datasets from neutron beam experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Exploring 𝛽 decay and 𝛽-delayed neutron emission in exotic 46,47 Cl isotopes

In this paper, 𝛽 − and 𝛽-delayed neutron decays of 46,47 Cl are reported from an experiment carried out at the National Superconducting Cyclotron Laboratory using the Beta Counting System. The half-lives of both 46 Cl and 47 Cl were extracted. Based on the delayed 𝛾-ray transitions observed, the level structure of 𝑁=28 46 Ar was determined. Completely different sets of excited states above the first 2 + state in 46 Ar were populated in the 46 Cl 𝛽⁢0⁢𝑛 and 47 Cl 𝛽⁢1⁢𝑛 decay channels. Two new 𝛾-ray transitions in 47 Ar were identified from the very weak 47 Cl 𝛽⁢0⁢𝑛 decay. Furthermore, 46 Cl 𝛽⁢1⁢𝑛 and 47 Cl 𝛽⁢2⁢𝑛 were also observed to yield different population patterns for levels in 45 Ar, including states of different parities. Here, the experimental results allow us to address some of the open questions related to the delayed neutron emission process. For isotopes with large neutron excess and high 𝑄 𝛽 values, delayed neutron emission remains an important decay mode and can be utilized as a powerful spectroscopic tool. Experimental results were compared with shell-model calculations using the FSU and 𝑉 MU effective interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Capacitive discharge in flash experiments: A discussion of the charged species

Abstract We report on the discharge of the capacitance formed at the electrodes in flash experiments with yttria stabilized zirconia. The experiments were carried out by disconnecting the current, and, instead, short circuiting the electrodes through a resistor. The time dependent voltage across the resistor was measured; the ratio yielded the discharge current. The current decayed exponentially with time, as expected in an RC circuit, which allowed the measurement of the capacitance. Experiments were carried out in two ways. In one case the specimens were flashed within a glove box filled with Ar (< 1 ppm O 2 ): these yielded electronic conductors. In the other case the specimens were flashed in ambient air: while electronically conducting in‐flash these specimens recover their prior insulating behavior as soon as the current is turned off. The stored charge was two orders of magnitude greater in the Ar experiments. Rather unexpectedly, the sign of the voltage expressed at the electrodes was opposite in experiments carried out in air and in Ar. The capacitance measured in the discharge experiments is attributed to the formation of space charge adjacent to the electrodes. In the case of Ar experiments, the capacitance is very large, approaching 1 F; in this case the space charge is expected to be constituted from ions. In the air experiments the specimen becomes insulating, trapping the electrons as a space charge. Hall effect measurements of the carrier density and X‐ray photoelectron spectroscopy characterization of electronically conducting single crystal specimen of cubic zirconia are reported.

Materials Science↗
DOE OSTI2021

Economic, Exergy, and Environmental Analyses of the Energy Assessments for U.S. Industries

Abstract This paper highlights the expected versus actual outcomes of 152 energy assessments that were performed between 2011 and 2020. The 1317 energy-assessment recommendations (ARs) are grouped into eight categories. This study adopted four measures per each category of recommendations: annual electricity savings, annual gas savings, annual cost savings, and annual CO2 emission reduction. The first part of the analysis compares the expected to the actually implemented values of the measures applied to each recommendation’s category. It was found that the percentages of the actual to the expected electricity, gas, and cost savings are 26.6%, 11.4%, and 17.1%, respectively, while the percentage of the actual to the expected CO2 reduction is 22%. Moreover, the second part of the analysis presents each category's implementation rate and the reasons for rejecting the unimplemented ARs. Cash flow and expensive initial investment resulted in rejecting 25% of ARs. Furthermore, the study proposes techniques and strategies to increase ARs’ implementation rate and improve private energy services companies’ implementation rate. Finally, exergy analysis is added to show the improvement that energy assessment achieves regarding exergy and exergy efficiencies of different industrial applications.

Energy & Fuels↗
DOE OSTI2020

Etching of Si 3 N 4 induced by electron beam plasma from hollow cathode plasma in a downstream reactive environment

An etching system based on the interaction of electrons extracted from a direct current (DC) hollow-cathode (HC) Ar plasma and injected towards a Si3N4 covered silicon substrate located in the downstream reactive environment created by an additional remote CF 4 /O 2 plasma source was developed and evaluated. By controlling the properties of the injected beam electrons, this approach allows to deliver energy to a surface functionalized by exposure to reactive species, and initiate surface etching. The energy of the primary beam electrons is controlled by the acceleration voltage relative to the HC discharge. Ar atoms flow from the high-pressure HC discharge into the low pressure downstream reactive environment in the process chamber. For an acceleration voltage greater than the ionization potential of Ar and/or process gas species, the energetic primary beam electrons produce a secondary plasma in the process chamber, and can also cause additional dissociation. We have characterized the properties of the secondary plasma and also surface etching of Si 3 N 4 as a function of process parameters, including acceleration voltage (0-80 V), discharge current of the HC discharge (1-2 A), pressure (3.5-20 mTorr), source to substrate distance (1.5 to 5 cm), and feed gas composition (20% and 80% O2 in CF 4 /O 2 ). The electron energy probability function (EEPF) measured with a Langmuir probe about 2.5 cm below the extraction ring suggests several major groups of electrons for this situation, including high energy primary beam electrons with an energy that varies as the acceleration voltage is changed and low energy electrons produced by beam electron-induced ionization of the Ar gas in the process chamber. When a remote CF 4 /O 2 plasma is additionally coupled to the process chamber, Si 3 N 4 surfaces can be functionalized, and by varying the energy of the beam electrons, Si 3 N 4 etching can be induced by electron-neutral synergy effect with plasma-surface interaction. For conditions without beam electron injection, the remote plasma etching rate (ER) of Si 3 N 4 depends strongly on the O 2 concentration in the CF 4 /O 2 processing gas mixture and can be suppressed for O 2 -rich process conditions by the formation of a SiONF passivation layer on the Si 3 N 4 surface. The combination of the HC electron beam source with the remote plasma source (HCEB-RP) makes it possible to induce Si 3 N 4 etching for O 2 -rich remote plasma conditions where remote plasma by itself produces negligible Si 3 N 4 etching. The electron enhanced etching of Si 3 N 4 depends strongly on O 2 /CF 4 mixing ratio reflecting changing arrival rates of O and F species at the surface. Optical emission spectroscopy was used to estimate the ratio of gas phase F and O densities and found to be controlled by the gas mixing ratio and independent of HC EB operating conditions. At this time the detailed sequence of events operative in the etching mechanism is unclear. While the increase of the electron energy is ultimately responsible for initiating surface etching, at present we cannot rule out a role of ions from the simultaneously produced secondary plasma in plasma surface interaction mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗