SEARCH · Engineering Papers

Results for “Weak form”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A BACUS S UMMIT N -body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Comparative Studies of Optoelectronic Properties, Structures, and Surface Morphologies for Phosphorus-Doped Poly-Si/SiOx Passivating Contacts

We investigated and compared optoelectronic properties, crystallographic structures, and nanoscale surface morphologies of ex-situ phosphorus-doped polycrystalline silicon (poly-Si)/SiO x passivating contacts, formed by different deposition methods (sputtering, plasma-enhanced chemical vapour deposition (PECVD), and low-pressure chemical vapour deposition (LPCVD)). Across all these deposition technologies, a similar trend is observed: higher diffusion temperatures yield films that are more crystalline but have rougher surface morphologies due to bigger surface crystal grains. Also, the recrystallization process of the as-deposited Si films starts from the SiO x interface, rather than from the film surface and bulk. However, there are some distinct differences among these technologies. Firstly, the LPCVD method yields the roughest surface and smallest degree of crystallinity on finished poly-Si films. In contrast, the PECVD method has the smoothest surface for both as-deposited Si and annealed poly-Si films. Secondly, as-deposited sputtered and PECVD Si films contain only an amorphous phase whereas as-deposited LPCVD films has already had some crystalline phase. Thirdly, the LPCVD phosphorus in-diffusion into the substrate depends strongly on the initial film thickness, whereas for the other two methods it is weakly dependent on thickness.

crystallographic structures↗

DOE OSTI2025

Measurement of the parity-violating asymmetry in the 𝑁 → Δ transition at low 𝑄 2

Here, we report the measurement of the parity-violating asymmetry in the 𝑁 → Δ transition via the 𝑒 − + 𝑝 → 𝑒 − + Δ + reaction at two different kinematic points with low four-momentum transfer 𝑄 2 . Measurements were made with incident electron beam energies of 0.877 and 1.16 GeV, corresponding to 𝑄 2 values of 0.0111 and 0.0208 (GeV/c) 2 , respectively. These measurements put constraints on a low-energy constant in the weak Lagrangian, 𝑑 Δ , corresponding to a parity-violating electric-dipole transition matrix element. This matrix element has been shown to be large in the strangeness-changing channel, via weak hyperon decays such as Σ + → 𝑝⁢𝛾. The measurements reported here constrain 𝑑 Δ in the strangeness-conserving channel. The final asymmetries were −0.65±1.00⁢(stat.)±1.02⁢(syst.) ppm (parts per million) for 0.877 GeV and −3.59±0.82⁢(stat.)±1.33⁢(syst.) ppm for 1.16 GeV. With these results we deduce a small value for 𝑑 Δ , consistent with zero, in the strangeness-conserving channel, in contrast to the large value for 𝑑 Δ previously reported in the strangeness-changing channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Qualification and Quantification of Porosity at the Top of the Fuel Pins in Metallic Fuels Using Image Processing

Approximately 130,000 metal fuel pins were irradiated in the Experimental Breeder Reactor II (EBR-II) during its 30 years of operation to develop and characterize existing and prospective fuels. For many of the metal fuel irradiation experiments, neutron radiography imaging was performed to characterize fuel behavior, such as fuel axial expansion. While several fuel expansion results obtained from neutron radiography imaging have been published, the analysis of neutron radiography for the purpose of describing statistical properties of porous matter formed on top of the fuel pins, also referred to as fluff in previous publications, is significantly less represented in the literature with just a single paper so far. This study aims to validate and augment results reported in previous publications using automated image processing. The paper describes the statistical properties of the porous matter in terms of nine parameters derived from radiography images and correlates those parameters with such fuel properties as composition, expansion, temperature, and burnup. The reported results are based on 1097 fuel pins of eight different fuel compositions. For three major fuel types, U-10Zr, U-8Pu-10Zr, and U-19Pu-10Zr, a clear negative correlation is found between the Pu content and five parameters describing the amount of porous matter generated. The parameters describing granularity properties, however, showed either negative correlation or nonlinear dependency from fuel composition. The parameters describing the amount showed a positive correlation with fuel axial expansion, while granularity parameters showed a negative correlation with axial expansion. The dependency on cladding temperature was found to be weak. A positive correlation is demonstrated for volume parameters and fuel burnup. In general, reported results confirm and validate findings published in previous studies using a much larger number of pins and automated processing techniques, which easily lend themselves to reproducibility, thus avoiding subjective bias.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Remote Epitaxy of III-V Solar Cells via Hydride Vapor Phase Epitaxy

The promise of remote epitaxy, in which epitaxial registry of a film can still be maintained through two-dimensional layers, has emerged from an advanced understanding of the interactions between two dimensional materials and bulk semiconductors. A promising application of this technology is to facilitate substrate reuse in epitaxial III-V systems, which could significantly reduce material cost and improve commercial viability of this high-efficiency photovoltaic technology. In this process the weak van der Waal bonds between the 2D layer and the III-V substrate make separation of epitaxially-grown layers simple and leaves a smooth substrate that can be reused for cost reduction. Here we present a study of growth parameters of GaAs on AlGaAs through a low dimensional carbon interlayer via hydride vapor phase epitaxy (HVPE). We Investigate the effect of various GaAs nucleation conditions such as V/III ratio, temperature and growth rate on film roughness, crystallinity, and degree of epitaxial alignment. We demonstrate growth of single crystal GaAs layer on top of the low dimensional carbon that follow the crystal orientation of the substrate. Surface roughness measured via AFM found that lower V/III ratios result in smoother III-V films grown on the low dimensional carbon, with a lowest RMS of ~60 nm obtained at a V/III ratio of 5. Growth temperatures are optimized at 650 C, as temperatures below 550 C show non-crystalline growth and higher temperature samples display very rough surfaces. Additionally, we investigate HVPE growth on various forms of low-dimensional carbon including transferred graphene, CVD-grown graphene, and MOCVD-growth amorphous carbon and find the latter to be most well suited for III-V growth as this method results in cleaner interface and potential for direct, large-area coverage. We also identify techniques to mitigate surface degradation in carbon layers as a function of both exposure time and growth environment by analyzing surface changes via x-ray reflectivity. We show growth of heteroepitaxial III-V solar cells on carbon interlayers and preliminary cell results, which are currently limited to 4% efficiency under one sun.

AlGaAs↗

DOE OSTI2021

Lyman continuum leakage from low-mass galaxies with M ⋆ < 108 M⊙

ABSTRACT We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of nine low-mass star-forming galaxies at redshifts, z, in the range 0.3179–0.4524, with stellar masses $M_\star \, \lt $ 108 M⊙ and very high specific star-formation rates sSFR ∼150–630 Gyr−1, aiming to study the dependence of leaking Lyman continuum (LyC) emission on stellar mass and some other characteristics of the galaxy. We detect LyC emission in four out of nine galaxies with escape fractions, fesc(LyC), in the range of 11–35 per cent, and establish upper limits for fesc(LyC) in the remaining five galaxies. We observe a narrow Ly α emission line with two peaks in seven galaxies and likely more complex Ly α profiles in the two remaining galaxies. The velocity separation between the peaks Vsep varies in the range from ∼229 to ∼512 km s−1. Our additional data on low-mass galaxies confirm and strengthen the tight anticorrelation between fesc(LyC) and Vsep found for previous low-redshift galaxy samples with higher stellar masses. Vsep remains the best indirect indicator of LyC leakage. It is better than O32 on which fesc(LyC) depends weakly, with a large scatter. Finally, contrary to expectations, we find no increase of fesc(LyC) with decreasing galaxy stellar mass M⋆.

Izotov, Y. I.↗

DOE OSTI2023

Unifying principles for catalytic hydrotreating processes (Final Technical Report)

This project builds on the hypothesis that the hydrotreating processes for the removal of oxygen and sulfur are fundamentally similar at the atomic-scale and existing knowledge from the treatment of petroleum derived feedstock can be leveraged for the design of novel catalysts for the upgrade of bio-oil. We tested this hypothesis by comparing computed potential energy diagrams for hydrodesulfurization (HDS) of thiophene over MoS 2 with hydrodeoxygenation (HDO) of furan over MoO 3 and concluded that certain aspects, such as catalyst promotion with transition metals, are valid strategies for both reactions. On the other hand, we also noticed significant differences in the mechanism for hydrogen (H 2 ) activation, which requires sites with metallic character. While MoS 2 is known to have metallic edge states that can catalyze H 2 dissociation, this elementary step is prohibitively slow on defect-free oxides. Only in the presence of vacancies or by creating metal/oxide interfaces can efficient H 2 activation sites during HDO be formed. The need for bifunctional catalyst when it comes to efficient and selective HDO or dehydrogenation reactions was further corroborated in joint experimental and theoretical studies of the Guerbet reaction for the coupling of biomass derived oxygenates over PdCu alloys, nitrate reduction over In-promoted Pd nanoparticles, and ethylene dehydroaromatization over Ga-exchanged ZSM-5 zeolites. All of these catalytic systems have in common that catalytic sites with distinct functional requirements are needed to create a working catalyst. Detailed computational studies were carried out for HDO of m -cresol and phenol on Ru-modified TiO 2 surfaces, which allowed us to attribute catalytic activity to the metal/oxide interface. A surprising finding was that heterolytic cleavage of the H-H bond across the Ru/TiO 2 interface was critically important, despite lower barriers for homolytic H 2 activation on Ru metal. The explanation lies in the high barriers for hydrogen spillover from Ru to TiO 2 , which becomes unnecessary in the heterolytic activation pathway. Moreover, we also reported that proton transfer steps between metal and oxide sites are mediated by weakly adsorbed surface water. During attempts to develop and validate a kinetic Monte Carlo (kMC) model for HDO reactions at the Ru/TiO 2 interface, it became clear that lateral interactions are paramount to describe realistic surface chemistry and without these interactions, the reduction and hydroxylation behavior from our simulations was inconsistent with reported experiments. To assess the importance of lateral interactions in popular computational catalyst design strategies relying on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends, we compared predicted trends with those obtained from descriptor-based kMC models. We critically evaluated the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. The mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, we concluded that the MKM approach does provide an overall better entry point for computational catalyst design. Overall, this project has led to 11 peer-reviewed publications, and their scientific is impact is well illustrated by their combined 702 citations.

08 HYDROGEN↗

DOE OSTI2021

Morphology, microstructure, and doping behaviour: A comparison between different deposition methods for poly-Si/SiO x passivating contacts

In this work, we study how crystallographic structures, optoelectronic properties, and nanoscale surface morphologies of ex situ phosphorus-doped polycrystalline silicon (poly-Si)/SiO x passivating contacts, formed by different deposition methods (sputtering, plasma-enhanced chemical vapour deposition [PECVD], and low-pressure chemical vapour deposition [LPCVD]), are investigated and compared. Across all these deposition technologies, we noted the same trend: higher diffusion temperatures yield films that are more crystalline but that have rougher surface morphologies due to bigger surface crystal grains. Also, the recrystallization process of the as-deposited Si films starts from the SiO x interface, rather than from the film surface and bulk. However, there are some distinct differences among these technologies. First, the LPCVD method yields the lowest deposition rate, roughest surfaces, and smallest degree of crystallinity on finished poly-Si films. In contrast, the PECVD method has the highest deposition rate and smoothest surfaces for both as-deposited Si and annealed poly-Si films. Second, as-deposited sputtered and PECVD Si films contain only an amorphous phase, whereas as-deposited LPCVD films already has some crystalline phase. Third, the LPCVD phosphorus in-diffusion into the substrate depends strongly on the initial film thickness, whereas for the other two methods, it is weakly dependent on thickness. Finally, the passivation quality of every poly-Si film type has different responses to the film thickness and diffusion temperature, suggesting that the ex situ doping optimization should be performed independently.

14 SOLAR ENERGY↗

DOE OSTI2020

The landscape of disc outflows from black hole–neutron star mergers

ABSTRACT We investigate mass ejection from accretion discs formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH–NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disc configurations expected from BH–NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion discs that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower limit to the fraction of the initial disc mass ejected. We find a nearly linear inverse dependence of this fraction on disc compactness (BH mass over initial disc radius). The dependence is related to the fraction of the disc mass accreted before the ouflow is launched, which depends on the disc position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disc mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disc masses. We estimate the radioactive luminosity from the disc outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH–NS parameter space, the disc outflow contribution is well below the kilonova mass upper limits for GW190814.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Elongated galactan side chains mediate cellulose–pectin interactions in engineered Arabidopsis secondary cell walls

SUMMARY The plant secondary cell wall is a thickened matrix of polysaccharides and lignin deposited at the cessation of growth in some cells. It forms the majority of carbon in lignocellulosic biomass, and it is an abundant and renewable source for forage, fiber, materials, fuels, and bioproducts. The complex structure and arrangement of the cell wall polymers mean that the carbon is difficult to access in an economical and sustainable way. One solution is to alter the cell wall polymer structure so that it is more suited to downstream processing. However, it remains difficult to predict what the effects of this engineering will be on the assembly, architecture, and properties of the cell wall. Here, we make use of Arabidopsis plants expressing a suite of genes to increase pectic galactan chain length in the secondary cell wall. Using multi‐dimensional solid‐state nuclear magnetic resonance, we show that increasing galactan chain length enhances pectin–cellulose spatial contacts and increases cellulose crystallinity. We also found that the increased galactan content leads to fewer spatial contacts of cellulose with xyloglucan and the backbone of pectin. Hence, we propose that the elongated galactan side chains compete with xyloglucan and the pectic backbone for cellulose interactions. Due to the galactan topology, this may result in comparatively weak interactions and disrupt the cell wall architecture. Therefore, introduction of this strategy into trees or other bioenergy crops would benefit from cell‐specific expression strategies to avoid negative effects on plant growth.

59 BASIC BIOLOGICAL SCIENCES↗

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 OSTI2024

Surface dependence of electronic growth of Cu(111) on MoS 2

Scanning tunneling microscopy shows that copper deposited at room temperature onto a freshly exfoliated MoS 2 surface forms Cu(111) clusters with periodic preferred heights of 5, 8, and 11 atomic layers. These height intervals correlate with Fermi nesting regions along the necks of the bulk Cu Fermi surface, indicating a connection between physical and electronic structures. Density functional theory calculations of freestanding Cu(111) films support this as well, predicting a lower density of states at the Fermi level for these preferred heights. This is consistent with other noble metals deposited on MoS 2 that exhibit electronic growth, in which the metal films self-assemble as nanostructures minimizing quantum electronic energies. Here, we have discovered that it is critical for the metal deposition to begin on a clean MoS 2 surface. If copper is deposited onto an already Cu coated surface, even if the original film displays electronic growth, the resulting Cu film lacks quantization. Instead, the preferred heights of the Cu clusters simply increase linearly with the amount of Cu deposited upon the surface. We believe this is due to different bonding conditions during the initial stages of growth. Newly deposited copper would bond strongly to the already present copper clusters, rather than the weak bonding, which exists to the van der Waals terminated surface of MoS 2 . The stronger bonding with previously deposited clusters hinders additional Cu atoms from reaching their lowest quantum energy state. The interface characteristics of the van der Waals surface enable surface engineering of self-assembled structures to achieve different applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2 ) n , n = 1,2

Pyrimidine has two in-plane CH(δ+)/N̈(δ–)/CH(δ+) binding sites that are complementary to the (δ–/2δ+/δ–) quadrupole moment of CO 2 . For this study, we recorded broadband microwave spectra over the 7.5–17.5 GHz range for pyrimidine-(CO 2 ) n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14 N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO 2 complex is planar, with CO 2 occupying one of the quadrupolar binding sites, forming a structure in which the CO 2 is stabilized in the plane by interactions with the C–H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO 2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191–195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C 2v symmetry in which the second CO 2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO 2 is –13.7 kJ mol –1 , including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO 2 ) 2 is almost exactly double that size, indicating little interaction between the two CO 2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.

abundances↗

DOE OSTI2021

Hydrogenation and C-S bond activation pathways in thiophene and tetrahydrothiophene reactions on sulfur-passivated surfaces of Ru, Pt, and Re nanoparticles

Thiophene-H2 reactions proceed via sulfur removal and hydrogenation routes on dispersed metal nanoparticles that become decorated by refractory S-adlayers during catalysis. The identity and kinetic relevance of the required elementary steps are described here based on rates measured at S-chemical potentials set by H2S/H2 ratios similar to those prevalent during practical catalysis on Re, ReSx, Ru, and Pt catalysts. Free energies for S adatom formation (from H2S decomposition and H2 evolution) are strongly exothermic (< -50 kJ mol-1 on Pt(111) and < -150 kJ mol-1 on Re and Ru(0001)), but strong repulsions between S adatoms cause adsorption free energies to increase significantly with coverage on all three surfaces, preventing complete monolayer formation. These adlayers, composed of unreactive S-atoms (S') that cover 1/3–2/3 ML leave residual interstitial spaces (*) that bind S-atoms (S*), intermediates, and transition states reversibly, as required for catalytic turnovers. The number and binding properties of these interstices depend on the identity and chemical state of the nanoparticle bulk phase, which influences S'-binding and coverages and cause large differences in direct desulfurization and hydrogenation turnover rates (per exposed metal atom) on dispersed Re, ReSx, Ru, and Pt. The identity and kinetic relevance of elementary steps for desulfurization (to C4¬ hydrocarbons) and hydrogenation (to tetrahydrothiophene; THT) are similar among these catalysts; they involve the kinetically-relevant formation of a thiophene-derived intermediate (monohydrothiophene on Re and ReSx; dihydrothiophene on Ru and Pt) that either cleaves its C-S bond or “over-hydrogenates” to THT in one surface sojourn. THT then undergoes C-S bond cleavage in secondary reactions that correct such over-hydrogenation to form the more unsaturated species that cleave C-S bonds. THT/C4 product ratios are insensitive to H2S/H2 ratios and thiophene pressure, even though active interstitial spaces are covered by kinetically-detectable coverages of S* and thiophene; therefore, primary and secondary reactions must involve the same active surfaces. The observed increase in THT/C4 ratios with H2 pressure shows that THT formation transition states involve a larger number of H-atoms than for C-S cleavage. The requirement for bound species with intermediate unsaturation (between THT and thiophene) for C-S bond cleavage is reminiscent of the H-shuttling required in C-C and C-O hydrogenolysis, reactions that involve the partial dehydrogenation of alkanes and alkanols, respectively, to weaken such bonds and to increase the formation entropy of the relevant transition states via the evolution of H2(g). These mechanistic details challenge prevalent paradigms about different site requirements for hydrogenation and desulfurization pathways and about how metal-sulfur bond energies act as descriptors of reactivity; in fact, such binding energies merely act to define the refractory S-adlayers that enable the formation of weakly-binding interstices that reversibly bind intermediates and transition states, thus allowing catalytic turnovers.

Yik, Edwin↗

DOE OSTI2024

Competing Intermolecular and Molecule–Surface Interactions: Dipole–Dipole-Driven Patterns in Mixed Carborane Self-Assembled Monolayers

Carboranedithiol isomers adsorbing with opposite orientations of their dipoles on surfaces are self-assembled together to form mixed monolayers where both lateral dipole–dipole and lateral thiol–thiolate (S–H···S) interactions provide enhanced stability over single-component monolayers. Here, we demonstrate the first instance of the ability to map individual isomers in a mixed monolayer using the model system carboranedithiols on Au{111}. The addition of methyl groups to one isomer provides both an enhanced dipole moment and extra apparent height for differentiation via scanning tunneling microscopy (STM). Associated computational investigations rationalize favorable interactions of mixed pairs and the associated stability changes that arise from these interactions. Both STM images and Monte Carlo simulations yield similarly structured mixed monolayers, where approximately 10% of the molecules have reversed dipole moment orientations but no direct chemical attachment to the surface, leading to homogeneous monolayers with no apparent phase separation. Deprotonating the thiols by depositing the molecules under basic conditions eliminates the lateral S–H···S interactions while accentuating the dipole–dipole forces. The molecular system investigated is composed of isomeric molecules with opposite orientations of dipoles and identical surface packing, which enables the mapping of individual molecules within the mixed monolayers and enables analyses of the contributions of the relatively weak lateral interactions to the overall stability of the assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Bulk viscosity from Urca processes: n p e μ matter in the neutrino-transparent regime

We study the bulk viscosity of moderately hot and dense, neutrino-transparent relativistic npeμ matter arising from weak-interaction direct Urca processes. This work parallels our recent study of the bulk viscosity of npeμ matter with a trapped neutrino component. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations—DDME2 (which does not allow for the low-temperature direct-Urca process at any density) and NL3 (which allows for low-temperature direct-Urca process above a low-density threshold). Here, we compute the equilibration rates of Urca processes of neutron decay and lepton capture, as well as the rate of the muon decay, and find that the muon decay process is subdominant to the Urca processes at temperatures T ≥ 3 MeV in the case of DDME2 model and T ≥ 1 MeV in the case of NL3 model. Thus, the Urca-process-driven bulk viscosity is computed with the assumption that pure leptonic reactions are frozen. As a result the electronic and muonic Urca channels contribute to the bulk viscosity independently and at certain densities the bulk viscosity of npeμ matter shows a double-peak structure as a function of temperature instead of the standard one-peak (resonant) form. In the final step, we estimate the damping timescales of density oscillations by the bulk viscosity. We find that, e.g., at a typical oscillation frequency f = 1~kHz, the damping of oscillation is most efficient at temperatures 3 ≤ T ≤ 5~MeV and densities n B ≤ 2n 0 where they can affect the evolution of the post-merger object.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗