Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dissociation”

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

An Experimental and Theoretical Study of the High Temperature Reactions of the Four Butyl Radical Isomers

The high temperature gas phase chemistry of the four butyl radical isomers (n-butyl,sec-butyl, iso-butyl, andtert-butyl) was investigated in a combined experimental and theoretical study. Organic nitrites were used as convenient and clean sources of each of the butyl radical isomers. Rate coefficients for dissociation of each nitrite were obtained experimentally and are at, or close to, the high pressure limit. Low pressure experiments were performed in a diaphragmless shock tube with laser schlieren densitometry at post-shock pressures of 65, 130, and 260 Torr and post-shock temperatures of 700-1000 K. Additional experiments were conducted with iso-butyl radicals at 805 K and 8.7 bar to elucidate changes in mechanism at higher pressures. These experiments were performed in a miniature shock tube with synchrotron-based photoionization mass spectrometry. The mass spectra confirmed that scission of the O-NO bond is the primary channel by which the precursors dissociate, but they also provided evidence of a minor channel (<7.7%) through HNO loss and formation of an aldehyde. These high pressure experiments were also used to determine the disproportionation/recombination ratio for iso-butyl radicals as 0.3. Reanalysis of the lower-temperature literature and the present data yielded rate constants for the disproportionation reaction, iso-butyl + iso-butyl = iso-butene + iso-butane. A chemical kinetics model was developed for the reactions of the butyl isomers that included new paths for highly energized adducts. These adducts are formed by the addition of H, CH 3 or C 2 H 5 to the butyl radicals. Accompanying theoretical investigations show that chemically activated pathways are competitive with stabilization of the adduct by collision under the conditions of the laser schlieren experiments. These calculations also show that at 10 bar and T< 1000 K stabilization is the only important reaction, but at higher temperatures, even at 10 bar, chemically activated product channels should also be considered. Branching fractions and rate coefficients are presented for these reactions. Furthermore, this study also highlights the importance of the radical structure for determining branching ratios for disproportionation and recombination of alkyl radicals, and these were facilitated by theoretical calculations of recombination rate coefficients for the four butyl radical isomers. Overall, the results reveal previously unknown features of butyl radical chemistry under conditions that are relevant to a wide range of applications and reaction mechanisms are presented that incorporate pressure dependent rate coefficients for the key steps.

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

Charge regulation mechanism in end-tethered weak polyampholytes

Weak polyampholytes, containing oppositely charged dissociable groups, are expected to be responsive to changes in ionic conditions. Here, we determine structural and thermodynamic properties, including the charged groups’ degrees of dissociation, of end-tethered weak polyampholyte layers as a function of salt concentration, pH, and the solvent quality. For diblock weak polyampholytes grafted by their acidic blocks, we find that the acidic monomers increase their charge while the basic monomers decrease their charge with decreasing salt concentration for pH values less than the pK a value of both monomers and vice versa when the pH > pK a . This complex charge regulation occurs because the electrostatic attraction between oppositely charged blocks is stronger than the repulsion between monomers with the same charge in both good and poor solvents when the screening by salt ions is weak. This is shown by the retraction of the top block into the bottom layer. In the case of poor solvent conditions to the basic block (the top block), we find lateral segregation of basic monomers into micelles, forming a two-dimensional hexagonal pattern on the surface at intermediate and high pH values for monovalent salt concentrations from 0.01 to 0.1 M. When the solvent is poor to both blocks, we find lateral segregation of the grafted acidic block into lamellae with longitudinal undulations of low and high acidic monomer density. By exploiting weak block polyampholytes, our work expands the parameter space for creating responsive surfaces stable over a wide range of pH and salt concentration.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Concentration-dependent ion correlations impact the electrochemical behavior of calcium battery electrolytes.

Ion interactions strongly determine the solvation environments of multivalent electrolytes even at concentrations below that required for practical battery-based energy storage. This statement is particularly true of electrolytes utilizing ethereal solvents due to their low dielectric constants. These solvents are among the most commonly used for multivalent batteries based on reactive metals (Mg, Ca) due to their reductive stability. Recent developments in multivalent electrolyte design have produced a variety of new salts for Mg2+ and Ca2+ that test the limits of weak coordination strength and oxidative stability. Such electrolytes have great potential for enabling full-cell cycling of batteries based on these working ions. However, the ion interactions in these electrolytes exhibit significant and non-intuitive concentration relationships. In this work, we investigate a promising exemplar, calcium tetrakis(hexafluoroisopropoxy)borate (Ca(BHFIP)(2)), in the ethereal solvents 1,2-dimethoxyethane (DME) and tetrahydrofuran (THF) across a concentration range of several orders of magnitude. Surprisingly, we find that effective salt dissociation is lower at relatively dilute concentrations (e.g. 0.01 M) than at higher concentrations (e.g. 0.2 M). Combined experimental and computational dielectric and X-ray spectroscopic analyses of the changes occurring in the Ca2+ solvation environment across these concentration regimes reveals a progressive transition from well-defined solvent-separated ion pairs to de-correlated free ions. This transition in ion correlation results in improvements in both conductivity and calcium cycling stability with increased salt concentration. Comparison with previous findings involving more strongly associating salts highlights the generality of this phenomenon, leading to important insight into controlling ion interactions in ether-based multivalent battery electrolytes.

Hahn, Nathan T.↗
DOE OSTI2021

Electrocatalytic nitrate reduction on rhodium sulfide compared to Pt and Rh in the presence of chloride

Chloride poisoning is a serious problem for the electrocatalytic reduction of aqueous nitrate (NO 3 – ) and improved electrocatalysts are needed. Here we study the electrocatalytic activity of rhodium sulfide supported on carbon (Rh x S y /C) for the reduction of nitrate and compare it against Pt/C and Rh/C in the presence of chloride. Between 0.05–0.15 V vs. RHE, Rh x S y /C has a steady-state nitrate reduction current density in 1 M H 2 SO 4 + 1 M NaNO 3 that is 1.6–5.6 times greater than Rh/C (the most active metal electrocatalyst) and 10–24 times greater than Pt/C. Current densities are decreased by 37% for Rh x S y /C, 62% for Rh/C, and 40% for Pt/C at 0.1 V vs. RHE in the presence of 1 mM chloride. The decrease in nitrate reduction activity for Pt, Rh, and Rh x S y is due to the competitive adsorption of chloride and nitrate on the surface. Density functional theory (DFT) modeling predicts that chloride poisoning will persistently inhibit nitrate reduction on metals due to linear adsorbate scaling relations between nitrate and chloride. DFT calculations and microkinetic modeling of our experimental measurements predict that nitrate converts to nitrite via an H-assisted dissociation mechanism on Pt and direct nitrate dissociation on Rh and Rh x S y . Pristine Rh x S y (i.e., Rh 3 S 4 , Rh 2 S 3 , and Rh 17 S 15 ) terraces are predicted to be inactive toward nitrate reduction. In contrast, sulfur vacancies in Rh 3 S 4 terraces are predicted to be active for nitrate reduction, but also bind chloride strongly. Furthermore, sulfur-defected Rh 3 S 4 rationalize the experimentally observed high activity but moderate chloride poison-resistance of Rh x S y /C for nitrate reduction.

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

Multichannel photodissociation dynamics in CS 2 studied by ultrafast electron diffraction

In this study, the structural dynamics of photoexcited gas-phase carbon disulfide (CS 2 ) molecules are investigated using ultrafast electron diffraction. The dynamics were triggered by excitation of the optically bright 1 B 2 ( 1 Σ u + ) state by an ultraviolet femtosecond laser pulse centred at 200 nm. In accordance with previous studies, rapid vibrational motion facilitates a combination of internal conversion and intersystem crossing to lower-lying electronic states. Photodissociation via these electronic manifolds results in the production of CS fragments in the electronic ground state and dissociated singlet and triplet sulphur atoms. The structural dynamics are extracted from the experiment using a trajectory-fitting filtering approach, revealing the main characteristics of the singlet and triplet dissociation pathways. Finally, the effect of the time-resolution on the experimental signal is considered and an outlook to future experiments provided.

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

Sub 20 cm -1 computational prediction of the CH bond energy – a case of systematic error in computational thermochemistry

Here, the bond dissociation energy of methylidyne, D 0 (CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller–Peterson–Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol -1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm -1 for HEAT and 27 956 ± 15 cm -1 for FPD), along with equivalent treatments of the CH ionization energy and the CH + dissociation energy (85 829 ± 15 cm -1 and 32 946 ± 15 cm -1 , respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D 0 (CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D 0 (CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm -1 for D0(CH), 32 946.7 ± 0.6 cm -1 for D 0 (CH + ), and 85 831.0 ± 6.0 cm -1 for IE(CH).

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

Voltage-dependent excitation dynamics in UV-absorbing organic photovoltaics with efficient charge transfer exciton emission

Intermolecular charge-transfer excitons play a central role in determining the performance of organic solar cells as their voltage-dependent formation, dissociation, and recombination dynamics contribute to photocurrent generation, radiative/nonradiative voltage losses, and photovoltaic fill factor. Here, we explore the properties of brightly-emitting wide energy gap (>2 eV) charge transfer excitons by measuring the voltage-dependent photoluminescence, photocurrent, and ultrafast pump–probe transient absorption spectra of organic solar cells employing five UV-absorbing donor molecules that differ only by the length of the oligophenylene or acene group at their core. We find that organic solar cells with a strong correlation between their voltage-dependent photocurrent and charge-transfer exciton photoluminescence have low photovoltaic fill factors as they require voltage to facilitate efficient charge-transfer exciton dissociation. In contrast, solar cells that are efficient can readily generate charges without an applied field and have a separate population of tightly-bound charge-transfer excitons that are responsible for emission. Furthermore, considering that the sum of all excitation loss rates (i.e., recombination and charge extraction) must be equal to the excitation generation rate in the steady state, these voltage-dependent data allow us to solve for the voltage-dependent fate of all excitations in the solar cells and estimate upper and lower bounds for geminate and non-geminate recombination, respectively.

14 SOLAR ENERGY↗
DOE OSTI2023

Prediction of the structures and heats of formation of MO 2 , MO 3 , and M 2 O 5 for M = V, Nb, Ta, Pa

Structures for the mono-, di-, and tri-bridge isomers of M 2 O 5 as well as those for the MO 2 and MO 3 fragments for M = V, Nb, Ta, and Pa were optimized at the density functional theory (DFT) level. Single point CCSD(T) calculations extrapolated to the complete basis set (CBS) limit at the DFT geometries were used to predict the energetics. Here, the lowest energy dimer isomer was the di-bridge for M = V and Nb and the tri-bridge for M = Ta and Pa. The di-bridge isomers were predicted to be composed of MO 2 + and MO 3 - fragments, whereas the mono- and tri-bridge are two MO 2 + fragments linked by an O 2 - . The heats of formation of M 2 O 5 dimers, as well as MO 2 and MO 3 neutral and ionic species were predicted using the Feller–Peterson–Dixon (FPD) approach. The heats of formation of the MF 5 species were calculated to provide additional benchmarks. Dimerization energies to form the M 2 O 5 dimers are predicted to become more negative going down group 5 and range from -29 to -45 kcal mol -1 . The ionization energies (IEs) for VO 2 and TaO 2 are essentially the same at 8.75 eV whereas the IEs for NbO 2 and PaO 2 are 8.10 and 6.25 eV, respectively. The predicted adiabatic electron affinities (AEAs) range from 3.75 eV to 4.45 eV for the MO 3 species and vertical detachment energies from 4.21 to 4.59 eV for MO 3 - . The calculated M[double bond, length as m-dash]O bond dissociation energies increase from 143 kcal mol -1 for M = V to ~170 kcal mol -1 for M = Nb and Ta to ~200 kcal mol -1 for M = Pa. The M–O bond dissociation energies are all similar ranging from 97 to 107 kcal mol -1 . Natural bond analysis provided insights into the types of chemical bonds in terms of their ionic character. Pa 2 O 5 is predicted to behave like an actinyl species dominated by the interactions of approximately linear PaO 2 + groups.

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

A weakly ion pairing electrolyte designed for high voltage magnesium batteries

High-voltage rechargeable magnesium batteries (RMBs) are potential alternatives to lithium-ion batteries owing to the low cost and high abundance of magnesium. However, the parasitic reactions of the latter with many electrolytes greatly hinders the stability and kinetics of Mg plating/stripping. Here we report a new and easily accessible solvent-designed electrolyte, which effectively solves the difficulty of ion pair dissociation and facilitates fast nanoscale Mg nucleation/growth using simple Mg(TFSI)2 as the salt, enabling a facile interfacial charge transfer process. Dendrite-free Mg plating/stripping is maintained for over 7000 hours (∼10 months) at a practical areal capacity of 2 mA h cm−2. The high-voltage stability of these electrolytes is demonstrated by benchmarking with polyaniline||Mg full cells with an operating voltage up to 3.5 V that exhibit stable cycling at a 2C rate with 99% coulombic efficiency after 400 cycles. This work opens up new frontiers in coupling low-cost electrolytes with next-generation high-voltage cathode materials for fast-charging RMBs with long life and high energy densities.

Li, Chang↗
DOE OSTI2024

First principles density functional theory study of tritium species adsorption on Ni(111) surface and diffusion in nickel-sublayer for tritium storage

The nickel-plated zircaloy-4 is used as a tritium ( 3 H) getter in the tritium-producing burnable absorber rods (TPBARs) to capture 3 H produced in the 6 Li-riched annular γ-LiAlO 2 pellet under neutron irradiation. The experimental data and our previous theoretical results showed that the 3 H species produced from the γ-LiAlO 2 pellet were mainly 3 H 2 and 3 H 2 O. These 3 H species diffuse from the surface of the LiAlO 2 pellet across vacuum to the nickel-plated zircaloy-4 getter and then further diffuse into the getter to chemically form metal hydrides. While a number of studies show that oxygen binds strongly as compared to 3 H on the nickel (Ni) layer, the detailed mechanism of 3 H species absorption and diffusion across the Ni plate and Ni/Zr interface are still unclear. By employing density functional theory calculations, here we explored the 3 H 2 and 3 H 2 O species adsorption and dissociation on the Ni(111) surface and diffusion into the Ni sublayer. Our results indicated that the 3 H 2 and 3 H 2 O dissociate on the Ni(111) surface. The NiO x and Ni(O 3 H) x could be formed in the Ni layer due to the higher oxygen (O) diffusion energy barrier and formation of Ni vacancy defects. The oxygen was found to be retained in the Ni layer from diffusing across the Ni–Zr interface. This was revealed by comparing the diffusion barriers for 3 H with O. 3 H was found to have nearly three times smaller diffusion barrier than for O, making 3 H comparatively easier to diffuse through the Ni layer. In conclusion, the obtained results provide guidelines for experimental measurements on 3 H retention behavior in TPBARs and may open further avenues to explore the impurity effects on 3 H diffusion and storage at the Ni/zircaloy interfaces.

Tafen, De Nyago [National Energy Technology Lab. (↗
DOE OSTI2025

Detection of Br and I as atomic anions using liquid sampling – atmospheric pressure glow discharge/Orbitrap mass spectrometry

The quantification and determination of halogen isotope composition is essential in many fields including geochemistry and nuclear forensics. Detection of the halogen elements is commonly attempted with inductively coupled plasma mass spectrometry (ICP-MS) however, there are multiple challenges faced. Most significantly this includes very poor ionization efficiency and the potential for isobaric interferences from either matrix or plasma species. Thus, sensitivity and accuracy are often limiting issues. Here a liquid sampling-atmospheric pressure glow discharge (LS-APGD) ionization source is coupled with an ultrahigh resolution Orbitrap MS to perform halogen detection of bromine and iodine as initial analytes. This facilitates simple and sensitive detection of these elements as atomic anions from simple aqueous salt solutions. Collision-induced dissociation (CID) and higher-energy collisional dissociation (HCD) energies were optimized to produce the maximum response of Br − and I − . The LS-APGD conditions were optimized using a design of experiments (DOE) approach for Br − and I − concurrently. Response curves for Br − and I − solutions determined limit of detection (LOD) values of 50 pg and 5 pg, respectively, in 20 μL aliquots. The curves indicated response factors of 0.67 for 79 Br − and 0.90 for 127 I − . Br isotope ratios were determined with precision between 0.7 and 7.3% RSD, with the isotope ratios determined with precision 0.8% RSD at concentrations above the limit of quantification. Preliminary tests were conducted to evaluate the effect of different cations (Na + , K + , and Mg 2+ ) on Br − responses, with no discernible impacts observed on the halogen signal responses. This study demonstrates the potential use of the LS-APGD-Orbitrap-MS for the detection of multiple halogens while avoiding interferences, minimizing sample preparation, and overcoming ionization barriers.

Lang, Dehlia A. [Clemson University, SC (United St↗
DOE OSTI2024

Ultraviolet H 2 luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B′, D, B″, D′, and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O 2 + , S2, SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Quarkonium suppression in strongly coupled plasmas

Suppression of open heavy quarks and quarkonia in heavy-ion collisions are among the most informative probes of quark-gluon plasma (QGP). Interpreting the full wealth of data obtained from the collision events requires a precise theoretical understanding of the evolution of heavy quarks and quarkonia as they propagate through a strongly coupled plasma. Such calculations require the evaluation of a gauge-invariant correlator of chromoelectric fields. This chromoelectric correlator encodes all the characteristics of QGP that the dissociation and recombination dynamics of quarkonium are sensitive to, which is to say can in principle measure. We review its distinctive qualitative features at weak coupling in QCD up to next-to-leading order and at strong coupling in N = 4 SYM using the AdS/CFT correspondence, as well as its formulation in Euclidean QCD. Furthermore, we report on recent progress in applying our results to the calculation of the final quarkonium abundances after propagating through a cooling droplet of QGP, which illustrates how we may learn about QGP from quarkonium measurements. We devote special attention to how the presence of a strongly coupled plasma modifies the transport description of quarkonium, in comparison to approaches that rely on weak coupling approximations to describe quarkonium dissociation and recombination.

Scheihing-Hitschfeld, Bruno [University of Califor↗
DOE OSTI2020

Improvements in the orbitalwise scaling down of Perdew–Zunger self-interaction correction in many-electron regions

The Perdew-Zunger (PZ) method offers a way to remove the self-interaction (SI) error from density functional approximations on an orbital by orbital basis. The PZ method provides significant improvements for the properties such as barrier heights or dissociation energies but results in over-correcting the properties well described by SI-uncorrected semi-local functional. One cure to rectify the over-correcting tendency is to scale down the magnitude of SI-correction of each orbital in the many electron region. We have implemented the orbitalwise scaled down SI-correction (OSIC) scheme of Vydrov et al. [J. Chem. Phys. 124, 094108 (2006)] using the Fermi-L owdin SI-correction method. After validating the OSIC implementation with previously reported OSICLSDA results, we examine its performance with the most successful non-empirical SCAN meta-GGA functional. Using different forms of scaling factors to identify one-electron regions, we assess the performance of OSIC-SCAN for a wide range of properties: total energies, ionization potentials and electron affinities for atoms, atomization energies, dissociation and reaction energies, and reaction barrier heights of molecules. Our results show that OSIC-SCAN provides superior results than the previously reported OSIC-LSDA, -PBE, and -TPSS results. Furthermore, we propose selective scaling of OSIC (SOSIC) to remove its major shortcoming that destroys the -1/r asymptotic behavior of the potentials. The SOSIC method gives the highest occupied orbital eigenvalues practically identical to those in PZSIC and unlike OSIC provides bound atomic anions even with larger powers of scaling factors. Furthermore, SOSIC compared to PZSIC or OSIC provides more balanced description of total energies and barrier heights.

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

Interaction of hydrogen with actinide dioxide (011) surfaces

The corrosion and oxidation of actinide metals, leading to the formation of metal-oxide surface layers with the catalytic evolution of hydrogen, impacts the management of nuclear materials. In this work, the interaction of hydrogen with actinide dioxide (AnO 2 , An = U, Np, or Pu) (011) surfaces by Hubbard corrected density functional theory (PBEsol+U) has been studied, including spin–orbit interactions and non-collinear 3k anti-ferromagnetic behavior. The actinide dioxides crystalize in the fluorite-type structure, and although the (111) surface dominates the crystal morphology, the (011) surface energetics may lead to more significant interaction with hydrogen. The dissociative adsorption of hydrogen on the UO 2 (0.44 eV), NpO 2 (-0.47 eV), and PuO 2 (-1.71 eV) (011) surfaces has been calculated. It is found that hydrogen dissociates on the PuO 2 (011) surface; however, UO 2 (011) and NpO 2 (011) surfaces are relatively inert. Recombination of hydrogen ions is likely to occur on the UO 2 (011) and NpO 2 (011) surfaces, whereas hydroxide formation is shown to occur on the PuO 2 (011) surface, which distorts the surface structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydration structure of flat and stepped MgO surfaces

We investigate the solvation structure of flat and stepped MgO(001) in neutral liquid water using ab initio molecular dynamics based on a hybrid density functional with dispersion corrections. Our simulations show that the MgO surface is covered by a densely packed layer of mixed intact and dissociated adsorbed water molecules in a planar arrangement with strong intermolecular H-bonds. The water dissociation fractions in this layer are >20% and >30% on the flat and stepped surfaces, respectively. Slightly above the first water layer, we observe metastable OH groups perpendicular to the interface, similar to those reported in low temperature studies of water monolayers on MgO. Here, these species receive hydrogen bonds from four nearby water molecules in the first layer and have their hydrophobic H end directed toward bulk water, while their associated protons are bound to surface oxygens. The formation of these OH species is attributed to the strong basicity of the MgO surface and can be relevant for understanding various phenomena from morphology evolution and growth of (nano)crystalline MgO particles to heterogeneous catalysis.

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

Phase field modeling of dislocations and obstacles in InSb

We present a phase-field dislocation dynamics (PFDD) model informed by first-principle calculations to elucidate the competitive dislocation nucleation and propagation between the glide and shuffle sets in InSb diamond cubic crystal. The calculations are directly informed with generalized stacking fault energy curves on the (111) slip plane for both the “glide set,” with the smaller interplanar spacing, and the “shuffle set,” with the larger interplanar spacing. The formulation also includes elastic anisotropy and the gradient term associated with the dislocation core. The PFDD calculations show that under no stress the equilibrium structure of screw glide set dislocations dissociates into Shockley partials, while those of the shuffle set dislocations do not dissociate, remaining compact. The calculated dislocation core widths of these InSb dislocations agree well with the measured values for other semiconductor materials, such as Si and GaN. We find that a shuffle set dislocation emits from a dislocation source at an applied stress about three times smaller than that needed to emit leading and trailing partials successively on the glide set plane. Once the partial dislocations in the glide set are emitted, they propagate faster than the shuffle set perfect dislocation at the same stress level.

36 MATERIALS SCIENCE↗