Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective core potentials”

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 181 records · Page 10

Spectroscopic signature of negative electronic compressibility from the Ti core-level of titanium carbonitride MXene

Two-dimensional transition metal carbides, carbonitrides, and nitrides called MXenes exhibit high metallic conductivity, ion intercalation capability and reversible redox activity, prompting their applications in energy storage and conversion, electromagnetic interference (EMI) shielding, and electronics, among many other fields. It has been shown that replacement of about 50% of carbon atoms in the most popular MXene family member, titanium carbide (Ti 3 C 2 T x ), by nitrogen atoms, forming titanium carbonitride (Ti 3 CNT x ), leads to drastically different properties, such as very high negative charge in solution and extreme EMI shielding effectiveness, exceeding all known materials, even metals at comparable thicknesses. Here, by using ultraviolet photoemission spectroscopy (UPS), the electronic structures of Ti 3 CNT x and Ti 3 C 2 T x are systematically investigated and compared as a function of charge carrier density. We observe that, in contrast to Ti 3 C 2 T x , the Ti 3p core-level of Ti 3 CNT x exhibits a counterintuitive shift to a lower binding energy of up to approximately 250 meV upon increasing the electron density, which is a spectroscopic signature of negative electronic compressibility (NEC). These experimentally measured chemical potential shifts are well-captured by the density functional theory (DFT) calculation. The DFT results also further suggest that the hybridization of titanium-nitrogen bonding in Ti 3 CNT x helps promoting the available states of Ti atoms for receiving more electron above the Fermi level and leads to the observed NEC. Furthermore, our findings explain the differences in electronic properties between the two very important and widely studied MXenes and also suggest a new strategy to apply the NEC effect of Ti 3 CNT x in energy and charge storage applications.

2D materials↗

Integrated experimental studies of pore structure and fluid uptake in the Bossier Shale in eastern Texas, USA

Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.

58 GEOSCIENCES↗

A theoretical study of solid solution strengthening in the refractory medium entropy alloy Nb 45 Ta 25 Ti 15 Hf 15

The refractory medium-entropy alloy (RMEA) Nb 45 Ta 25 Ti 15 Hf 15 exhibits exceptional tensile ductility and fracture toughness at ambient temperature, but its engineering applications are limited by a lack of high temperature strength. Using a machine-learning interatomic potential (MLIP) with near-density functional theory (DFT) accuracy, we conducted molecular dynamics (MD) and statics simulations of the behavior of dislocations with both screw and edge characters. We also analyze experimentally measured yield strengths using the Rao-Suzuki model and the Maresca-Curtin model modified to include a temperature-dependent shear modulus and a bulk modulus-dependent misfit volume, thereby uncovering the mechanisms underlying the yielding of this RMEA. Compared with the published experimental yield strength, the models parameterized by the MLIP effectively reproduce the experimental results over a wide temperature range. The models and MD simulations indicate that yielding is governed by screw dislocations, with dipole dragging as the dominant mechanism. In MD simulations, we observed a potential softening mechanism not considered by the Rao-Suzuki screw model: slow migration of interstitial jogs along the dislocation core, which could lead to the annihilation of vacancy and interstitial jog pairs by their combination.

BCC complex concentrated alloys↗

Anisotropic surface potentials induced by competitive ion adsorption enable the synthesis of branched cubic Pt mesocrystals

Creation of complex nanostructured materials through oriented attachment (OA) requires the manipulation of interparticle forces, including electrostatic repulsion, which depends strongly on surface potentials and can be modified through the effect of solution environment on interfacial chemistry. Here we show that time-dependent anisotropies in surface potential driven by competitive ion adsorption can alter facet-selectivity during OA. This phenomenon enables the synthesis of branched cubic Pt mesocrystals. Initially, Pt nanoparticles attach preferentially at their {100} facets to form a well-defined cubic core. Over time, changes in ion adsorption shift the attachment preference to the {111} facets, promoting branch formation. In both stages, anisotropic surface potentials generate electrostatic torques that align the particles prior to attachment. These findings demonstrate a generalizable strategy for directing the architecture of nanomaterials through time-resolved control of interfacial chemistry during OA, offering new pathways for the design of complex mesoscale structures.

Bae, Yuna [Pacific Northwest National Laboratory (↗

Conceptual design of inverted core lead bismuth eutectic fast reactor for marine applications

The development of an inverted core fast reactor aims to generate 60 MWth for about 30 Effective Full Power Years without refueling. The reactor design is a transportable reactor using UO{sub 2} fuel and lead-bismuth-eutectic cooled designed for marine applications and is intended to improve the reactor performances compared to the normal core design: better condition for passive cooling system capability by lower core pressure drop, taking advantage of potential power uprate from the lower maximum fuel temperature. Systematic design processes are presented in this work: fuel pin geometry selection, fuel assembly (FA) design, and core design. A relationship between pressure drops, coolant velocity, maximum fuel temperature, coolant channel diameter, and fuel volume fraction was introduced in a single graph used as a tool to select fuel pin geometry. Fuel fabrication capability also took place in consideration of FA design which led to 7 holes per FA, and two-dimensional temperature distribution studies were also carried out. Core design processes including radial zoning, axial zoning, and core optimization were conducted using Monte Carlo code MCS, which is UNIST CORE laboratory in-house code. The current core design uses 3 fuel enrichment levels and 3 FA types to control the local power distribution and power shift during its lifetime. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗

Determining the Structure-Property Relationships of Quasi-Two-Dimensional Semiconductor Nanoplatelets.

We report a theoretical study of CdSe nanoplatelets aimed at identifying the main factors determining their photophysical properties. Using atomic configurations optimized with density functional theory calculations, we computed quasiparticle and exciton binding energies of nanoplatelets with two to seven monolayers. We employed many body perturbation theory at the GW level and solved the Bethe-Salpeter equation to obtain absorption spectra and excitonic properties. Our results, which agree well with recent experiments, were then used to design a model that allows us to disentangle the effects of quantum confinement, strain induced by passivating ligands, and dielectric environment on the electronic properties of nanoplatelets. We found that, for the model to accurately reproduce our first principle results, it is critical to account for surface stress and consider a finite potential barrier and energy-dependent effective masses when describing quantum confinement. Our findings call into question previous assumptions on the validity of an infinite barrier to describe carrier confinement in nanoplatelets, suggesting that it may be possible to optimize interfacial charge transfer and extraction by appropriately choosing passivating ligands. The model developed here is generalizable to core-shell platelets and enables the description of system sizes not yet directly treatable by first-principles calculations.

Greenwood, Arin R.↗

An attractive model: simulating fuzzy dark matter with attractive self-interactions

ABSTRACT Fuzzy dark matter (FDM), comprised of ultralight ($m \sim 10^{-22}\,{\rm eV}$) boson particles, has received significant attention as a viable alternative to cold dark matter (CDM), as it approximates CDM on large scales (${\gtrsim}1$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (the boson mass), is subject to a tension: small boson masses yield the desired cores of dwarf galaxies but underpredict structure in the Lyman-α forest, while large boson masses render FDM effectively identical to CDM. This Catch-22 problem may be alleviated by considering an axion-like particle with attractive particle self-interactions. We simulate an idealized FDM halo with self-interactions parametrized by an energy decay constant $f \sim 10^{15}~\rm {GeV}$ related to the axion symmetry-breaking conjectured to solve the strong-CP problem in particle physics. We observe solitons, a hallmark of FDM, condensing within a broader halo envelope, and find that the density profile and soliton mass depend on self-interaction strength. We propose generalized formulae to extend those from previous works to include self-interactions. We also investigate a critical mass threshold predicted for strong interactions at which the soliton collapses into a compact, unresolved state. We find that the collapse happens quickly, and its effects are initially contained to the central region of the halo.

Painter, Connor A. (ORCID:0000000235314806)↗

Accelerating the Structure Exploration of Diverse Bi–Pt Nanoclusters via Physics‐Informed Machine Learning Potential and Particle Swarm Optimization

Bimetallic Bi–Pt nanoclusters exhibit diverse structural motifs, including core-shell, Janus, and mixed alloy configurations, due to the unique bonding characteristics between Bi and Pt atoms. Using density functional theory refinements from ChIMES physically machine-learned potential and CALYPSO particle swarm optimization global searches, 34 Bi20-Pt20 nanoclusters are systematically classified. The results reveal that Bi atoms predominantly occupy surface sites, driven by charge transfer effects. Cohesive energy trends alone prove insufficient for structure differentiation, necessitating a data-driven approach employing principal component analysis and K-means clustering. Furthermore, vibrational, electronic, and infrared spectral analyses provide additional insights into structure-property relationships. The findings offer an original framework for the automated classification and analysis of bimetallic nanoclusters, enhancing the understanding of their stability and functional properties.

bimetallic nanoparticles↗

Multi-GPU porting of a phase-change cascaded lattice Boltzmann method for three-dimensional pool boiling simulations

The Lattice Boltzmann method (LBM) has proven effective in simulating phase-change phenomena, such as melting, solidification, evaporation, and boiling. In this work, we develop a highly parallelized multi-GPU implementation of LBM for three-dimensional pool boiling simulations. The code is based on the OpenACC programming model, which enables the code to be deployed efficiently on multi-core CPUs, GPUs, and potentially other accelerators, without the need for architecture-specific rewrites. To support large-scale simulations, the domain is decomposed and distributed across multiple compute nodes using MPI. We demonstrate that the code exhibits excellent scaling properties, with ideal strong-scaling running with up to 256 GPUs on the MareNostrum5 cluster.

97 MATHEMATICS AND COMPUTING↗

Ironsand (Titanomagnetite-Titanohematite): Chemistry, Magnetic Properties and Direct Applications for Wireless Power Transfer

Ironsand is an abundant and inexpensive magnetic mineral resource. However, the magnetic properties of unprocessed ironsand are often inadequate for any practical applications. In this work, the applicability of ironsand for use as a component in a soft magnetic composite for large-scale inductive power transfer applications was investigated. After magnetic separation, the chemical, structural and magnetic properties of ironsand sourced from different locations were compared. Differences observed in the DC magnetic properties were consistent with changes in the chemical compositions obtained from X-ray Absorption Near-Edge Spectroscopy (XANES), which suggests varying the titanohematite to titanomagnetite content. Increased content in titanomagnetite and magnetic permeability correlated well with the total Fe content in the materials. The best-performing ironsand with the highest permeability and lowest core losses was used alongside Mn,Zn-Ferrite particles (ranging from ~100 μm to 2 mm) to fabricate toroid cores with varying magnetic material loading. It was shown that ironsand can be used to replace up to 15 wt.% of the magnetic materials with minimal impact on the composite magnetic performance, thus reducing the cost. Ironsand was also used as a supporting material in a single-rail wireless power transfer system, effectively increasing the power transfer, demonstrating potential applications to reduce flux leakage.

36 MATERIALS SCIENCE↗

Concentric core-shell tracks and spectroscopic properties of SrTiO 3 under intense electronic excitation

A deeper understanding of the intrinsic link between the irradiation-induced microstructures and the corresponding spectroscopic properties is becoming increasingly attractive for the research fields of materials science, physics, and information technology. In this work, the structural damage response, spectroscopic features, and associated physical mechanisms of SrTiO 3 single crystals under swift heavy ion irradiation were comparatively analyzed using a combination of experimental and theoretical approaches. Corresponding to 0.11–5.00 MeV/u ion irradiation with electronic energy loss ranging from 4.0 to 29.3 keV/nm, the inelastic thermal spike calculations combined with molecular dynamics simulations are compared with the experimental observations, revealing the track fine structures (individual spherical defects with a disordered region, and discontinuous and continuous tracks consisting of an amorphous core and a disordered outer shell) and demonstrating the dominant effects of the deposition energy and lattice temperature on track damage formation and evolution; thus, two essential thresholds for defect formation of ~0.60 eV/atom and amorphous region formation of ~1.81 eV/atom were identified to better describe the concentric core-shell track structure. The enhanced nanohillock formation in the surface region is attributed to the combined action of kinetic and potential energy depositions, and the fluence dependence in regulating the hillock dimensions is also presented. The measured refractive index profiles along different crystal axes further indicate the anisotropy of irradiation-induced lattice expansion. With increasing ion fluence and damage level, the intrinsic bandgap gradually decreased, and the concentration of additional radiative recombination centers (1.81 eV, 2.31 eV, and 2.42 eV) accordingly increased in the SrTiO 3 crystal, providing the possibility of regulating related defects to achieve tunable/selective photon emission for more critical technological applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Influence of living grass Roots and endophytic fungal hyphae on soil hydraulic properties

Soil hydraulic properties are often estimated based on laboratory data or pedotransfer functions dependent on soil physical properties, which often do not consider potential impacts of soil roots or fungal hyphae. Here, we first review current knowledge of how these soil biotic components affect hydraulic properties, then we conducted laboratory experiments to specifically test if the presence of roots and mycorrhizal fungi had a significant effect on the hydraulic properties of two soils with contrasting textures: Flint sand and Hamblen silt loam. Soil cores were seeded with (Panicum virgatum) and grown in a greenhouse over three separate growth periods. The endophytic fungus Serendipita indica was injected as liquid inoculant into designated mycorrhizal cores. Saturated hydraulic conductivity (K sat ) measurements were made with a constant head permeameter, and soil water retention curves were obtained by the evaporation method, supplemented at the dry end for Hamblen silt loam with water activity meter data. Retention curve parameters were obtained by fitting the van Genuchten equation to the resulting measurements. Mean root volume ratios were higher in the mycorrhizal inoculated treatment than in the uninoculated treatment for both soils. For Flint sand, analysis of variance revealed that K sat was reduced by the presence of roots as compared to bare soil. This was likely due to roots clogging soil pores. Results also indicated the presence of roots changed the shape of the water retention curve for Flint sand by increasing water content at saturation and by reducing the slope of the curve. These changes suggested roots created additional porosity and broadened the pore-size distribution. The presence of mycorrhizal fungi accentuated the root effects. The influence of roots and mycorrhizal fungi on hydraulic properties was less obvious for the Hamblen silt loam, as none of the treatments differed from each other at p < 0.05. The results highlight the necessity to consider the impact of root and fungal structures on models of soil hydraulic properties.

59 BASIC BIOLOGICAL SCIENCES↗

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of Postulated Accident Scenarios for the Transformational Challenge Reactor

The Transformational Challenge Reactor (TCR) is a helium-cooled, yttrium-hydride-moderated reactor that was designed for the U.S. Department of Energy Office of Nuclear Energy. A key objective of the TCR was to employ advanced manufacturing techniques in a nuclear system and demonstrate their potential for revolutionizing the nuclear reactor design process. One purpose of the present work is to demonstrate the safety of the TCR under postulated accidents. Based on RELAP5-3D and COMSOL analyses, the TCR remained below all current safety limits and far below the expected failure limits for the core materials. Another purpose of this work is to provide useful insights and recommendations regarding the application of RELAP5-3D to gas-cooled or other advanced reactors. A novel approach was implemented for simultaneously modeling conduction and radiation in RELAP5-3D, which was found to provide reasonable predictions of radial core, vessel, and ex-vessel heat transfer during postulated events. A multicode approach was also applied, in which high-fidelity COMSOL calculations were used to tune the radial heat transfer parameters in RELAP5-3D. The tuned RELAP5-3D model demonstrated comparable peak temperature predictions as COMSOL, despite a coarse treatment of the core in RELAP5-3D consisting of only two lumped heat structures. This high-fidelity tuning approach enabled enhanced accuracy as well as minimal complexity within the RELAP5-3D model, even for complex fuel geometric designs as in the TCR. Finally, investigations were made into the potential for flow reversal during a pressurized loss-of-forced-flow event in the TCR. The TCR is designed with downward helium flow through the core during normal operation. The RELAP5-3D model predicted that this downward flow would persist, without flow reversal, up to several days after the circulator trip. This was attributed to natural circulation hysteresis effects as have been noted in similar thermofluidic systems. Although flow stagnation and eventual reversal did not lead to unsafe TCR conditions, interesting spatial effects were observed which may have safety relevance for other reactor system designs and coolant types that are designed for downward core flow during normal operation, warranting closer investigation of the flow reversal phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CELLULAR COFFERDAMS AS PERMANENT HYDROPOWER DAM STRUCTURES

This paper presents the results of a comprehensive study on the potential use of cellular cofferdams as basis for the design and construction of water retaining structures to sustainably and cost-effectively harness hydropower. Previously, cellular cofferdams have been widely used mainly as temporary water exclusion devices to permit dry construction of in-water structures such as dams, locks, bridge footings and piers, and hydroelectric power plants. Design and construction requirements for cellular cofferdams are less stringent than for hydropower dams. To make cellular cofferdams suitable for permanent hydropower use, different design concepts that utilize cellular cofferdams as the main or core element of the water-retaining dam structure are proposed. One particular key design concept is the so-called “dry construction technique” in which the granular fill in cofferdam cells and the downstream berm are permanently kept dry in contrast to the wet construction technique for temporary use of cellular cofferdams. The viability of the proposed permanent cellular cofferdam design concepts is demonstrated using well-established structural and geotechnical design procedures and computational modeling. The improved performance of the proposed design concepts, particularly in combination with the dry construction technique, show cellular cofferdams have the potential to be used as basis for the construction of permanent hydropower dam structures that are versatile, with less impact on the environment, and will cost less to build than conventional hydropower dams.

13 HYDRO ENERGY↗

Using Substituted [Fe 4 N(CO) 12 ] – as a Platform To Probe the Effect of Cation and Lewis Acid Location on Redox Potential

The impact of cationic and Lewis acidic functional groups installed in the primary or secondary coordination sphere (PCS or SCS) of an (electro)catalyst is known to vary depending on the precise positioning of those groups. However, it is difficult to systematically probe the effect of that position. In this report we probe the effect of the functional group position and identity on the observed reduction potentials (E p,c ) using substituted iron clusters, [Fe 4 N(CO) 11 R] n , where R = NO + , PPh 2 -CH 2 CH 2–9 BBN, ( Me PTA + ) 2 , ( Me PTA + ) 4 , and H + , where n = 0, -1, +1, or +3 (9-BBN is 9-borabicyclo(3.3.1)nonane and Me PTA + is 1-methyl-1 azonia-3,5-diaza-7-phosphaadamantane). The cationic NO + and H + ligands cause anodic shifts of 700 and 320 mV, respectively, in E p,c relative to unsubstituted [Fe 4 N(CO) 12 ] – . Infra-red absorption band data, νCO, suggests that some of the 700 mV shift by NO + results from electronic changes to the cluster core. Furthermore, this contrasts with the effects of cationic MePTA + and H + which cause primarily electrostatic effects on E p,c . Lewis acidic 9-BBN in the SCS had almost no effect on E p,c .

14 SOLAR ENERGY↗

Tidal disruption of solitons in self-interacting ultralight axion dark matter

We report ultralight axions (ULAs) are promising dark matter candidates that can have a distinct impact on the formation and evolution of structure on nonlinear scales relative to the cold, collisionless dark matter (CDM) paradigm. However, most studies of structure formation in ULA models do not include the effects of self-interactions, which are expected to arise generically. Here, we study how the tidal evolution of solitons is affected by ULA self-interaction strength and sign. Specifically, using the pseudospectral solver UltraDark.jl, we simulate the tidal disruption of self-interacting solitonic cores as they orbit a 10 11 M ⊙ Navarro-Frenk-White CDM host halo potential for a range of orbital parameters, assuming a fiducial ULA particle mass of 10 -22 eV . We find that repulsive (attractive) self-interactions significantly accelerate (decelerate) soliton tidal disruption. We also identify a degeneracy between the self-interaction strength and soliton mass that determines the efficiency of tidal disruption, such that disruption timescales are affected at the ~ 50 % level for variations in the dimensionless ULA self-coupling from λ = -10 -92 to λ = 10 -92 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗