Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sim”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,765 records · Page 98

Modeling O and N Alloying in Nb for SRF Applications

N and O-alloyed superconducting radio frequency cavities exhibit extraordinary quality factors. Developing diffusion models that describe interstitial N and O in Nb is important for optimizing alloyed cavity quality factors and accelerating gradients. N and O-alloyed Nb samples are examined with SEM AND SIMS and their diffusion profiles modeled.

Lechner, Eric↗

RADIATION DOSE SIMULATIONS FOR JEFFERSON LAB’S PERMANENT MAGNET RESILIENCY LDRD STUDY

In late 2023, Thomas Jefferson National Accelerator Facil- ity (Jefferson Lab) funded a Laboratory Directed Research and Development (LDRD) grant dedicated to investigating the impact of radiation on permanent magnet materials. This research initiative is specifically geared towards assessing materials slated for use in the CEBAF energy upgrade. The experimental approach involves strategically placing perma- nent magnet samples throughout the accelerator, exposing them to varying radiation doses. The simulation code BD- SIM is used to first validate the data and then to simulate the effects on future higher energy passes to study the degra- dation effects on the permanent magnets. In this paper we present the progress of that work.

Deitrick, K.↗

Dispersion Suppression for Wedge-Based Final Cooling at a 10 TeV Muon Collider

Achieving a luminosity of $\gtrsim 10^{34} cm^{-2} s^{-1}$ in a $10 \text{ } TeV$ Muon Collider, given the short lifetime of a muon, requires reducing the 6D emittance of the muon beam through a process known as ionization cooling. In the final stage of this cooling process, the transverse emittance must be reduced to $22 \text{ } μm$, typically by allowing longitudinal emittance growth up to downstream acceptance limits. While the current International Muon Collider Collaboration designs involve $40 \text{ } T$ solenoids to reach the transverse emittance target, such high-field solenoids come with several challenges, including mechanical stress management, quench protection, and potential limitations in relying on High Temperature Superconductor technology. Designed as an alternative to using such solenoids while simultaneously reaching target transverse emittance, the previously proposed wedge-based, reverse emittance-exchange cooling scheme requires excellent dispersion suppression. In this study, we design and simulate a dispersion suppressor channel for the wedge-based final cooling design that reduces dispersion in the target direction to a target value of $D_x \sim 0.001 \text{ } m$.

Karaaslan, I. [Chicago U.] (ORCID:0000000337849879↗

The Effect of Splashback on Weak Lensing Mass Estimates of Galaxy Clusters and Groups

The splashback radius of a dark matter halo, which corresponds to the first apocenter radius reached by infalling matter and substructures, has been detected around galaxy clusters using a multitude of observational methods, including weak lensing measurements. In this manuscript, we present how the splashback feature in the halo density profile affects galaxy cluster masses derived through weak lensing measurements if it is not accounted for. We find that the splashback radius has an increasingly large effect on group-sized halos towards $M_{200m} \sim 10^{13.5} \mathrm{M_\odot}$. Depending on the model and the radial scale used, the cluster/group masses can be biased low by more than 0.1 dex. This bias, in turn, would result in a slightly lower $\Omega_m$ value if propagated into a cluster cosmology analysis. The splashback effect with group-sized dark matter halos may become important to consider, given the increasingly stringent cosmological constraints coming from optical wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

pysrim: Automation, Analysis, and Plotting of SRIM Calculations

The behavior of ions traveling through a material is of great interest to many fields. For instance, radiation damage in nuclear materials and ion beam modification of materials are most concerned with understanding the formation and evolution of defects along the ion path (Was 2016). Some of the important properties that can be gleaned from investigating the initial damage from the ion include: numbers of interstitials and vacancies produced, energy deposited per unit length to electrons and atomic nuclei, track diameter, and implanted ion profile. These properties enable further simulations and allow computing a common unit of radiation damage dose, displacements per atom (dpa), which can be used to compare irradiation experiments with different ions (Backman et al. 2013) (Stoller et al. 2013). Additionally, the ejection of materials near the surface due to incident ions is important to sputtering and ion-beam analysis techniques, such as elastic recoil detection analysis (ERDA) and secondary ion mass spectrometry (SIMS) (Vickerman and Gilmore 2011). The interaction of ions within a material can be broken into two parts: electronic and nuclear stopping. Electronic stopping is the energy lost from the ion due to inelastic collisions with electrons along its path. Nuclear stopping is the energy lost due to elastic collisions between the ion and atomic nuclei within the material. The Stopping and Range of Ions in Matter (SRIM) code is a well known software in the radiation damage and ion-beam communities that allows the simulation, via Monte Carlo, of ions through a material by modeling the energy transfer through electronic and nuclear stopping (J. F. Ziegler, Ziegler, and Biersack 2010). SRIM was originally developed in 1985 and has had numerous updates on the electronic stopping powers since then, with the nuclear stopping well explained by the ZBL potential (J. F. Ziegler 1988). The executable SRIM is free to use for non-commercial use but the source code is not available to the community despite requests. While SRIM is a scientifically accurate code it does not get updated frequently and has many bugs from a usability standpoint.

97 MATHEMATICS AND COMPUTING↗

Oxygen isotope ratios in zircon and garnet: A record of assimilation and fractional crystallization in the Dinkey Dome peraluminous granite, Sierra Nevada, California

Abstract The 119 Ma Dinkey Dome pluton in the central Sierra Nevada Batholith is a peraluminous granite and contains magmatic garnet and zircon that are complexly zoned with respect to oxygen isotope ratios. Intracrystalline SIMS analysis tests the relative importance of magmatic differentiation processes vs. partial melting of metasedimentary rocks. Whereas δ18O values of bulk zircon concentrates are uniform across the entire pluton (7.7‰ VSMOW), zircon crystals are zoned in δ18O by up to 1.8‰, and when compared to late garnet, show evidence of changing magma chemistry during multiple interactions of the magma with wall rock during crustal transit. The evolution from an early high-δ18O magma [δ18O(WR) = 9.8‰] toward lower values is shown by high-δ18O zircon cores (7.8‰) and lower δ 18O rims (6.8‰). Garnets from the northwest side of the pluton show a final increase in δ18O with rims reaching 8.1‰. In situ REE measurements show zircon is magmatic and grew before garnets. Additionally, δ18O in garnets from the western side of the pluton are consistently higher (avg = 7.3‰) relative to the west (avg = 5.9‰). These δ18O variations in zircon and garnet record different stages of assimilation and fractional crystallization whereby an initially high-δ18O magma partially melted low-δ18O wallrock and was subsequently contaminated near the current level of emplacement by higher δ18O melts. Collectively, the comparison of δ18O zoning in garnet and zircon shows how a peraluminous pluton can be constructed from multiple batches of variably contaminated melts, especially in early stages of arc magmatism where magmas encounter significant heterogeneity of wall-rock assemblages. Collectively, peraluminous magmas in the Sierran arc are limited to small <100 km2 plutons that are intimately associated with metasedimentary wall rocks and often surrounded by later and larger metaluminous tonalite and granodiorite plutons. The general associations suggest that early-stage arc magmas sample crustal heterogeneities in small melt batches, but that with progressive invigoration of the arc, such compositions are more effectively blended with mantle melts in source regions. Thus, peraluminous magmas provide important details of the nascent Sierran arc and pre-batholithic crustal structure.

Geochemistry & Geophysics↗

Oxygen isotope evidence for input of magmatic fluids and precipitation of Au-Ag-tellurides in an otherwise ordinary adularia-sericite epithermal system in NE China

Abstract Tellurium-rich (Te) adularia-sericite epithermal Au-Ag deposits are an important current and future source of precious and critical metals. However, the source and evolution of ore-forming fluids in these deposits are masked by traditional bulk analysis of quartz oxygen isotope ratios that homogenize fine-scale textures and growth zones. To advance understanding of the source of Te and precious metals, herein, we use petrographic and cathodoluminescence (CL) images of such textures and growth zones to guide high spatial resolution secondary ion mass spectroscopy (SIMS) oxygen isotope analyses (10 µm spot) and spatially correlated fluid inclusion microthermometric measurements on successive quartz bands in contemporary Te-rich and Te-poor adularia-sericite (-quartz) epithermal Au-Ag vein deposits in northeastern China. The results show that large positive oxygen isotope shifts from –7.1 to +7.7‰ in quartz rims are followed by precipitation of Au-Ag telluride minerals in the Te-rich deposit, whereas small oxygen isotope shifts of only 4‰ (–2.2 to +1.6‰) were detected in quartz associated with Au-Ag minerals in the Te-poor deposits. Moreover, fluid-inclusion homogenization temperatures are higher in comb quartz rims (avg. 266.4 to 277.5 °C) followed by Au-Ag telluride minerals than in previous stages (~250 °C) in the Te-rich deposit. The Te-poor deposit has a consistent temperature (~245 °C) in quartz that pre- and postdates Au-Ag minerals. Together, the coupled increase in oxygen isotope ratios and homogenization temperatures followed by precipitation of Au-Ag tellurides strongly supports that inputs of magmatic fluid containing Au, Ag, and Te into barren meteoric water-dominated flow systems are critical to the formation of Te-rich adularia-sericite epithermal Au-Ag deposits. In contrast, Te-poor adularia-sericite epithermal Au-Ag deposits show little or no oxygen isotope or fluid-inclusion evidence for inputs of magmatic fluid.

Geochemistry & Geophysics↗

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗

Constructing the infrared conformal generators on the fuzzy sphere

We investigate the conformal algebra on the fuzzy sphere, and in particular the generators of translations and special conformal transformations which are emergent symmetries in the infinite IR but are broken along the RG flow. We show how to extract these generators using the energy momentum tensor, which is complicated by the fact that one does not have a priori access to the energy momentum tensor of the CFT limit but rather must construct it numerically. We discuss and quantitatively analyze the main sources of corrections to the conformal generators due to the breaking of scale-invariance at finite energy, and develop efficient methods for removing these corrections. The resulting generators have matrix elements that match CFT predictions with accuracy varying from sub-percent level for the lowest-lying states up to several percent accuracy for states with dimension \sim 5 ∼ 5 with N=16 N = 16 fermions. We show that the generators can be used to accurately identify primary operators vs descendant operators in energy ranges where the spectrum is too dense to do the identification solely based on the approximate integer spacing within conformal multiplets.

Fardelli, Giulia (ORCID:0000000217998124)↗

Thermal decay without information loss in horizonless microstate geometries

We develop a new hybrid WKB technique to compute boundary-to-boundary scalar Green functions in asymptotically-AdS backgrounds in which the scalar wave equation is separable and is explicitly solvable in the asymptotic region. We apply this technique to a family of six-dimensional \frac{1}{8} 1 8 -BPS asymptotically AdS _3\,\times\, 3 × S ^3 3 horizonless geometries that have the same charges and angular momenta as a D1-D5-P black hole with a large horizon area. At large and intermediate distances, these geometries very closely approximate the extremal-BTZ \,\times\, × S ^3 3 geometry of the black hole, but instead of having an event horizon, these geometries have a smooth highly-redshifted global-AdS _3\,\times\, 3 × S ^3 3 cap in the IR. We show that the response function of a scalar probe, in momentum space, is essentially given by the pole structure of the highly-redshifted global-AdS _3 3 modulated by the BTZ response function. In position space, this translates into a sharp exponential black-hole-like decay for times shorter than N_1 N_5 N 1 N 5 , followed by the emergence of evenly spaced “echoes from the cap,” with period \sim N_1 N_5 ∼ N 1 N 5 . Our result shows that horizonless microstate geometries can have the same thermal decay as black holes without the associated information loss.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring unsupervised top tagging using Bayesian inference

Recognizing hadronically decaying top-quark jets in a sample of jets, or even its total fraction in the sample, is an important step in many LHC searches for Standard Model and Beyond Standard Model physics as well. Although there exists outstanding top-tagger algorithms, their construction and their expected performance rely on Montecarlo simulations, which may induce potential biases. For these reasons we develop two simple unsupervised top-tagger algorithms based on performing Bayesian inference on a mixture model. In one of them we use as the observed variable a new geometrically-based observable \tilde{A}_{3} A ̃ 3 , and in the other we consider the more traditional \tau_{3}/\tau_{2} τ 3 / τ 2 N N -subjettiness ratio, which yields a better performance. As expected, we find that the unsupervised tagger performance is below existing supervised taggers, reaching expected Area Under Curve AUC \sim 0.80-0.81 ∼ 0.80 − 0.81 and accuracies of about 69% - − 75% in a full range of sample purity. However, these performances are more robust to possible biases in the Montecarlo that their supervised counterparts. Our findings are a step towards exploring and considering simpler and unbiased taggers.

Alvarez, Ezequiel↗

The status of the galactic center gamma-ray excess

The Galactic Center Gamma-Ray Excess has a spectrum, angular distribution, and overall intensity that agree remarkably well with that expected from annihilating dark matter particles in the form of a m_X \sim m X ∼ 50 GeV thermal relic. Previous claims that these photons are clustered on small angular scales or trace the distribution of known stellar populations once appeared to favor interpretations in which this signal originates from a large population of unresolved millisecond pulsars. More recent work, however, has overturned these conclusions, finding that the observed gamma-ray excess does not contain discernible small scale power, and is distributed with approximate spherical symmetry, not tracing any known stellar populations. In light of these results, it now appears significantly more likely that the Galactic Center Gamma-Ray Excess is produced by annihilating dark matter.

Hooper, Dan↗

Testing the mean field theory of scalar field dark matter

Scalar field dark matter offers an interesting alternative to the traditional WIMP dark matter picture. Astrophysical and cosmological simulations are useful to constraining the mass of the dark matter particle in this model. This is particularly true at low mass where the wavelike nature of the dark matter particle manifests on astrophysical scales. These simulations typical use a classical field approximation. In this work, we look at extending these simulations to include quantum corrections. We look into both the ways in which large corrections impact the predictions of scalar field dark matter, and the timescales on which these corrections grow large. Corrections tend to lessen density fluctuations and increase the effect of “quantum pressure”. During collapse, these corrections grow exponentially, quantum corrections would become important in about \sim 30 ∼ 30 dynamical times. This implies that the predictions of classical field simulations may differ from those with quantum corrections for systems with short dynamical times.

Eberhardt, Andrew↗

Exploring the unknown Lambda-neutron interaction

No published \Lambda Λ n scattering data exist. A relativistic heavy-ion experiment has suggested that a \Lambda Λ nn bound state was seen. However, several theoretical analyses have cast serious doubt on the bound-state assertion. Nevertheless, there could exist a three-body \Lambda Λ nn resonance. Such a resonance could be used to constrain the \Lambda Λ n interaction. We discuss \Lambda Λ nn calculations using nn and \Lambda Λ n pairwise interactions of rank-one, separable form that fit effective range parameters of the nn system and those hypothesized for the as yet unobserved \Lambda Λ n system based upon four different \Lambda Λ N potentials. The use of rank-one separable potentials allows one to analytically continue the \Lambda Λ nn Faddeev equations onto the second complex energy plane in search of resonance poles, by examining the eigenvalue spectrum of the kernel of the Faddeev equations. Although each of the potential models predicts a \Lambda Λ nn sub-threshold resonance pole, scaling of the \Lambda Λ n interaction by as little as \sim ∼ 5% does produce a physical resonance. This suggests that one may use photo-(electro-)production of the \Lambda Λ nn system from tritium as a tool to examine the strength of the \Lambda Λ n interaction.

Gibson, Benjamin↗

A model for $F_L$ structure function at low values of $Q^2$ and $x$ - revisited

A reanalysis of the model for the longitudinal structure function F_L (x,Q^2) F L ( x , Q 2 ) at low x x and low Q^2 Q 2 was undertaken, in view of the advent of the EIC. The model includes the kinematic constraint F_L\sim Q^4 F L ∼ Q 4 as Q^2\rightarrow Q 2 → 0. It is based on the photon-gluon fusion mechanism suitably extrapolated to the region of low Q^2 Q 2 . Revised model was critically updated, extended to the EIC kinematic region, and e.g. contains new parameterisations of the parton distribution functions.

Badelek, Barbara↗

Single-Logarithmic Corrections to Small-$x$ Helicity Evolution

The small- x x quark helicity evolution equations at double-logarithmic order, with the kernel \sim\alpha_s\ln^2(1/x) ∼ α s ln 2 ( 1 / x ) , have been derived previously. In this work, we derive the single-logarithmic corrections to the equations, to order \alpha_s\ln(1/x) α s ln ( 1 / x ) of the evolution kernel. The new equations include the effects of the running coupling and the unpolarized small- x x evolution, both of which are parametrically significant at single-logarithmic order. The large- N_c N c and large- N_c\ N_f N c N f approximations to the equation are computed. (Here, N_c N c and N_f N f are the numbers of quark colors and flavors, respectively.) Their solutions will provide more precise estimates of the quark helicity distribution at small x x , contributing to the resolution of the proton spin puzzle.

Tawabutr, Yossathorn↗

ICP-MS Workshop

This is a short document that explains the materials that will be transmitted to LLNL and DNN HQ regarding the ICP-MS Workshop held at PNNL June 17-19th. The goal of the information is to pass on to LLNL information regarding the planning and preparations for the Workshop at PNNL in preparation of the SIMS workshop at LLNL.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Sterile Neutrino Searches in MINOS/MINOS+ Experiment

This dissertation presents the searches on sterile neutrinos using the data collected in MINOS+ Experiment from September 2013 to September 2014, and the full data set of MINOS Experiment collected from 2005 to 2012. Anomalies in short baseline experiments, such as LSND and MiniBooNE, showed hints of sterile neutrinos, a type of neutrino that does not interact with the Standard Model particles. In this work, two models are considered: 3+1 and large extra dimension (LED). In the 3+1 model, one sterile neutrino state is added into the standard oscillation scheme consisting of three known active neutrino states $\nu_e$, $\nu_\mu$ and $\nu_\tau$. In the LED model, sterile neutrinos arise as Kaluza-Klein (KK) states due to assumed large extra dimensions. Mixing between sterile and active neutrino states may modify the oscillation patterns observed in the MINOS detectors. Both searches yield null results. For 3+1, a combined fit of MINOS and MINOS+ data gives a stronger limit on $\theta_{24}$ in the range of $10^{-2}\text{ eV}^2 < \Delta m^2_{41} < 1\text{ eV}^2$ than previous experiments. For LED, with the complete MINOS data set, the size of extra dimensions is constrained to be smaller than $\sim 0.35\ \mu\text{m}$ at 90\% C.L. in the limit of a vanishing lightest neutrino mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗