Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH CALCULATIONS”

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

Surface Reactivity Analysis of the Crude Oil–Brine–Limestone Interface for a Comprehensive Understanding of the Low-Salinity Waterflooding Mechanism

Low-salinity waterflooding (LSWF) has proven to improve oil recovery in carbonate formations through rock wettability alteration, although the underlying mechanism remains elusive. Multivalent ionic exchange and calcite dissolution have usually been investigated using geochemical analysis in secondary coreflooding. In this work, coreflooding, in tertiary mode, coupled with a surface reactivity analysis approach was employed to investigate the interplay of wettability alteration mechanisms such as mineral dissolution, electrostatic bond attraction, and the effect of pH at in situ conditions. Improved oil recovery (IOR) in tertiary mode observed by coreflooding in Indiana limestone rocks showed an ionic strength dependence, that is, reducing brine ionic strength resulted in an increase in oil recovery. Coreflooding results showed that the seawater and low-salinity brines deprived of Mg 2+ ions resulted in the lowest IOR in tertiary mode, indicating the significance of Mg 2+ on IOR in limestone rocks. Similar results were observed through the contact angle measurement showing the limestone rock wettability state dependence on ionic strength and the effect of Mg2+ ions. Surface reactivity analysis showed an increase in solution pH, Ca 2+ and Mg 2+ ions concentration in the effluent solution from the coreflooding in tertiary mode using low salinity brines (about 40 and 20% increase in the effluent composition for Ca 2+ and Mg 2+ , respectively). These changes in solution composition were used to calculate the in situ oil–brine and rock–brine zeta potential using a validated surface complexation model, showing the changes of zeta potential as brine is injected into limestone rocks. The results show that using seawater-like brine in tertiary mode resulted in no mineral dissolution or ionic exchange. However, improved oil recovery (IOR) using such seawater-like brine was due to wettability alteration caused by reduced electrostatic bond attraction associated with Mg 2+ ions [from 2.6 × 10 –13 (mol/m 2 ) 2 for formation water salinity to 1.5 × 10 –13 (mol/m 2 ) 2 for seawater salinity]. Using low-salinity brines in tertiary mode improved oil recovery by mineral dissolution, resulting in oil desorption and an increase in solution pH. Finally, the increase in solution pH also resulted in reduced electrostatic bond attraction which lead to rock wettability alteration using low-salinity brines.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

Electronic nature of chiral charge order in the kagome superconductor CsV 3 Sb 5

Kagome superconductors with T C up to 7 K have been discovered for over 40 y. Recently, unconventional chiral charge order has been reported in kagome superconductor KV 3 Sb 5 , with an ordering temperature of one order of magnitude higher than the T C . However, the chirality of the charge order has not been reported in the cousin kagome superconductor CsV 3 Sb 5 , and the electronic nature of the chirality remains elusive. In this paper, we report the observation of electronic chiral charge order in CsV 3 Sb 5 via scanning tunneling microscopy (STM). We observe a 2 × 2 charge modulation and a 1 × 4 superlattice in both topographic data and tunneling spectroscopy. 2 × 2 charge modulation is highly anticipated as a charge order by fundamental kagome lattice models at van Hove filling, and is shown to exhibit intrinsic chirality. We find that the 1 × 4 superlattices form various small domain walls, and can be a surface effect as supported by our first-principles calculations. Crucially, we find that the amplitude of the energy gap opened by the charge order exhibits real-space modulations, and features 2 × 2 wave vectors with chirality, highlighting the electronic nature of the chiral charge order. STM study at 0.4 K reveals a superconducting energy gap with a gap size 2Δ = 0.85 meV, which estimates a moderate superconductivity coupling strength with 2Δ/k B T C = 3.9. Furthermore, when further applying a c-axis magnetic field, vortex core bound states are observed within this gap, indicative of clean-limit superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗
DOE OSTI2025

Temperature-Dependent Mechanical Properties of Ni-Based Concentrated Alloys: Insights from First-Principles Calculations

The present work focuses on predicting temperature-dependent mechanical properties of Ni-based concentrated alloys Ni 18 Cr 10 Co 10 Fe 6 M 4 (abbreviated by X 44 M 4 , with M = Al, V, Mn, Fe, Nb, Mo, and W) using density functional theory (DFT). These predictions are based on shear (plastic) and elastic deformations, utilizing the special quasirandom structure (SQS), the phonon-based quasiharmonic approach (QHA), and the quasistatic approach. The resulting properties include coefficient of thermal expansion via QHA, ideal shear strength (τ IS ), and stable and unstable stacking fault energies (γ SF and γ US ) through pure alias shear deformation, and elastic constants (c ij ), bulk modulus (B 0 ), and shear modules (G 0 ) via elastic deformation. Notably, predicting accurate γ SF is challenging due to uncertainties that can exceed the γ SF values. τ IS and γ US exhibit a strong linear relationship, enabling the accurate prediction of γ US based on the precisely determined τ IS . All mechanical properties of X 44 M 4 decrease with increasing temperature, except for some γ SF cases such as X 44 M 4 with M = V, Mn, Fe, Mo, and W. Among the X 44 M 4 alloys, X 44 Nb 4 exhibits the lowest τ IS , γ US , and G 0 values, and the highest B 0 /G 0 ratio, while X 44 Mn 4 has the lowest B 0 and B 0 /G 0 ratio. We found that volume is a crucial descriptor for understanding and modeling mechanical properties (except B0 and maybe also γ SF ) affected by alloying elements and temperature. Ni-based dilute alloys (e.g., Ni 11 M 1 and Ni 31 M 1 ) and concentrated alloys (e.g., X 44 M 4 ) show similar trends in mechanical properties influenced by alloying elements and temperature, simplifying the analysis and design of Ni-based alloys.

Elastic properties↗
DOE OSTI2021

Towards stacking fault energy engineering in FCC high entropy alloys

Stacking Fault Energy (SFE) is an intrinsic alloy property that governs much of the plastic deformation mechanisms observed in fcc alloys. While SFE has been recognized for many years as a key intrinsic mechanical property, its inference via experimental observations or prediction using, for example, computationally intensive first-principles methods is challenging. This difficulty precludes the explicit use of SFE as an alloy design parameter. In this work, we combine DFT calculations (with necessary configurational averaging), machine-learning (ML) and physics-based models to predict the SFE in the fcc CoCrFeMnNiV-Al high-entropy alloy space. The best-performing ML model is capable of accurately predicting the SFE of arbitrary compositions within this 7-element system. Finally, this efficient model along with a recently developed model to estimate intrinsic strength of fcc HEAs is used to explore the strength–SFE Pareto front, predicting new-candidate alloys with particularly interesting mechanical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultrahigh thermal conductivity in three-dimensional covalent organic frameworks

Covalent organic frameworks (COFs) are a new, growing class of crystalline polymers with exceptional characteristics such as low density, mechanical strength, and flexibility. Taking three-dimensional (3D) COF-300 derivatives as an example, we apply molecular dynamics simulations to show that ultrahigh thermal conductivity (>15 W/(mK)) is achieved at a pore size of d = 0.63 nm. This is the first observation of such an outstanding intrinsic thermal conductivity of a 3D polymer in all directions. Further temperature-dependent thermal conductivity calculations and X-ray diffraction (XRD) patterns highlight the importance of maintaining microscopic order for a large thermal conductivity. We also demonstrate that the tunability and anisotropy in COFs provide sufficient room to engineer and obtain desired thermal conductivity. These findings provide a fundamental understanding of structure-thermal conductivity relation in COFs and shine a light on high-thermal-conductivity COFs for thermal management applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Adsorption-Induced Surface Magnetism

We report the emergence of adsorption-induced magnetism from heterohelicene molecules on a nonmagnetic Cu(100) surface. Spin-polarized low-energy electron microscopy measurements reveal spin-dependent electron reflectivity for enantiopure 7,12,17-trioxa[11]helicene (TO[11]H) monolayers, indicating the formation of a spin-polarized state localized in the topmost copper layer. Control experiments on clean Cu(100) and TO[11]H on highly oriented pyrolytic graphite show no such effect, excluding artifacts and chirality-induced spin selectivity as origins. Spin-polarized density functional theory calculations with hybrid functionals attribute the magnetism to strong chemisorption, which induces hybridization between the molecular HOMO and copper s- and d-states, driving asymmetric spin-polarized charge redistribution at the interface. An extended Newns-Anderson-Grimley model incorporating on-site Coulomb repulsion in Cu d-orbitals reproduces the emergence of interfacial spin polarization above a threshold interaction strength, highlighting the key roles of hybridization parameters and Coulomb correlation. These findings reveal a mechanism for inducing magnetism at molecule-metal interfaces without inherently magnetic components, offering avenues for engineering spin-polarized states in organic-inorganic hybrid systems.

Hubbard model↗
DOE OSTI2020

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Exploring the Photogenerated Charge Transfer Mechanism in Cu 2 O@MoS 2 Heterojunction Photocatalyst Using Transient Absorption Spectroscopy

Both heterojunction and core–shell photocatalysts have demonstrated promising performance in photocatalytic CO 2 conversions to fuels. However, fundamental knowledge of heterojunctions in core–shell structures is highly desired to facilitate the design of future photocatalysts. By combining advanced experimental characterizations and density functional theory (DFT) calculations, the role of the Cu 2 O@MoS 2 heterojunction in photocatalytic CO 2 conversions to fuels was investigated. We discovered that the charge dynamics and electron transfer properties of Cu2O@MoS 2 photocatalysts are altered by the heterojunction and Cu 2 O underlayer due to the electron transfer from Cu 2 O to MoS 2 and the change in CO 2 adsorption strength on the hybrid catalyst surface. Consequently, more electrons can travel to the surrounding liquid environment to be consumed by CO 2 reduction. This study provides experimental and theoretical investigations of the fundamental mechanisms of heterojunction core–shell photocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Pressure-induced concomitant topological and metal-insulator quantum phase transitions in Ce3Pd3Bi4

Abstract The electronic property and magnetic susceptibility of Ce 3 Pd 3 Bi 4 were systemically investigated from 18 to 290 K for varying values of cell volume using dynamic mean-field theory coupled with density functional theory. By extrapolating to zero temperature, the ground state of Ce 3 Pd 3 Bi 4 at ambient pressure is found to be a correlated semimetal due to insufficient hybridization. Upon applying pressure, the hybridization strength increases and a crossover to the Kondo insulator is observed at finite temperatures. The characteristic temperature signaling the formation of Kondo singlet, as well as the characteristic temperature associated with f -electron delocalization–localization change, simultaneously vanishes around a critical volume of 0.992 V 0 , suggesting that such metal–insulator transition is possibly associated with a quantum critical point. Finally, Wilson’s loop calculations indicate that the Kondo insulating side is topologically trivial, thus a topological transition also occurs across the quantum critical point.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Solving disorder in (3D) real space: a comparative study of the three-dimensional difference pair distribution function and atomic resolution holography reconstructions

The quantitative analysis of local ordering principles in disordered crystalline systems has gained much attention over the past few years, as it is often considered crucial for optimizing material functionality. This development has been driven by significant advancements in computational and experimental methods, which have led to the establishment and widespread use of various analytical techniques. In this study, we perform model calculations to compare the effectiveness of atomic resolution holography and three-dimensional difference pair distribution function analysis (3D-ΔPDF). Using Cu 3 Au as a model system, we demonstrate an approach to derive local order parameters quantitatively and show that both techniques are well suited to quantifying chemical short-range order correlations and local bond-distance variations. By evaluating the strengths and limitations of both techniques, we advocate for their combined use to solve complex short-range order problems accurately.

3D-ΔPDF↗
DOE OSTI2020

Connecting Nuclear Structure to Stellar Astrophysics: Neutron Skin in Tin Isotopes

The first observation of a neutron star merger by the LIGO-Virgo collaboration in 2017 highlights the need to improve our fundamental understanding of the equation of state of dense, neutron rich matter. The origin of heavy elements in the r-process and the structure of neutron stars are governed by the properties of neutron rich matter, for which experimental data is limited. Further analysis of this historic event and all future neutron star mergers relies on constraining the nuclear equation of state with experimental observables. We propose a novel method for systematically studying the evolution of the neutron skin in stable tin isotopes, by measuring the low-energy nuclear dipole strength over the broadest possible range of neutron-to-proton ratios in a single element. Nuclear resonance fluorescence with 100% linearly polarized photons from the High Intensity Gamma-ray Source (HIGS) facility was used to selectively measure the E1 photoabsorption strength of 112 Sn and 124 Sn at excitation energies from 3.5 MeV up to neutron separation, where the Pygmy Dipole Resonance dominates. The dipole polarizabilities of 112 Sn and 124 Sn were measured to be 9.03 ± 0.23 fm 3 and 9.11 ± 0.24 fm 3 , respectively. These uniquely systematically consistent measurements provide highly accurate experimental data for improving microscopic nuclear models as well as calculations of astrophysical nucleosynthesis and neutron star structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Space–time structure of weak magnetohydrodynamic turbulence

The two-time energy spectrum of weak magnetohydrodynamic turbulence is found by applying a wave-turbulence closure to the cumulant hierarchy constructed from the dynamical equations. Solutions are facilitated via asymptotic expansions in terms of the small parameter $\varepsilon$ , describing the ratio of time scales corresponding to Alfvénic propagation and nonlinear interactions between counter-propagating Alfvén waves. The strength of nonlinearity at a given spatial scale is further quantified by an integration over all possible delta-correlated modes compliant in a given set of three-wave interactions that are associated with energy flux through the said scale. The wave-turbulence closure for the two-time spectrum uncovers a secularity occurring on a time scale of order $\varepsilon ^{-2}$ , and the asymptotic expansion for the spectrum is reordered in a manner comparable to the one-time case. It is shown that for the regime of stationary turbulence, the two-time energy spectrum exponentially decays on a lagged time scale $(\varepsilon ^2 \gamma _k^s)^{-1}$ in proportion to the strength of the associated three-wave interactions, characterized by nonlinear decorrelation frequency $\gamma _k^s$ . The scaling of the form $k_{\perp } v_0 \chi _0$ exhibited by this frequency is reminiscent of random sweeping by the outer scale with characteristic fluctuation velocity $v_0$ that is modified due to competition with Alfvénic propagation (characterized by $\chi _0$ ) at the said scale. A brief calculation of frequency broadening of the power spectrum due to nonlinear interactions is also presented.

Physics↗
DOE OSTI2021

Pulse width dependence of magnetic field generation using laser-powered capacitor coils

Megagauss magnetic fields were generated by a current flowing through a U-shaped coil connecting two parallel copper foils. Two kJ-class lasers at various pulse widths from 2 ns to 9.9 ns passed through holes in the front foil and were focused on the back foil with an intensity of ~1.7×10 16 W/cm 2 . The coil current and resulting magnetic fields were characterized using ultrafast proton radiography, timed at the end of the laser pulses. Additionally, the measurements show that magnetic field strength decays with increasing laser pulse width. A lumped-circuit model was developed and showed consistency with the experimental measurements, demonstrating an ion shorting effect: as the ion current neutralizes the electron current contribution to interplate voltage, the coil current peaks on a timescale close to the ion transit time t i = d/v ion . FLASH simulations of the coil current are performed, and the calculated resistance values are used to constrain ion speed as a function of hot electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Absolute magneto-optical measurement of magnetic field on a long pulse high electrical current driven solenoid

We describe a direct absolute magneto-optical approach to measuring the magnetic field driven by long pulse width (>5 ms) and high electrical current (∼kA) in a solenoidal coil. The magnetic field and electrical current are among the most important operating parameters in lightning research and electrical power distribution and are essential to understanding the properties of the radiation output in pulsed energy accelerators as well as other research platforms. However, accurately measuring these magnetic fields and electrical currents using conventional electrical diagnostics is difficult due to the strength of ionizing radiation and electromagnetic interference. Our approach uses a fiber coupled laser beam with a rare earth element material sensing crystal that is highly resistant to electromagnetic interference, allows for absolute measurement, and does not require an external calibration. Here, we focus on device theory, operating parameters, magneto-optical results from an experiment on a high current solenoid, and comparison with a calculation of the magnetic field and reference electrical current flowing downstream of the solenoid.

Owens, Israel [Sandia National Laboratories (SNL-N↗
DOE OSTI2025

Improved pressure-gradient sensor for the prediction of separation onset in RANS models

Here, we improve upon two key aspects of the Menter shear stress transport (SST) turbulence model: (1) We propose a more robust adverse pressure gradient sensor based on the strength of the pressure gradient in the direction of the local mean flow; (2) We propose two alternative eddy viscosity models to be used in the adverse pressure gradient regions identified by our sensor. Direct numerical simulations of the Boeing Gaussian bump are used to identify the terms in the baseline SST model that need correction, and a posteriori Reynolds-averaged Navier-Stokes calculations are used to calibrate coefficient values, leading to a model that is both physics driven and data informed. The two sensor-equipped models are applied to two thick airfoils representative of modern wind turbine applications, the FFA-W3-301 and the DU00-W-212. The proposed models improve the prediction of stall (onset of separation) with respect to the prediction of the baseline SST model.

17 WIND ENERGY↗
DOE OSTI2020

Influence of hydrostatic pressure on hidden order, the Kondo lattice, and magnetism in URu 2 Si 2-x P x

Within the chemical substitution series URu 2 Si 2-x P x , there is an evolution in the ground-state behavior from hidden ordered (HO) for x ≲ 0.03, to Kondo lattice behavior with no ordering (NO) for 0.03 ≲ x ≲ 0.26, to antiferromagnetism (AFM-2) for 0.26 ≲ x ≲ 0.5 [A. Gallagher et al., Nat. Commun. 7, 10712 (2016); A. Gallagher et al., J. Phys.: Condens. Matter 29, 024004 (2016)]. To better understand what factors control this behavior, temperature-dependent electrical resistivity measurements are performed for this series under applied pressures P up to 20.5 kbar. Specimens in the HO x region show similarities to the parent compound, where HO transforms into antiferromagnetism (AFM-1) at a critical pressure (P c ). P c decreases with increasing x and collapses towards P = 0 near x ≈ 0.03, suggesting that AFM-1 occurs at ambient pressure for this concentration. No pressure-induced phase transitions are observed in the NO x region and the AFM-2 state is only weakly suppressed by P. Measurements further reveal that AFM-1 and AFM-2 are distinct from each other. Calculations of the wave functions using the tight-binding Hartree-Fock approximation are performed and show (i) that the radial probability distributions for the phosphorus ions are more tightly bound than those for the silicon and (ii) that the energy difference between the orbitals decreases with increasing x. The cumulative effect of these two factors is that Si → P substitution decreases the hybridization strength, which correlates with the weakening of HO. At large x, additional effects such as electrical charge tuning also play an important role in determining the ground-state

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗