Search NASA⌕ Search

SEARCH · Search NASA

Results for “STABILITY CALCULATIONS”

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

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Active Learning Accelerated Discovery of Stable Iridium Oxide Polymorphs for the Oxygen Evolution Reaction

The discovery of high-performing and stable materials for sustainable energy applications is a pressing goal in catalysis and materials science. Understanding the relationship between a material’s structure and functionality is an important step in the process, such that viable polymorphs for a given chemical composition need to be identified. Machine-learning-based surrogate models have the potential to accelerate the search for polymorphs that target specific applications. Herein, we report a readily generalizable active-learning (AL) accelerated algorithm for identification of electrochemically stable iridium oxide polymorphs of IrO 2 and IrO 3 . The search is coupled to a subsequent analysis of the electrochemical stability of the discovered structures for the acidic oxygen evolution reaction (OER). Structural candidates are generated by identifying all 956 structurally unique AB2 and AB3 prototypes in existing materials databases (more than 38000). Next, using an active learning approach, we find 196 IrO 2 polymorphs within the thermodynamic amorphous synthesizability limit and reaffirm the global stability of the rutile structure. We find 75 synthesizable IrO 3 polymorphs and report a previously unknown FeF 3 -type structure as the most stable, termed α-IrO 3 . To test the algorithms performance, we compare to a random search of the candidate space and report at least a 2-fold increase in the rate of discovery. Additionally, the AL approach can acquire the most stable polymorphs of IrO 2 and IrO 3 with fewer than 30 density functional theory optimizations. Analysis of the structural properties of the discovered polymorphs reveals that octahedral local coordination environments are preferred for nearly all low-energy structures. Subsequent Pourbaix Ir–H 2 O analysis shows that α-IrO3 is the globally stable solid phase under acidic OER conditions and supersedes the stability of rutile IrO 2 . Calculation of theoretical OER surface activities reveal ideal weaker binding of the OER intermediates on α-IrO 3 than on any other considered iridium oxide. We emphasize that the proposed AL algorithm can be easily generalized to search for any binary metal oxide structure with a defined stoichiometry.

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

Single Ru Atoms Stabilized by Hybrid Amorphous/Crystalline FeCoNi Layered Double Hydroxide for Ultraefficient Oxygen Evolution

Abstract In view of the sluggish kinetics suppressing the oxygen evolution reaction (OER), developing efficient and robust OER catalysts is urgent and essential for developing efficient energy conversion technologies. Herein, hybrid amorphous/crystalline FeCoNi layered double hydroxide (LDH)‐supported single Ru atoms (Ru SAs/AC‐FeCoNi) are developed for enabling a highly efficient electrocatalytic OER. The amorphous outer layer in Ru SAs/AC‐FeCoNi is composed of abundant defect sites and unsaturated coordination sites, which can serve as anchoring sites to stabilize single Ru atoms. The crystalline inner has a highly symmetric rigid structure, thereby strengthening the stability of support for a long‐lasting OER. The synergistic effects endow this hybrid catalyst with extremely low overpotential (205 mV at 10 mA cm −2 ). Density functional theory calculation indicates that single Ru atoms stabilized by hybrid amorphous/crystalline FeCoNi LDH facilitate the formation of Ru–O* (rate‐determining step), thus accelerating the OER process.

Hu, Yidong↗
DOE OSTI2025

Polymorphs of the n–Type Polymer P(NDI2OD–T2): A Comprehensive Description of the Impact of Processing on Crystalline Morphology and Charge Transport

A systematic study of the polymorphs emerging in P(NDI2OD-T2) (also commercially known as N2200), a prototypical organic semiconducting n-type polymer, is presented. Using a tightly integrated experimental and computational approach, detailed atomistic-level descriptions are provided investigating the three known P(NDI2OD-T2) polymorphs observed at room temperature as a function of thin-film processing. Importantly, over the course of the work, a missing link is uncovered, a fourth polymorph referred to here as Form I-β; this new form is a morphological intermediary observed upon thermal annealing, which evolves from Form I but tends to disappear upon full polymer chain melting. The computationally derived polymorph structures show excellent agreement with experimental X-ray scattering characterization. The relative stabilities of each polymorph are calculated in terms of both the bulk material and the polymorph-air interface. An energy landscape is then constructed to qualitatively compare the thermodynamic versus kinetic origins of each polymorph, and the factors driving (supra)assembly and associated transformations among polymorphs using an approach generalizable to other organic semiconducting polymers. Lastly, the relationships among preferential polymorphic crystallinity, relative chain orientations, and directional charge transport properties in P(NDI2OD-T2) are explored. Altogether, this work provides unprecedented insights into complex structure-processing-transport relationships in a representative semiconducting organic polymer.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of the Nature of Both Cation and Anion Substitution on the Structural Symmetry of Li‐Rich 3 d ‐Metal Chalcogenide Electrodes

Abstract Li‐rich layered chalcogenides have recently led to better understanding of the anionic redox process and its associated high capacity while providing ways to overcome its practical limitations of voltage fade and irreversibility. This study reports on the feasibility of triggering anionic activity in Li 2 TiS 3 , through anionic substitution (Se for S) or cationic substitution (Fe for Ti). Herein, the chalcogenide chemical space is further explored to prepare mono‐substituted Li 1.7 Ti 0.85 Mn 0.45 Ch 3 (Ch = S/Se) and doubly substituted cationic and anionic phases (Li 1.7 Ti 0.85 Fe 0.45 S 3‐z Se z ) which crystallize either in the O3‐ or O1‐type structures depending upon substituents. All series show a bell‐shape capacity variation as function of the transition metal (TM) substitution degree with values up to 240 mAh g −1 . For specific compositions, a structural O3 to O1 phase transition is observed upon Li removal, which is not reversible upon Li re‐insertion due to kinetic limitations and negatively affects long‐term cycling performance. Density functional theory (DFT) calculations confirm the O3/O1 relative stability along the different series and point subtle electronic differences in the TM‐doping, rationalizing the structural and electrochemical behaviors of these phases upon cycling. These findings provide further insights into the link between structural and electronic stability, which is of key importance for designing chalcogenide‐based anionic redox compounds.

Chemistry↗
DOE OSTI2021

Role of Fluorine in Chemomechanics of Cation-Disordered Rocksalt Cathodes

The role of nanoscale chemomechanical behavior in the macroscopic performance of functional materials is well recognized. For lithium-ion battery cathodes, tremendous effort has been devoted to the development of new chemistry and structure, yet fundamental understanding of the correlation between redox processes and mechanical properties of the novel materials lags behind. In the present study, we prepare large discrete single grains of Li-excess cation-disordered rocksalts (DRX) and investigate their chemomechanical behavior at the particle level, using nanoresolution X-ray and electron-based spectro-imaging and chemical mapping techniques. While irregular cracking upon lithium extraction leads to the eventual breakdown of the baseline DRX oxide (Li 1.2 Ti 0.4 Mn 0.4 O 2 ) particles at a high delithiation state, the fluorinated-DRX (Li 1.3 Ti 0.3 Mn 0.4 O 1.7 F 0.3 ) clearly displays aligned cracking along the <001> direction. The resulting periodicity in the cracking pattern enables the particles to retain their integrity and, consequently, improved electrochemical stability. Density functional theory (DFT) calculations showed that fluorination leads to increased concentration of Li+ on the (001) planes and preferential Li movements along the <001>-family directions, revealing the underlying mechanism for directional cracking. Overall, our study demonstrates the unique role of fluorine in modulating nanoscale chemomechanics, which in turn influences the evolution of charge and strain heterogeneity at the particle level. These insights provide important design guidelines in further improving DRX cathode materials.

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

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗
DOE OSTI2021

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond (resistance) and recover from (resilience) to rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the ED-2.2 model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900, and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes, and interacted with one another. We then calculated metrics of functional ecosystem stability. Modeled gross primary production exhibited a mean resistance of -0.13, -0.23, and -0.42 (unitless) in response to increasing disturbance severity of 45%, 65%, and 85% mortality, respectively; and resilience of 0.002, 0.004, 0.008 (yr-1) respectively. This pattern held true for other carbon fluxes, meaning that in general simulated GPP declined more initially but recovered more quickly from higher mortality. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Explainable tokamak-agnostic forecasting of fusion plasma instability via megahertz turbulent fluctuations

Scientific applications of artificial intelligence (AI) often remain limited by device-specific training and unexplained “black-box” approaches, creating fundamental barriers to cross-system generalization. This challenge is critical for nuclear fusion, where future reactors will have limited operational data for AI training. Here, we demonstrate that our neural network, trained solely on megahertz-scale turbulence measurements from one machine (DIII-D), forecasts Type-I edge localized mode (ELM) onsets in a different tokamak (KSTAR) through zero-shot weight transfer following physics-consistent preprocessing without device-specific retraining. Through an explainable AI framework combining gradient-weighted class activation mapping with physics validation, we reveal that our network can internalize physics relationships governing the ELM instabilities rather than memorizing device-specific patterns. The network perceives spatiotemporal features that correlate consistently with independently calculated instability growth rates, magnetohydrodynamic stability limits, and pedestal structure dynamics. Statistical analyses of dimensionally-reduced saliency features reveal the identical triangular features between the saliency representations, instability growth rates, and prediction probability across tokamaks, providing evidence that our forecasting system can show tokamak-agnostic generalization. This work contributes to a foundation for explainable scientific AI systems, where cross-system developments are essential for transcending traditional domain-specific constraints.

AI↗
DOE OSTI2020

Determinants for proton location and electron coupled proton transfer in hydrogen bonded pentafluorophenol-anion clusters

This work aims to reveal the determinant factors for proton locations and electron coupled proton transfer (ECPT) in biologically relevant hydrogen bonded systems. Pentafluorophenol-anion clusters [C6F5O-•H+•A-]- are chosen to model active sites of biological functions, with the anion A- being systematically varied to ensure the proton affinities (PAs) of the anions well cover the referenced PA of C6F5O- from being appreciably smaller, similar, and to significantly larger. Negative ion photoelectron spectroscopy of these clusters provides spectroscopic evidence showing that proton location in the anionic state is largely following the PA prediction, while ECPT is observed only for the clusters with the anion possessing the electron binding energy (EBE) significantly larger than that of the referenced C6F5O-. Theoretical calculations suggest these clusters are stabilized by forming single strong hydrogen bond between donor and acceptor, and the associated charge and MO analyses fully support the experimental observations. The current holistic cluster model study indicates that PA is the right determinant that can be used to predict the proton location and describe hydrogen bonding structures, while both PA and EBE of anionic groups play important roles in facilitating ECPT process. This work was supported by U.S. Department of Energy (DOE), Office of Basic Energy Sciences, Division of Chemical Science, Geosciences, and Biosciences (X.-B.W.), and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. Part of this work was supported by the National Natural Science Foundation of China (No-11604046 for JZ; No. 11727810 for ZRS), as well as “Chenguang Program” by Shanghai Education Development Foundation, Shanghai Municipal Education Commission (No. 16CG38) and Natural Science Foundation of Shanghai (No. 16ZR1448100) (JZ).

Zhang, Jian↗
DOE OSTI2023

Solvent effects on extractant conformational energetics in liquid–liquid extraction: a simulation study of molecular solvents and ionic liquids

Extractant design in liquid–liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant–metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. Here, we calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

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

Reliable operation of Cr 2 O 3 :Mg/ $β$-Ga 2 O 3 p–n heterojunction diodes at 600 °C

Beta gallium oxide (β-Ga 2 O 3 )-based semiconductor heterojunctions have recently demonstrated improved performance at high voltages and elevated temperatures and are, thus, promising for applications in power electronic devices and harsh environment sensors. However, the long-term reliability of these ultra-wideband gap (UWBG) semiconductor devices remains barely addressed and may be strongly influenced by chemical reactions at the p–n heterojunction interface. Here, we experimentally demonstrate operation and evaluate the reliability of Cr 2 O 3 :Mg/β-Ga 2 O 3 p–n heterojunction diodes during extended operation at 600 °C, as well as after 30 repeated cycles between 25 and 550 °C. The calculated pO 2 -temperature phase stability diagram of the Ga-Cr-O material system predicts that Ga 2 O 3 and Cr 2 O 3 should remain thermodynamically stable in contact with each other over a wide range of oxygen pressures and operating temperatures. The fabricated Cr 2 O 3 :Mg/β-Ga 2 O 3 p–n heterojunction diodes show room-temperature on/off ratios >104 at ±5 V and a breakdown voltage (V Br ) of -390 V. The leakage current increases with increasing temperature up to 600 °C, which is attributed to Poole–Frenkel emission with a trap barrier height of 0.19 eV. Over the course of a 140-h thermal soak at 600 °C, both the device turn-on voltage and on-state resistance increase from 1.08 V and 5.34 mΩ cm 2 to 1.59 V and 7.1 mΩ cm 2 , respectively. This increase is attributed to the accumulation of Mg and MgO at the Cr 2 O 3 /Ga 2 O 3 interface as observed from the time-of-flight secondary ion mass spectrometry analysis. These findings inform future design strategies of UWBG semiconductor devices for harsh environment operation and underscore the need for further reliability assessments for β-Ga 2 O 3 -based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Metallicity and anomalous Hall effect in epitaxially strained, atomically thin RuO 2 films

The anomalous Hall effect (AHE), a hallmark of time-reversal symmetry breaking, has been reported in rutile RuO 2 , a debated metallic altermagnetic candidate. Previously, AHE in RuO 2 was observed only in strain-relaxed thick films under extremely high magnetic fields (~50 T). Yet, in ultrathin strained films with distinctive anisotropic electronic structures, there are no reports, likely due to disorder and defects suppressing metallicity thus hindering its detection. Here, we demonstrate that ultrathin, fully strained 2 nm TiO 2 /t nm RuO 2 /TiO 2 (110) heterostructures, grown by hybrid molecular beam epitaxy, retain metallicity and exhibit a sizeable AHE at a significantly lower magnetic field (< 9 T). Density functional theory calculations reveal that epitaxial strain stabilizes a noncompensated magnetic ground state and reconfigures magnetic ordering in RuO 2 (110) thin films. These findings establish ultrathin RuO 2 as a platform for strain-engineered magnetism and underscore the transformative potential of epitaxial design in advancing spintronic technologies.

MBE↗
DOE OSTI2023

PFPO plasma scenarios for exploration of long pulse operation in ITER

Long Pulse Scenarios (LPS) in ITER foreseen during the Pre-Fusion Power Operation (PFPO) phase of the ITER Research Plan (IRP) are assessed using 1.5D transport simulations within the ASTRA framework. Such assessment is required to predict the operational space for LPS operation in PFPO, as well as to evaluate which physics processes for LPS operation during Fusion Power Operation (FPO) could be studied during PFPO. An important aspect in the development of LPSs in PFPO is to minimize lifetime consumption of the Central Solenoid (CS) for these scenarios. The maximum pulse length achievable for LPSs in PFPO with no consumption of CS lifetime (currents in CS coils $&LessSlantEqual;$30 kA per turn) has been assessed for a range of heating schemes and heating mixes, confinement regimes (L-mode and H-mode) and for helium and hydrogen plasmas. The operational space of LPS and pulse length has been explored through density scans with the Heating and Current Drive mix required for the FPO Q $&GreaterSlantEqual;$ 5 steady-state plasma scenario (namely Neutral Beam Injection and Electron Cyclotron Heating) including acceptable shine through losses on the first wall for both helium and hydrogen plasmas. Fast particle physics aspects that are common between FPO plasmas and LPS PFPO H-mode plasmas at low densities are studied including MHD stability analysis with the KINX code and non-perturbative critical gradient model based on high-n Toroidal Alfven Eigenmodes (TAE) stability kinetic ballooning code HINST calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Design principles and physical properties of two-dimensional heterostructured borides

Principles of design to create dynamically stable transition metal-, lanthanide-, and actinide-based low-dimensional borides are presented. A charge transfer analysis of donor metal atoms to electron deficient honeycombed B lattices allows to predict complex covalent heterostructures hosting Dirac states. The applicable guidelines are supported with the analysis of phonon spectra computed with first-principles calculations to demonstrate the physical stability of nanometer-thick heterostructures. Similar or dissimilar layered borides can be stacked on top of each other in a layer-by-layer fashion creating compounds that can be fundamentally different from the individual layers, opening a rich playground to explore novel physical properties and new materials. Functionalities such as multiple Dirac states, highly dispersive electronic bands, and decoupled acoustic-optical phonon branches are studied. The combination of appealing electronic properties and physical realization makes of predicted layered borides promising materials to integrate a new generation of two-dimensional materials.

2-dimensional systems↗
DOE OSTI2021

Silicon Nanosheets as Candidates for Silicon-Based Optoelectronics

The growth in computational ability over the past decades has positively impacted global development, the economy, healthcare, and science. As on-chip components are approaching the atomic scale, alternative paradigms are needed to address the thermal and electronic issues that impose bottlenecks for computing. One approach to address this is with optoelectronics. However, silicon—the backbone of microelectronics—is a poor choice due to its indirect bandgap, while existing optoelectronic materials are incompatible with CMOS infrastructure. Monolayer silicon nanosheets (SiNSs) are an intriguing material that exhibit photoluminescence, and are compositionally-compatible with the CMOS process. Here, we synthesize and characterize monolayer SiNSs, and show spectroscopic evidence that they exhibit a quasi-direct bandgap, which is corroborated by DFT calculations. We probe their thermal stability, demonstrating their structure and photoluminescence are stable beyond the required operating temperatures for computing applications. These optoelectronic properties, CMOS-compatibility, and stability make SiNSs a viable candidate for silicon-based photonics.

Ryan, Bradley↗
DOE OSTI2024

FAR3d

The FAR3d model calculates the linear and nonlinear stability properties of energetic particle driven Alfven instabilities for both tokamak and stellarator plasma confinement devices using gyro-landau closure methods. This is an important fundamental physics problem for existing fusion energy experiments and for future fusion reactors.

Varela, Jacobo↗
DOE OSTI2020

Revealing the Formation Energy–Exfoliation Energy–Structure Correlation of MAB Phases Using Machine Learning and DFT

MAB phases are popular as high-temperature materials with high damage tolerance and excellent electrical conductivity that are used to exfoliate two-dimensional (2D) transition metal borides (MBenes); a promising material for developing next- generation nanodevices. In this report, we explore the correlation between the formation energy, exfoliation energy and structural factors of MAB phases using DFT and machine learning. In order to do so, we developed three different machine learning models based on the support vector machine, deep neural network, and random forest regressor to study the stability of the MAB phases by calculating their formation energies. Furthermore, our support vector machine and deep neural network models are capable of predicting the formation energies with mean absolute errors less than 0.1 eV/atom. We demonstrated that the stability of a MAB phase for a given transition metal decreases when the A element changes from Al to Tl. The density functional theory (DFT) re- vealed that M-A and M-B bond strength strongly correlates with the stability of MAB phases. In addition, the exfoliation possibility of 2D MBenes becomes higher when the A element changes from Al to Tl due to weakening of M-A and M-B bonds.

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

Scalar dark matter vortex stabilization with black holes

Galaxies and their dark-matter halos are commonly presupposed to spin. But it is an open question how this spin manifests in halos and soliton cores made of scalar dark matter (SDM, including fuzzy/wave/ultralight-axion dark matter). One way spin could manifest in a necessarily irrotational SDM velocity field is with a vortex. But recent results have cast doubt on this scenario, finding that vortices are generally unstable except with substantial repulsive self-interaction. Here, in this paper, we introduce an alternative route to stability: in both (non-relativistic) analytic calculations and simulations, a black hole or other central mass at least as massive as a soliton can stabilize a vortex within it. This conclusion may also apply to AU-scale halos bound to the sun and stellar-mass-scale Bose stars.

79 ASTRONOMY AND ASTROPHYSICS↗