Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atmospheric physics parameters”

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

Neutrino oscillation phenomenology in the standard model and beyond

Over the last 20 years measurements of neutrino oscillation parameters have become very precise. In the standard neutrino oscillation picture most of the parameters are measured at the percent level. In this thesis we study neutrino oscillations in the standard picture and beyond. We analyze data from all types of neutrino oscillation experiments to obtain a global picture of neutrino oscillations. The remaining unknowns in the standard picture are the value of the CP-violating phase $\delta$, the octant of the atmospheric angle $\theta_{23}$ and the neutrino mass ordering. We discuss the current status of these unknowns and also comment on how well future experiments will do in measuring these quantities, where we discuss new facilities which will observe accelerator, atmossperic and reactor neutrinos. Using oscillation data it is possible to bound a possible violation of the CPT symmetry in the neutrino sector. This is done by performing a fit to neutrino and antineutrino oscillation data separately. We also address the capability of the future experiment DUNE to this scenario. We show that DUNE could improve some of the current bounds considerably. We discuss the sensitivity of DUNE together with the next generation reactor neutrino experiment JUNO to measure quasi-Dirac neutrino oscillations. We find that the future experiments will be able to bound the important quantities, mostly unbounded by current data. Finally, we discuss atmospheric neutrinos at the next generation neutrino telescope ORCA. We show that ORCA will put the best bounds on the invisible neutrino decay of $\nu_3$, coming from neutrino oscillation experiments. We also find that, a scenario like this would not affect the standard measurements of ORCA, namely the determination of the atmospheric neutrino oscillation parameters and the measurement of the neutrino mass ordering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Stability Analysis of Interface Conditions for Ocean-Atmosphere Coupling

In this paper, we analyze the stability of different coupling strategies for multidomain PDEs that arise in general circulation models used in climate simulations. We focus on fully coupled ocean–atmosphere models that are needed to represent and understand the complicated interactions of these two systems, becoming increasingly important in climate change assessment in recent years. Numerical stability issues typically arise because of different time-stepping strategies applied to the coupled PDE system. In particular, the contributing factors include using large time steps, lack of accurate interface flux, and single-iteration coupling. We investigate the stability of the coupled ocean–atmosphere models for various interface conditions such as the Dirichlet–Neumann condition and the bulk interface condition, which is unique to climate modeling. By analyzing a simplified model, we demonstrate here how the parameterization of the bulk condition and other numerical and physical parameters affect the coupling stability and establish stability conditions for different coupling strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Analysis of aerosol cloud interactions with a consistent signal of meteorology and other influencing parameters

Quantifying the impact of aerosols on cloud micro/macro physical properties and estimating the signature of Aerosol Cloud Interactions (ACI) is one of the challenging tasks in atmospheric sciences. The Moderate Resolution Imaging Spectroradiometer and the European Centre for Medium-Range Weather Forecasts ERA-5 reanalysis data are employed to systematically study the ACI over the monsoon region in Pakistan. Based on the monsoon occurrence and rainfall intensity, the whole region is divided into three sub-regions labeled as highly intensive (R1), moderately intensive (R2) and weak (R3) monsoon region. The results indicate that the monthly mean Aerosol Optical Depth (AOD) peaks in the summer monsoon months (Jun, Jul, Aug, Sep). Here, the well-known Twomey effect whereby the Cloud Droplet Radius (CDR) decreases with increasing AOD holds only over R3; the opposite effects (Anti-Twomey effect) are found over R1 and R2, all passing the test of statistical significance (p<0.05). The multi-year AOD is found to be positively correlated with Cloud Liquid Water Path (CLWP) and Cloud Optical Depth (COD) over R1 and R2, suggesting that thicker clouds containing more water droplets are formed in polluted atmosphere. Over R3, decreases in CLWP and COD are found with increasing AOD only when AOD is less than~0.325. The analysis of ACI over R1 and R2 during the winter months shows similar but stronger responses of CDR, CLWP and COD to the variation in AOD. The weaker responses during the summer monsoon season may attributed to the occurrence of high level cloud and unstable atmospheric condition. Further investigation of the influences of Relative Humidity and pressure vertical velocity on the CDR-AOD relationships shows that although the magnitude of the CDR-AOD correlations change with meteorological conditions, the sign of correlations remain unchanged with meteorological conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Physical Parameters of the Multiplanet Systems HD 106315 and GJ 9827

HD 106315 and GJ 9827 are two bright, nearby stars that host multiple super-Earths and sub-Neptunes discovered by K2 that are well suited for atmospheric characterization. We refined the planets’ ephemerides through Spitzer transits, enabling accurate transit prediction required for future atmospheric characterization through transmission spectroscopy. Through a multiyear high-cadence observing campaign with Keck/High Resolution Echelle Spectrometer and Magellan/Planet Finder Spectrograph, we improved the planets’ mass measurements in anticipation of Hubble Space Telescope transmission spectroscopy. For GJ 9827, we modeled activity-induced radial velocity signals with a Gaussian process informed by the Calcium II H&K lines in order to more accurately model the effect of stellar noise on our data. We measured planet masses of M {sub b} = 4.87 ± 0.37 M {sub ⊕}, M {sub c} = 1.92 ± 0.49 M {sub ⊕}, and M {sub d} = 3.42 ± 0.62 M {sub ⊕}. For HD 106315, we found that such activity radial velocity decorrelation was not effective due to the reduced presence of spots and speculate that this may extend to other hot stars as well (T {sub eff} > 6200 K). We measured planet masses of M {sub b} = 10.5 ± 3.1 M {sub ⊕} and M {sub c} = 12.0 ± 3.8 M {sub ⊕}. We investigated all of the planets’ compositions through comparison of their masses and radii to a range of interior models. GJ 9827 b and GJ 9827 c are both consistent with a 50/50 rock-iron composition, GJ 9827 d and HD 106315 b both require additional volatiles and are consistent with moderate amounts of water or hydrogen/helium, and HD 106315 c is consistent with a ∼10% hydrogen/helium envelope surrounding an Earth-like rock and iron core.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

The role of source geometry and atmospheric propagation in global bolide infrasound detectability

Global infrasound monitoring provides a persistent means of detecting energetic bolide atmospheric entries, complementing optical observations and extending coverage over remote regions. We present a global assessment of the physical factors governing bolide infrasound detectability by correlating 623 bolide events reported by the Center for Near-Earth Object Studies between 2007 and 2025 with waveform data from the International Monitoring System. We identify 311 events with confirmed infrasound detections, corresponding to a detection rate of approximately 50%, substantially higher than inferred from earlier surveys, reflecting both the maturation of the global infrasound network and advances in automated, multi-frequency array processing. Analysis of flight parameters shows that infrasound detectability is selective rather than uniform across the bolide population. Detected events are preferentially associated with steeper entry angles and lower-altitude energy deposition, while shallow, high-altitude trajectories are less consistently observed. Very high-energy events remain detectable regardless of geometry, but for the more common lower-energy regime, observability depends on specific combinations of entry parameters and propagation conditions. This geometric dependence persists across comparable energy ranges and atmospheric conditions, indicating that entry angle exerts a primary control on detectability, with energy and propagation acting as secondary modulating factors. Furthermore, these results provide new physical constraints on bolide-atmosphere interactions and improve interpretation of global infrasound observations for planetary defense and atmospheric-entry studies.

Bolides↗
DOE OSTI2022

New Particle Formation and Growth Dynamics for α-Pinene Ozonolysis in a Smog Chamber and Implications for Ambient Environments

Secondary organic aerosol (SOA) plays an important role in the Earth’s energy balance and air quality. Understanding the formation of SOA accompanying the oxidation of volatile organic compounds (VOCs) remains incomplete because of the complex reactions and partitioning mechanisms. In this study, we investigated new particle formation (NPF) of SOA derived from α-pinene ozonolysis in the Harvard Environmental Chamber using the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the condensing species as two products, one of lower volatility (LVOC, 0.01 μg m –3 ) and one of higher volatility (SVOC, 1.0 μg m –3 ). A module for the classical nucleation theory was included. The surface tension (σ), bulk diffusivity (D b ), and mass accommodation coefficient (α) were varied until convergence between the simulated and observed particle number size distributions. α of LVOC was constrained as 0.03–0.1 with the reported σ of α-pinene-ozonolysis-derived SOA and Tolman surface correction (~23.0–27.5 mN m –1 ) and Db from 10 –15.5 to >10 –11 cm 2 s –1 . The observed particle number size distribution of α-pinene ozonolysis was well-simulated to derive the physical parameters, including σ, Db, and α, suggesting the possible application to group the reported complex oxidation products and further application for other atmospheric-related oxidations for quantifying NPF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multiscale Reduced Order Modeling and Parameter Estimation for Climate Sciences

Several problems in earth system modeling are dependent on highly multiscale phenomena, such as turbulence, where computational modeling is challenging and expensive. This issue is exacerbated in atmospheric and oceanic domains, due to inherent high-dimensionality of the problem. One approach to this problem has been reduced order modeling (ROM); which aims to represent the key physics of the phenomena as a low-dimensional system. AI methods have huge potential in building accurate, stable ROMs and parameter estimation for these ROMs, as it requires extracting nonlinearities and patterns from simulation and/or observational data. Developing physics-based AI approaches specialized for the complexities of multiscale data, along with strategies to account for uncertainties, will revolutionize rapid modeling, analysis and decision making for earth system problems of practical interest.

58 GEOSCIENCES↗
DOE OSTI2023

Supporting Advancement in Weather and Water Prediction in the Upper Colorado River Basin: The SPLASH Campaign

Water is a critical resource that causes significant challenges to inhabitants of the western United States. These challenges are likely to intensify as the result of expanding population and climate-related changes that act to reduce runoff in areas of complex terrain. To better understand the physical processes that drive the transition of mountain precipitation to streamflow, the National Oceanic and Atmospheric Administration has deployed suites of environmental sensors throughout the East River watershed of Colorado as part of the Study of Precipitation, the Lower Atmosphere, and Surface for Hydrometeorology (SPLASH). This includes surface-based sensors over a network of five different observing sites, airborne platforms, and sophisticated remote sensors to provide detailed information on spatiotemporal variability of key parameters. With a 2-yr deployment, these sensors offer detailed insight into precipitation, the lower atmosphere, and the surface, and support the development of datasets targeting improved prediction of weather and water. Initial datasets have been published and are laying a foundation for improved characterization of physical processes and their interactions driving mountain hydrology, evaluation and improvement of numerical prediction tools, and educational activities. SPLASH observations contain a depth and breadth of information that enables a variety of atmospheric and hydrological science analyses over the coming years that leverage collaborations between national laboratories, academia, and stakeholders, including industry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Aircraft in-situ measurements from SOCRATES constrain the anthropogenic perturbations of cloud droplet number

Aerosol-cloud interactions (ACI) in warm clouds alter reflected shortwave radiation by influencing cloud microphysical and macrophysical properties. The variable of state controlling ACI is the cloud droplet number concentration (N d ). Here, we examine the perturbations in N d due to anthropogenic aerosols (ΔN d, PD-PI ) using a perturbed parameter ensemble (PPE) hosted in the sixth Community Atmosphere Model (CAM6). Surrogate models are created for the CAM6 PPE outputs and are used to generate 1 million model variants of CAM6 by sampling 45 sources of parameter uncertainty. The range of uncertain physical parameters related to ACI are constrained with observations of aerosol and cloud properties from SOCRATES. The likely range of uncertain parameters and the associated range of ΔN d, PD-PI are more strongly constrained with observations of N d relative to observations of cloud condensation nuclei. We conduct sensitivity tests of how constraints on ΔN d, PD-PI are affected by systematic uncertainties in observations and our limitations in our surrogate models created for CAM6 PPE outputs. Based on this, we provide guidance on the impact of reducing systematic uncertainty in airborne microphysical observations and in surrogate models.

Song, Ci [Univ. of Wyoming, Laramie, WY (United St↗
DOE OSTI2021

Fully quantum calculations of O 2 –N 2 scattering using a new potential energy surface: Collisional perturbations of the oxygen 118 GHz fine structure line

A proper description of the collisional perturbation of the shapes of molecular resonances is important for remote spectroscopic studies of the terrestrial atmosphere. Of particular relevance are the collisions between the $\text{O}_2$ and $\text{N}_2$ molecules—the two most abundant atmospheric species. In this work, we report a new highly accurate $\text{O}_2(\text{X}^3Σ^–_g)–\text{N}_2(\text{X}^1Σ^+_g)$ potential energy surface and use it for performing the first quantum scattering calculations addressing line shapes for this system. We use it to model the shape of the 118 GHz fine structure line in $\text{O}_2$ perturbed by collisions with $\text{N}_2$ molecules, a benchmark system for testing our methodology in the case of an active molecule in a spin triplet state. The calculated collisional broadening of the line agrees well with the available experimental data over a wide temperature range relevant for the terrestrial atmosphere. This work constitutes a step toward populating the spectroscopic databases with ab initio line shape parameters for atmospherically relevant systems.

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

Hunting extreme BL Lacertae blazars with Fermi-Large Area Telescope

ABSTRACT The emission of very-high-energy (VHE) photons ($E\ \gt\ 100\ \mathrm{GeV}$) in blazars is closely connected to the production of ultra-relativistic particles and the role of these γ-ray sources as cosmic particle accelerators. This work focuses on a selection of 22 γ-ray objects from the 2BIGB catalogue of high-synchrotron-peaked sources, which are classified as blazar candidates of uncertain type in the 4FGL-DR2 catalogue. We study these sources by means of a re-analysis of the first 10 yr of γ-ray data taken with the Fermi-Large Area Telescope, including the attenuation by the extragalactic background light. Their broad-band spectral energy distributions are also evaluated, using multiwavelength archival data in the radio, optical, and X-ray bands, in terms of one-zone synchrotron-self-Compton models, adding an external Compton component when needed. Out of this analysis, we identify 17 new extreme high-synchrotron-peaked (EHSP) candidates and compare their physical parameters with those of prototypical EHSP blazars. Finally, the resulting models are used to assess their detectability by the present and future generation of ground-based Imaging Atmospheric Cherenkov Telescopes. We find two VHE candidates within the reach of the current and next generation of Cherenkov telescopes: J0847.0−2336 and J1714.0−2029.

Nievas Rosillo, M. (ORCID:0000000283219168)↗
DOE OSTI2021

Spectroscopic and Dynamical Studies of Highly Energized Small Polyatomic Molecules (End-of-Grant Report)

The aim of this research program has been to develop models for qualitative and quantitative understanding of complex processes in small molecules, such as vibronic coupling, predissociation, photodissociation, and isomerization. Simple models capture physical insights and lead to intuitive ideas that can enlighten a broad array of problems. Numerical fits can provide predictive power, but rarely impart insights into the underlying mechanisms responsible for the observed phenomena. This DOE supported program has prioritized a quantum mechanical understanding over a perfect numerical fit of spectra. Finding few-parameter models, when those parameters have physical meaning, allows the insights gained from one problem to be easily applied to related problems in other systems. The systems that have been under investigation include atmospherically relevant sulfur-containing compounds, various unimolecular dynamical processes of electronically excited acetylene, and the rovibrational dynamics of formaldehyde and SO 2 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Development of a triple-moment ice-phase cloud microphysics scheme and its application to the Single Column Atmosphere Model

Parameterization of cloud microphysics is critical for accurate simulation of weather and climate, in which the characteristics of cloud particle spectrum inevitably further affect climate simulation by changing the cloud evolution and cloud radiation effects. The popular currently used double-moment cloud microphysics schemes in numerical models can predict the intercept (N 0 ) and slope (λ) parameters of cloud particle spectrum but cannot predict the spectral shape parameter (μ), which hinders accurate description of cloud physical processes in climate models. Therefore, in the present study, we built upon the ideas of previously developed triple-moment cloud microphysics scheme, considered radar reflectivity factor as the third predictor in addition to number and mass concentration, and introduced its related prediction equation to the Single Column Atmosphere Model Version 5.3 (SCAM5.3). Moreover, the relevant microphysical process formulas in the model were revised, and a triple-moment ice-cloud microphysics scheme was constructed to predict the μ of ice particle spectrum. Based on this model, a 29-day case of the Atmospheric Radiation Measurement Program in the summer of 1997 (ARM97) was simulated, and the differences in cloud fraction and radiation simulations between the double- and triple-moment schemes were analysed. The μ of ice particle spectrum predicted by the triple-moment scheme mainly ranged from 0 to 4; the peak value was around 2, and at least 85% of the μ values were greater than the default 0. Therefore, the μ=0 setting in the double-moment scheme is unreasonable. Moreover, the developed triple-moment ice-cloud microphysics scheme yielded a narrower ice particle spectrum, which was closer to the observation results of previous studies, than the double-moment cloud microphysics scheme (e.g., height=233 hPa). Furthermore, compared with the double-moment scheme, the triple-moment scheme achieved closer simulations of the cloud fraction observations, particularly for ice clouds in the upper levels. Considering the close association between cloud fraction and radiation, the error with the triple-moment scheme was smaller than that with the double-moment scheme regardless of the shortwave (surface downward shortwave radiation, net downward shortwave radiation at the top of the atmosphere, and net shortwave radiation at the surface) or longwave (surface downward longwave radiation and net longwave radiation at the surface) flux density. Finally, the improvement mechanism of the triple-moment scheme on cloud fraction and radiation simulations was explored. The major reason is that the triple-moment scheme weakens the autoconversion of ice crystal to form snow and enhances the growth process of ice crystal deposition, thereby increasing the mass concentration of ice crystals and ice cloud fraction. Overall, the developed triple-moment ice-cloud microphysics scheme can help improve the simulation ability of models for cloud-climate feedback and other processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

CPT violation sensitivity of NoVA, T2K and INO experiments using $\nu$ and $\bar{\nu}$ oscillation parameters

Charge-Parity-Time (CPT) symmetry allows only identical oscillation parameters for $\nu$ and $\bar{\nu}$. But,different mass and mixing parameters for $\nu$ and $\bar{\nu}$ can give us possible hint for CPT violation or new physics. Using, different oscillation parameters for $\nu$ and $\bar\nu$, we find sensitivity for ($\Delta m^{2}_{32}-\Delta\bar{m}^{2}_{32}$) and ($\sin^{2}\theta_{23}-\sin^{2}\bar{\theta}_{23}$) for long-baseline (T2K and NOvA) and atmospheric neutrino (INO) experiments in different possible combinations of octant for neutrinos and anti-neutrinos. We present the joint sensitivity of the T2K, NOvA and INO experiments to such CPT violating observables.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Physics of parameter correlations around the solar-scale enhancement in neutrino theory with unitarity violation

We discuss the physics of the three neutrino flavor transformation with non-unitary mixing matrix, with particular attention to the correlation between the $\nu$SM- and the $\alpha$ parameters which represent the effect of unitarity-violating (UV) new physics. Towards this goal, a new perturbative framework is created to illuminate the effect of non-unitarity in the region of the solar-scale enhanced oscillations. We refute the skepticism about the physical reality of the $\nu$Standard Model CP phase $\delta$–$\alpha$ parameter phase correlation by analysis with the SOL convention of $U_{{\tiny MNS}}$, in which $e^{\pm i \delta}$ is attached to $s_{12}$. Then, a comparative study between the solar- and atmospheric-scale oscillation regions allowed by the framework reveals a dynamical $\delta$–(blobs of the $\alpha$ parameters) correlation in the solar oscillation region, in sharp contrast to the “chiral”-type phase correlation $[e^{- i \delta} \bar{\alpha}_{\mu e},\ e^{- i \delta} \bar{\alpha}_{\tau e},\ \bar{\alpha}_{\tau \mu}]$ in the Particle Data Group convention seen in the atmospheric oscillation region. An explicit perturbative calculation to the first order in the $\nu_{\mu} \rightarrow \nu_{e}$ channel allows us to decompose the UV related part of the probability into the unitary evolution part and the genuine non-unitary part. We observe that the effect of non-unitarity tends to cancel between these two parts, as well as between the different $\alpha_{\beta \gamma}$ parameters.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Tidal Love numbers of novel and admixed celestial objects

A sub-fraction of dark matter or new particles trapped inside celestial objects can significantly alter their macroscopic properties. Here, we investigate the new physics imprint on celestial objects by using a generic framework to solve the Tolman-Oppenheimer-Volkoff (TOV) equations for up to two fluids. We test the impact of populations of new particles on celestial objects, including the sensitivity to self-interaction sizes, new particle mass, and net population mass. Applying our setup to neutron stars and boson stars, we find rich phenomenology for a range of these parameters, including the creation of extended atmospheres. These atmospheres are detectable by their impact on the tidal love number, which can be measured at upcoming gravitational wave experiments such as Advanced LIGO, the Einstein Telescope, and LISA. We release our calculation framework as a publicly available code at this URL, allowing the TOV equations to be generically solved for arbitrary new physics models in novel and admixed celestial objects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗