Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LaTe”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Blast from the past: constraining progenitor models of SN 1972E

ABSTRACT We present a novel technique to study Type Ia supernovae (SNe Ia) by constraining surviving companions of historical extragalactic SN by combining archival photographic plates and Hubble Space Telescope(HST) imaging. We demonstrate this technique for Supernova 1972E, the nearest known SN Ia in 125 yr. Some models of SNe Ia describe a white dwarf with a non-degenerate companion that donates enough mass to trigger thermonuclear detonation. Hydrodynamic simulations and stellar evolution models show that these donor stars will survive the explosion, and show increased luminosity for at least a 1000 yr. Thus, late-time observations of the exact location of a supernova can constrain the presence of a surviving donor star and progenitor models. We find the explosion site of SN 1972E by analysing 17 digitized photographic plates taken with the European Southern Observatory 1-m Schmidt and 1 plate taken with the Cerro Tololo Inter-American Observatory 1.5-m telescope. Using the Gaia eDR3 catalogue to determine Supernova 1972E’s location yields: α = 13h39m52${_{.}^{\rm s}}$708 ± 0${_{.}^{\rm s}}$004 and δ = −31°40’9${_{.}^{\prime\prime}}$00 ± 0${_{.}^{\prime\prime}}$04 (ICRS). In 2005, HST/ACS imaged the host galaxy of SN 1972E with the F435W, F555W, and F814W filters covering the explosion site. The nearest detected source is offset by 3.0 times our positional precision, and is inconsistent with the colours expected of a surviving donor star. Thus, the limiting magnitude of the HST observation (F555W > 28 mag) rules out all He star companion models and the most luminous main-sequence companion model currently in the literature. The remaining main-sequence companion models could be tested with a 10 orbit HST exposure in the F606W filter.

Do, Aaron↗
DOE OSTI2021

Progenitor and close-in circumstellar medium of type II supernova 2020fqv from high-cadence photometry and ultra-rapid UV spectroscopy

ABSTRACT We present observations of SN 2020fqv, a Virgo-cluster type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 d post-explosion; ground-based spectroscopic observations starting 1.1 d post-explosion; along with extensive photometric observations. Massive stars have complicated mass-loss histories leading up to their death as CCSNe, creating circumstellar medium (CSM) with which the SNe interact. Observations during the first few days post-explosion can provide important information about the mass-loss rate during the late stages of stellar evolution. Model fits to the quasi-bolometric light curve of SN 2020fqv reveal 0.23 M⊙ of CSM confined within 1450 R⊙ (1014 cm) from its progenitor star. Early spectra (<4 d post-explosion), both from HST and ground-based observatories, show emission features from high-ionization metal species from the outer, optically thin part of this CSM. We find that the CSM is consistent with an eruption caused by the injection of ∼5 × 1046 erg into the stellar envelope ∼300 d pre-explosion, potentially from a nuclear burning instability at the onset of oxygen burning. Light-curve fitting, nebular spectroscopy, and pre-explosion HST imaging consistently point to a red supergiant (RSG) progenitor with $M_{\rm ZAMS}\approx 13.5\!-\!15 \, \mathrm{M}_{\odot }$, typical for SN II progenitor stars. This finding demonstrates that a typical RSG, like the progenitor of SN 2020fqv, has a complicated mass-loss history immediately before core collapse.

Tinyanont, Samaporn (ORCID:0000000214814676)↗
DOE OSTI2021

SDSS-IV DR17: final release of MaNGA PyMorph photometric and deep-learning morphological catalogues

ABSTRACT We present the MaNGA PyMorph photometric Value Added Catalogue (MPP-VAC-DR17) and the MaNGA Deep Learning Morphological VAC (MDLM-VAC-DR17) for the final data release of the MaNGA survey, which is part of the SDSS Data Release 17 (DR17). The MPP-VAC-DR17 provides photometric parameters from Sérsic and Sérsic+Exponential fits to the two-dimensional surface brightness profiles of the MaNGA DR17 galaxy sample in the g, r, and i bands (e.g. total fluxes, half-light radii, bulge-disc fractions, ellipticities, position angles, etc.). The MDLM-VAC-DR17 provides deep-learning-based morphological classifications for the same galaxies. The MDLM-VAC-DR17 includes a number of morphological properties, for example, a T-Type, a finer separation between elliptical and S0, as well as the identification of edge-on and barred galaxies. While the MPP-VAC-DR17 simply extends the MaNGA PyMorph photometric VAC published in the SDSS Data Release 15 (MPP-VAC-DR15) to now include galaxies that were added to make the final DR17, the MDLM-VAC-DR17 implements some changes and improvements compared to the previous release (MDLM-VAC-DR15): Namely, the low end of the T-Types is better recovered in this new version. The catalogue also includes a separation between early or late type, which classifies the two populations in a complementary way to the T-Type, especially at the intermediate types (−1 < T-Type < 2), where the T-Type values show a large scatter. In addition, k-fold-based uncertainties on the classifications are also provided. To ensure robustness and reliability, we have also visually inspected all the images. We describe the content of the catalogues and show some interesting ways in which they can be combined.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Pushing automated morphological classifications to their limits with the Dark Energy Survey

ABSTRACT We present morphological classifications of ∼27 million galaxies from the Dark Energy Survey (DES) Data Release 1 (DR1) using a supervised deep learning algorithm. The classification scheme separates: (a) early-type galaxies (ETGs) from late-type galaxies (LTGs); and (b) face-on galaxies from edge-on. Our convolutional neural networks (CNNs) are trained on a small subset of DES objects with previously known classifications. These typically have mr ≲ 17.7 mag; we model fainter objects to mr < 21.5 mag by simulating what the brighter objects with well-determined classifications would look like if they were at higher redshifts. The CNNs reach 97 per cent accuracy to mr < 21.5 on their training sets, suggesting that they are able to recover features more accurately than the human eye. We then used the trained CNNs to classify the vast majority of the other DES images. The final catalogue comprises five independent CNN predictions for each classification scheme, helping to determine if the CNN predictions are robust or not. We obtain secure classifications for ∼87 per cent and 73 per cent of the catalogue for the ETG versus LTG and edge-on versus face-on models, respectively. Combining the two classifications (a) and (b) helps to increase the purity of the ETG sample and to identify edge-on lenticular galaxies (as ETGs with high ellipticity). Where a comparison is possible, our classifications correlate very well with Sérsic index (n), ellipticity (ϵ), and spectral type, even for the fainter galaxies. This is the largest multiband catalogue of automated galaxy morphologies to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The influence of Hickson-like compact group environment on galaxy luminosities

Compact groups of galaxies are devised as extreme environments where interactions may drive galaxy evolution. In this work, we analysed whether the luminosities of galaxies inhabiting compact groups differ from those of galaxies in loose galaxy groups. We computed the luminosity functions of galaxy populations inhabiting a new sample of 1412 Hickson-like compact groups of galaxies identified in the Sloan Digital Sky Survey Data Release 16. We observed a characteristic absolute magnitude for galaxies in compact groups brighter than that observed in the field or loose galaxy systems. We also observed a deficiency of faint galaxies in compact groups in comparison with loose systems. Our analysis showed that the brightening is mainly due to galaxies inhabiting the more massive compact groups. In contrast to what is observed in loose systems, where only the luminosities of Red (and Early) galaxies show a dependency with group mass, luminosities of Red and Blue (also Early and Late) galaxies in compact groups are affected similarly as a function of group virial mass. When using Hubble types, we observed that elliptical galaxies in compact groups are the brightest galaxy population, and groups dominated by an elliptical galaxy also display the brightest luminosities in comparison with those dominated by spiral galaxies. Moreover, we show that the general luminosity trends can be reproduced using a mock catalogue obtained from a semi-analytical model of galaxy formation. These results suggest that the inner extreme environment in compact groups prompts a different evolutionary history for their galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

iMaNGA: mock MaNGA galaxies based on IllustrisTNG and MaStar SSPs – I. Construction and analysis of the mock data cubes

ABSTRACT Galaxy formation and evolution simulations are essential tools to probe poorly known astrophysics processes, but particular care is needed to compare simulations with galaxy observations, as observed data need to be modelled as well. We present a method to generate mock galaxies from the hydro-dynamical IllustrisTNG simulations which are suited to compare with integral field spectroscopic observation of galaxies from the SDSS-IV/Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. First, we include the same instrumental effects and procedures as adopted in the acquisition and analysis of real data. Furthermore, we generate the galaxy spectra from the simulations using new stellar population models based on the MaNGA stellar library (MaStar). In this way, our mock data cubes have the same spatial sampling, cover the same wavelength range (3600–10 300 Å), and share the same spectral resolution (R ≈ 1800) and flux calibration of real MaNGA galaxy spectra. In this first paper, we demonstrate the method over an early- and a late-type simulated galaxy from TNG50. We analyse the correspondent mock MaNGA-like data cubes with the same full spectral fitting code, firefly, which was used for the observed spectra. We find that the intrinsic and recovered age and metallicity gradients are consistent within 1σ, with residuals over all tassels consistent with 0 at the 68 ${{\ \rm per\ cent}}$ confidence level. We also perform the challenging test at comparing intrinsic and recovered star formation histories, finding a close resemblance between input and output. In follow-up papers, we will present a full simulated MaNGA-like catalogue (≈10 000 galaxies) with a comprehensive comparison of TNG50 simulations to MaNGA observational results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Why are we still using 3D masses for cluster cosmology?

ABSTRACT The abundance of clusters of galaxies is highly sensitive to the late-time evolution of the matter distribution, since clusters form at the highest density peaks. However, the 3D cluster mass cannot be inferred without deprojecting the observations, introducing model-dependent biases and uncertainties due to the mismatch between the assumed and the true cluster density profile and the neglected matter along the sightline. Since projected aperture masses can be measured directly in simulations and observationally through weak lensing, we argue that they are better suited for cluster cosmology. Using the Mira–Titan suite of gravity-only simulations, we show that aperture masses correlate strongly with 3D halo masses, albeit with large intrinsic scatter due to the varying matter distribution along the sightline. Nonetheless, aperture masses can be measured ≈2–3 times more precisely from observations, since they do not require assumptions about the density profile and are only affected by the shape noise in the weak lensing measurements. We emulate the cosmology dependence of the aperture mass function directly with a Gaussian process. Comparing the cosmology sensitivity of the aperture mass function and the 3D halo mass function for a fixed survey solid angle and redshift interval, we find the aperture mass sensitivity is higher for Ωm and $w_a$, similar for σ8, ns, and $w_0$, and slightly lower for h. With a carefully calibrated aperture mass function emulator, cluster cosmology analyses can use cluster aperture masses directly, reducing the sensitivity to model-dependent mass calibration biases and uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Sub-percentage measure of distances to redshift of 0.1 by a new cosmic ruler

ABSTRACT Distance-redshift diagrams probe expansion history of the Universe. We show that the stellar mass-binding energy (massE) relation of galaxies proposed in our previous study offers a new distance ruler at cosmic scales. By using elliptical galaxies in the main galaxy sample of the Sloan Digital Sky Survey Data Release 7, we construct a distance-redshift diagram over the redshift range from 0.05 to 0.2 with the massE ruler. The best-fit dark energy density is 0.675 ± 0.079 for flat Λ-cold dark matter (ΛCDM) model, consistent with those by other probes. At the median redshift of 0.11, the median distance is estimated to have a fractional error of 0.34 per cent, much lower than those by supernova (SN) Ia and baryonic acoustic oscillation (BAO) and even exceeding their future capability at this redshift. The above low-$\mathit{ z}$ measurement is useful for probing dark energy that dominates at the late Universe. For a flat dark energy equation of state model (flat wCDM), the massE alone constrains w to an error that is only a factor of 2.2, 1.7, and 1.3 times larger than those by BAO, SN Ia, and cosmic microwave background (CMB), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A detailed study of the barium central star of the planetary nebula Abell 70

ABSTRACT We present a detailed study of the barium star at the heart of the planetary nebula Abell 70. Time-series photometry obtained over a period of more than 10 yr demonstrates that the barium-contaminated companion is a rapid rotator with temporal variability due to spots. The amplitude and phasing of the photometric variability change abruptly; however, there is no evidence for a change in the rotation period (P = 2.06 d) over the course of the observations. The co-addition of 17 high-resolution spectra obtained with Ultraviolet and Visual Échelle Spectrograph mounted on the Very Large Telescope allows us to measure the physical and chemical properties of the companion, confirming it to be a chromospherically active, late G-type sub-giant with more than +1 dex of barium enhancement. We find no evidence of radial velocity variability in the spectra, obtained over the course of approximately 130 d with a single additional point some 8 yr later, with the radial velocities of all epochs approximately −10 km s −1 from the previously measured systemic velocity of the nebula. This is perhaps indicative that the binary has a relatively long period (P ≳ 2 yr) and high eccentricity (e ≳ 0.3), and that all the observations were taken around radial velocity minimum. However, unless the binary orbital plane is not aligned with the waist of the nebula or the systemic velocity of the binary is not equal to the literature value for the nebula, this would imply an unfeasibly large mass for the nebular progenitor.

Jones, David (ORCID:0000000339475946)↗
DOE OSTI2022

Chemical abundances of the Typhon stellar stream

ABSTRACT We present the first high-resolution chemical abundances of seven stars in the recently discovered high-energy stream Typhon. Typhon stars have apocentres r ≳ 100 kpc, making this the first detailed chemical picture of the Milky Way’s very distant stellar halo. Though the sample size is limited, we find that Typhon’s chemical abundances are more like a dwarf galaxy than a globular cluster, showing a metallicity dispersion and no presence of multiple stellar populations. Typhon stars display enhanced α-element abundances and increasing r-process abundances with increasing metallicity. The high-α abundances suggest a short star formation duration for Typhon, but this is at odds with expectations for the distant Milky Way halo and the presence of delayed r-process enrichment. If the progenitor of Typhon is indeed a new dwarf galaxy, possible scenarios explaining this apparent contradiction include a dynamical interaction that increases Typhon’s orbital energy, a burst of enhanced late-time star formation that raises [α/Fe], and/or group pre-processing by another dwarf galaxy before infall into the Milky Way. Alternatively, Typhon could be the high-energy tail of a more massive disrupted dwarf galaxy that lost energy through dynamical friction. We cannot clearly identify a known low-energy progenitor of Typhon in the Milky Way, but 70 per cent of high-apocentre stars in cosmological simulations are from high-energy tails of large dwarf galaxies. Typhon’s surprising combination of kinematics and chemistry thus underscores the need to fully characterize the dynamical history and detailed abundances of known substructures before identifying the origin of new substructures.

Astronomy & Astrophysics↗
DOE OSTI2022

Actinide opacities for modelling the spectra and light curves of kilonovae

ABSTRACT We extend previous ab initio calculations of lanthanide opacities to include a complete set of actinide opacities for use in the modelling of kilonova (KN) light curves and spectra. Detailed, fine-structure line features are generated using the configuration-interaction approach. These actinide opacities display similar trends to those observed for lanthanide opacities, such as the lighter actinides producing higher opacity than the heavier ones for relevant conditions in the dynamical ejecta. A line-binned treatment is employed to pre-compute opacity tables for 14 actinide elements (89 ≤ Z ≤ 102) over a grid of relevant temperatures and densities. These tabular opacities will be made publicly available for general usage in KN modelling. We demonstrate the usefulness of these opacities in KN simulations by exploring the sensitivity of light curves and spectra to different actinide abundance distributions that are predicted by different nuclear theories, as well as to different choices of ejecta mass and velocity. We find very little sensitivity to the two considered distributions, indicating that opacities for actinides with Z ≥ 99 do not contribute strongly. On the other hand, a single actinide element, protactinium, is found to produce faint spectral features in the far-infrared at late times (5–7 d post merger). More generally, we find that the choice of ejecta mass and velocity have the most significant effect on KN emission for this study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The whisper of a whimper of a bang: 2400 d of the Type Ia SN 2011fe reveals the decay of 55Fe

ABSTRACT We analyse new multifilter Hubble Space Telescope (HST) photometry of the normal Type Ia supernova (SN Ia) 2011fe out to ≈2400 d after maximum light, the latest observations to date of a SN Ia. We model the pseudo-bolometric light curve with a simple radioactive decay model and find energy input from both 57Co and 55Fe are needed to power the late-time luminosity. This is the first detection of 55Fe in a SN Ia. We consider potential sources of contamination such as a surviving companion star or delaying the deposition time-scale for 56Co positrons but these scenarios are ultimately disfavored. The relative isotopic abundances place direct constraints on the burning conditions experienced by the white dwarf (WD). Additionally, we place a conservative upper limit of <10−3 M⊙ on the synthesized mass of 44Ti. Only two classes of explosion models are currently consistent with all observations of SN 2011fe: (1) the delayed detonation of a low-ρc, near-MCh (1.2–1.3 M⊙) WD, or (2) a sub-MCh (1.0–1.1 M⊙) WD experiencing a thin-shell double detonation.

Astronomy & Astrophysics↗
DOE OSTI2022

Identification of Grand-design and Flocculent spirals from SDSS using deep convolutional neural network

ABSTRACT Spiral galaxies can be classified into the Grand-designs and Flocculents based on the nature of their spiral arms. The Grand-designs exhibit almost continuous and high contrast spiral arms and are believed to be driven by stationary density waves, while the Flocculents have patchy and low-contrast spiral features and are primarily stochastic in origin. We train a deep convolutional neural network model to classify spirals into Grand-designs and Flocculents, with a testing accuracy of $\mathrm{97.2{{\ \rm per\ cent}}}$. We then use the above model for classifying 1354 spirals from the SDSS. Out of these, 721 were identified as Flocculents, and the rest as Grand-designs. Interestingly, we find the mean asymptotic rotational velocities of our newly classified Grand-designs and Flocculents are 218 ± 86 and 146 ± 67 km s−1, respectively, indicating that the Grand-designs are mostly the high-mass and the Flocculents the intermediate-mass spirals. This is further corroborated by the observation that the mean morphological indices of the Grand-designs and Flocculents are 2.6 ± 1.8 and 4.7 ± 1.9, respectively, implying that the Flocculents primarily consist of a late-type galaxy population in contrast to the Grand-designs. Finally, an almost equal fraction of bars ∼0.3 in both the classes of spiral galaxies reveals that the presence of a bar component does not regulate the type of spiral arm hosted by a galaxy. Our results may have important implications for formation and evolution of spiral arms in galaxies.

Astronomy & Astrophysics↗
DOE OSTI2022

Lessons learned from the two largest Galaxy morphological classification catalogues built by convolutional neural networks

ABSTRACT We compare the two largest galaxy morphology catalogues, which separate early- and late-type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of ∼21 mag. The methodologies used for the construction of the catalogues include differences such as the cutout sizes, the labels used for training, and the input to the CNN – monochromatic images versus gri-band normalized images. In addition, one catalogue is trained using bright galaxies observed with DES (i < 18), while the other is trained with bright galaxies (r < 17.5) and ‘emulated’ galaxies up to r-band magnitude 22.5. Despite the different approaches, the agreement between the two catalogues is excellent up to i < 19, demonstrating that CNN predictions are reliable for samples at least one magnitude fainter than the training sample limit. It also shows that morphological classifications based on monochromatic images are comparable to those based on gri-band images, at least in the bright regime. At fainter magnitudes, i > 19, the overall agreement is good (∼95 per cent), but is mostly driven by the large spiral fraction in the two catalogues. In contrast, the agreement within the elliptical population is not as good, especially at faint magnitudes. By studying the mismatched cases, we are able to identify lenticular galaxies (at least up to i < 19), which are difficult to distinguish using standard classification approaches. The synergy of both catalogues provides an unique opportunity to select a population of unusual galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Colour gradients of low-redshift galaxies in the DESI Legacy Imaging Survey

ABSTRACT Radial colour gradients within galaxies arise from gradients of stellar age, metallicity, and dust reddening. Large samples of colour gradients from wide-area imaging surveys can complement smaller integral-field spectroscopy data sets and can be used to constrain galaxy formation models. Here, we measure colour gradients for low-redshift galaxies (z < 0.1) using photometry from the DESI Legacy Imaging Survey DR9. Our sample comprises ∼93 000 galaxies with spectroscopic redshifts and ∼574 000 galaxies with photometric redshifts. We focus on gradients across a radial range 0.5Reff to Reff, which corresponds to the inner disc of typical late-type systems at low redshift. This region has been the focus of previous statistical studies of colour gradients and has recently been explored by spectroscopic surveys such as MaNGA. We find that the colour gradients of most galaxies in our sample are negative (redder towards the centre), consistent with the literature. We investigate empirical relationships between colour gradient, average g − r and r − z colour, Mr, M⋆, and sSFR. Trends of gradient strength with Mr (M⋆) show an inflection around Mr ∼ −21 ($\log _{10} \, M_\star /\mathrm{M_\odot }\sim 10.5$). Below this mass, colour gradients become steeper with increasing M⋆, whereas colour gradients in more massive galaxies become shallower. We find that positive gradients (bluer stars at smaller radii) are typical for galaxies of $M_{\star }\sim 10^{8}\, \mathrm{M_\odot }$. We compare our results to age and metallicity gradients in two data sets derived from fits of different stellar population libraries to MaNGA spectra, but find no clear consensus explanation for the trends we observe. Both MaNGA data sets seem to imply a significant contribution from dust reddening, in particular, to explain the flatness of colour gradients along the red sequence.

Astronomy & Astrophysics↗
DOE OSTI2022

The X-ray view of optically selected dual AGN

We present a study of optically selected dual Active Galactic Nuclei (AGN) with projected separations of 3–97 kpc. Using multiwavelength (MWL) information (optical, X-ray, mid-IR), we characterized the intrinsic nuclear properties of this sample and compared them with those of isolated systems. Among the 124 X-ray-detected AGN candidates, 52 appear in pairs and 72 as single X-ray sources. Through MWL analysis, we confirmed the presence of the AGN in >80 per cent of the detected targets in pairs (42 out of 52). X-ray spectral analysis confirms the trend of increasing AGN luminosity with decreasing separation, suggesting that mergers may have contributed to triggering more luminous AGN. Through X-ray/mid-IR ratio versus X-ray colours, we estimated a fraction of Compton-thin AGN (with 10 22 cm −2 < N H < 10 24 cm −2 ) of about 80 per cent, while about 16 per cent are Compton-thick sources (with N H > 10 24 cm −2 ). These fractions of obscured sources are larger than those found in samples of isolated AGN, confirming that pairs of AGN show higher obscuration. This trend is further confirmed by comparing the de-reddened [O III ] emission with the observed X-ray luminosity. However, the derived fraction of Compton-thick sources in this sample at the early stages of merging is lower than that reported for late-merging dual-AGN samples. Comparing N H from X-rays with that derived from E(B − V ) from narrow-line regions, we found that the absorbing material is likely to be associated with the torus or broad-line regions. We also explored the X-ray detection efficiency of dual-AGN candidates, finding that, when observed properly (at on-axis positions and with long exposures), X-ray data represent a powerful way to confirm and investigate dual-AGN systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Accuracy of power spectra in dissipationless cosmological simulations

Abstract We exploit a suite of large N-body simulations (up to N = 40963) performed with abacus, of scale-free models with a range of spectral indices n, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At 1 per cent level it starts at early times from a fraction of $k_\Lambda$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, ${\sim}2{-}3 k_\Lambda$ at the very latest times we evolve to. At the $5{{\ \rm per\ cent}}$ level, accuracy extends up to wavenumbers of order $5k_\Lambda$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple n-dependence, allowing a extrapolation to place conservative bounds on the accuracy of N-body simulations of non-scale-free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small-scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_\Lambda$ may need revision in light of our convergence analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Long-term 3D MHD simulations of black hole accretion discs formed in neutron star mergers

ABSTRACT We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in flash 4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disc simulations that vary the initial magnetic field geometry and disc compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modelling the BH’s gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively inefficient state in the disc. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disc compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively inefficient state at late times. The electron fraction distribution of the disc outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding ${\sim} 50{{\ \rm per\ cent}}$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

Fahlman, Steven (ORCID:0000000348759940)↗