Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cents”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Analysing the effects of heating and gas puffing in Proto-MPEX helicon and auxiliary heated plasmas

The Material Plasma Exposure Experiment (MPEX) is being constructed at Oak Ridge National Laboratory to investigate critical fusion reactor issues, such as plasma–material interactions (PMIs) under reactor-relevant conditions and time scales. The linear device Proto-MPEX was used as a test bed to address anticipated research and development issues associated with heating scenarios and establish the physics basis for MPEX. The SOLPS-ITER code suite has been applied to understand plasma and neutral transport in Proto-MPEX and to increase confidence in predictive simulations for MPEX. Coupling between COMSOL and SOLPS is performed to implement a 2D electron heating profile of the helicon source. The simulations show reasonable agreement with the experimental data for plasma with helicon and auxiliary electron cyclotron heating (ECH). Both Bohm and constant diffusion (D$_\perp$: 0.5 m 2 s –1 and $\chi_\perp$: 1 m 2 s –1 ) simulations show similar levels of agreement with respect to the sparse experimental data available, assuming a few per cent impurity concentration. ECH significantly increases the target electron temperature, however, the target electron density is reduced compared to helicon-only heated plasmas due to an increase in flow velocity and radial losses. The simulations show that further increasing the ECH power results in an increase in the target electron density due to increased recycling flux and ionization. The results indicate that ECH significantly enhances the target heat flux, with ECH power of 50 kW increasing the target heat flux from 0.4 to 17 MW m –2 . It is found that a small amount of gas puffing (GP) near the target plate can further increase the target heat fluxes at the higher ECH power cases, but the target heat flux is reduced at higher GP conditions due to a significant reduction in electron temperature via radiation. ECH and GP scenarios can generate a higher target flux, facilitating improved PMI studies with more reactor-relevant plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A comparison of template vs. direct model fitting for redshift-space distortions in BOSS

The growth of large-scale structure, as revealed in the anisotropic of clustering of galaxies in the low redshift Universe at z < 2, provides a stringent test of our cosmological model. The strongest current constraints come from the BOSS and eBOSS surveys, with uncertainties on σ 8 , the amplitude of clustering on an 8 h –1 Mpc scale, of less than 10 per cent. A number of different approaches have been taken to fitting this signal, leading to discrepancies of up to 1σ in the measurements of the amplitude of fluctuations at late times. We compare in some detail two of the leading approaches, one based on fitting a template cosmology whose amplitude and length scales are allowed to float with one based on varying the underlying parameters of a cosmological model directly, when fitting to the BOSS DR12 data. Holding the input data, scale cuts, window functions and modeling framework fixed we are able to isolate the cause of the differences and discuss the implications for future surveys.

redshift surveys↗
DOE OSTI2020

Internal structure of the San Jacinto fault zone at the Ramona Reservation, north of Anza, California, from dense array seismic data

SUMMARY We image the internal structure of the San Jacinto fault zone (SJFZ) near Anza, California, with seismic data recorded by two dense arrays (RA and RR) from ∼42 000 local and ∼180 teleseismic events occurring between 2012 and 2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analysed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km × 1.4 km area) that operated for about a month in September–October 2016. The SJFZ at the site contains three near-parallel surface traces F1, F2 and F3 from SW to NE that have accommodated several Mw > 6 earthquakes in the past 15 000 yr. Waveform changes in the fault normal direction indicate structural discontinuities that are consistent with the three fault surface traces. Relative slowness from local events and delay time analysis of teleseismic arrivals in the fault normal direction suggest a slower SW side than the NE with a core damage zone between F1 and F2. This core damage zone causes ∼0.05 s delay at stations RR26–31 in the teleseismic P arrivals compared with the SW-most station, and generates both P- and S-type fault zone trapped waves. Inversion of S trapped waves indicates the core damaged structure is ∼100 m wide, ∼4 km deep with a Q value of ∼20 and 40 per cent S-wave velocity reduction compared with bounding rocks. Fault zone head waves observed at stations SW of F3 indicate a local bimaterial interface that separates the locally faster NE block from the broad damage zone in the SW at shallow depth and merges with a deep interface that separates the regionally faster NE block from rocks to the SW with slower velocities at greater depth. The multiscale structural components observed at the site are related to the geological and earthquake rupture history at the site, and provide important information on the preferred NW propagation of earthquake ruptures on the San Jacinto fault.

Geochemistry & Geophysics↗
DOE OSTI2020

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗
DOE OSTI2021

Characterization of damage processes in Montney siltstone under triaxial compression using acoustic emission and diagnostic imaging

SUMMARY Crack nucleation and rock failure processes in a fine-grained siltstone (Montney Formation) under triaxial compression are investigated using combined diagnostic techniques, including ultrasonic-wave measurement, acoustic-emission (AE) monitoring, computed tomography (CT) scanning, and thin-section imaging. The sample displays a weak-to-moderate inherent seismic anisotropy and noticeable stress-induced anisotropy prior to failure. No AE event was detected until the applied axial stress reached 95 per cent of the peak value. The signal-to-noise ratio is relatively low, however, and detectable AE events are more diffuse than those observed in highly brittle rocks. The AE locations correlate with a shear fracture zone imaged by CT scanning. AE moment-tensor analysis reveals that events with larger relative magnitudes are characterized by high volumetric (tensile or compressive) components, and the initiation of the failure zone is dominated by combined shear-tensile (ST) failure. Stress inversion of the AE events with high tensile components is in good agreement with the known applied stress. Microscopic imaging of thin sections from the failed sample shows that the failure zone is an en echelon structure consisting of a major fracture with branching micro and minor cracks. This failure mechanism is consistent with a ST source mechanism and is interpreted to be associated with the fine granular structure and mineral composition of Montney siltstone.

Jia, Suzie Q. (ORCID:0000000326496190)↗
DOE OSTI2022

Limited Dynamic Triggering in the Utah Region, USA

The state of Utah, USA, experiences around 3800 catalogued earthquakes per year, highlighting that the region is seismically active and susceptible to earthquakes. Following the 2002 Denali Fault (M7.9) earthquake in Alaska, the region showed an elevated seismicity rate for 3 weeks following the passage of high amplitude surface waves, suggesting that the region may be particularly susceptible to dynamic triggering. With over 23,396 faults and each fault presenting a potential fault for triggering, we systematically search for dynamic triggering throughout the state of Utah caused by large, global earthquakes with M ≥ 7. Specifically, we analyse earthquake catalogues and all available waveform data to determine statistical increases of seismicity rate following the passage of seismic arrivals. While we find instances of dynamic triggering, our results show that these events occur sparsely in the region. In total, less than 20 per cent of the 273 main shocks that occur from 2000 to the end of 2017 show a statistical indication of dynamic triggering throughout the Utah region, highlighting that dynamic triggering is limited for stresses created by transient signals from global M ≥ 7 earthquakes, with the exception being the Denali Fault (M7.9), Alaska earthquake (i.e. an instance of significant triggering).

58 GEOSCIENCES↗
DOE OSTI2022

Source scaling comparison and validation in Central Italy: data intensive direct S waves versus the sparse data coda envelope methodology

SUMMARY Robustness of source parameter estimates is a fundamental issue in understanding the relationships between small and large events; however, it is difficult to assess how much of the variability of the source parameters can be attributed to the physical source characteristics or to the uncertainties of the methods and data used to estimate the values. In this study, we apply the coda method by Mayeda et al. using the coda calibration tool (CCT), a freely available Java-based code (https://github.com/LLNL/coda-calibration-tool) to obtain a regional calibration for Central Italy for estimating stable source parameters. We demonstrate the power of the coda technique in this region and show that it provides the same robustness in source parameter estimation as a data-driven methodology [generalized inversion technique (GIT)], but with much fewer calibration events and stations. The Central Italy region is ideal for both GIT and coda approaches as it is characterized by high-quality data, including recent well-recorded seismic sequences such as L'Aquila (2009) and Amatrice–Norcia–Visso (2016–2017). This allows us to apply data-driven methods such as GIT and coda-based methods that require few, but high-quality data. The data set for GIT analysis includes ∼5000 earthquakes and more than 600 stations, while for coda analysis we used a small subset of 39 events spanning 3.5 < Mw < 6.33 and 14 well-distributed broad-band stations. For the common calibration events, as well as an additional 247 events (∼1.7 < Mw < ∼5.0) not used in either calibration, we find excellent agreement between GIT-derived and CCT-derived source spectra. This confirms the ability of the coda approach to obtain stable source parameters even with few calibration events and stations. Even reducing the coda calibration data set by 75 per cent, we found no appreciable degradation in performance. This validation of the coda calibration approach over a broad range of event size demonstrates that this procedure, once extended to other regions, represents a powerful tool for future routine applications to homogeneously evaluate robust source parameters on a national scale. Furthermore, the coda calibration procedure can homogenize the Mw estimates for small and large events without the necessity of introducing any conversion scale between narrow-band measures such as local magnitude (ML) and Mw, which has been shown to introduce significant bias.

Morasca, Paola (ORCID:0000000265254867)↗
DOE OSTI2022

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗
DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗
DOE OSTI2024

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics↗
DOE OSTI2025

Data-informed grid refinement to improve traveltime accuracy in the regional seismic traveltime (RSTT) model

The regional seismic traveltime (RSTT) model predicts traveltimes of regional seismic phases accounting for 3-D structure of the crust and the upper mantle on a global scale. Previous versions of the RSTT model have been implemented using nodes separated by ∼1° spacing across the globe. A regional-scale study using regional Pn and Pg traveltimes across Israel and the Middle East demonstrated that data driven, systematic grid refinement reduces traveltime residuals and enhances resolution of smaller tectonic features in regions having dense ray coverage. High density Pn ray coverage in the western US, Europe, Middle East and East Asia can likewise provide the resolution that allows systematic global grid refinement of the RSTT model. In this study, we use a large number of Pn ray paths originating from events located with an epicentral location uncertainty of 25 km (GT25) or better. We conduct targeted grid refinements at 1.0°, 0.5°, 0.25° and 0.125° on a global scale, producing a refined RSTT model that yields a 21.6 per cent reduction in median event location error in Europe and the Middle East, when compared with the original global RSTT model presented in Begnaud et al. The new model also resolves finer tectonic structures in regions with high Pn ray density.

58 GEOSCIENCES↗
DOE OSTI2020

Reconstruction with velocities

ABSTRACT Reconstruction is becoming a crucial procedure of galaxy clustering analysis for future spectroscopic redshift surveys to obtain subper cent level measurement of the baryon acoustic oscillation scale. Most reconstruction algorithms rely on an estimation of the displacement field from the observed galaxy distribution. However, the displacement reconstruction degrades near the survey boundary due to incomplete data and the boundary effects extend to ${\sim}100\, \mathrm{Mpc}/h$ within the interior of the survey volume. We study the possibility of using radial velocities measured from the cosmic microwave background observation through the kinematic Sunyaev–Zeldovich effect to improve performance near the boundary. We find that the boundary effect can be reduced to ${\sim}30-40\, \mathrm{Mpc}/h$ with the velocity information from Simons Observatory. This is especially helpful for dense low redshift surveys where the volume is relatively small and a large fraction of total volume is affected by the boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Galaxy mergers in eagle do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN

ABSTRACT We investigate the connection between galaxy–galaxy mergers and enhanced black hole (BH) growth using the cosmological hydrodynamical eagle simulation. We do this via three methods of analysis, investigating: the merger fraction of AGN, the AGN fraction of merging systems, and the AGN fraction of galaxies with close companions. In each case, we find an increased abundance of AGN within merging systems relative to control samples of inactive or isolated galaxies (by up to a factor of ≈3 depending on the analysis method used), confirming that mergers are enhancing BH accretion rates for at least a subset of the galaxy population. The greatest excess of AGN triggered via a merger are found in lower mass (M* ∼ 1010 M⊙) gas rich (fgas > 0.2) central galaxies with lower mass BHs (MBH ∼ 107 M⊙) at lower redshifts (z < 1). We find no enhancement of AGN triggered via mergers in more massive galaxies (M* ≳ 1011 M⊙). The enhancement of AGN is not uniform throughout the phases of a merger, and instead peaks within the early remnants of merging systems (typically lagging ≈300 Myr post-coalescence of the two galaxies at z = 0.5). We argue that neither major ($M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \ge \frac{1}{4}$) nor minor mergers ($\frac{1}{10} \lt M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \lt \frac{1}{4}$) are statistically relevant for enhancing BH masses globally. Whilst at all redshifts the galaxies experiencing a merger have accretion rates that are on average 2–3 times that of isolated galaxies, the majority of mass that is accreted on to BHs occurs outside the periods of a merger. We compute that on average no more than 15 per cent of a BHs final day mass comes from the enhanced accretion rates triggered via a merger.

McAlpine, Stuart↗
DOE OSTI2020

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first ~1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass haloes. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant, and scaling with stellar mass) from both star formation ( |$\langle f_{\rm esc}^{\rm sf} \rangle$| ) and AGN ( |$f_{\rm esc}^{\rm bh}$| ) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ( |$M_*\lt 10^9 \, {\rm \rm M_\odot }$| ) galaxies at z > 6 with the AGN contribution (driven by |$M_{bh}\gt 10^6 \, {\rm \rm M_\odot }$| black holes in |$M_* \gtrsim 10^9\, {\rm \rm M_\odot}$| galaxies) dominating at lower redshifts; (ii) AGN only contribute |$10-25{{\ \rm per\ cent}}$| to the cumulative ionizing emissivity by z = 4 for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of |$\langle f_{\rm esc}^{\rm sf} \rangle$| is shallower than |$f_{\rm esc}^{\rm bh}$| , at z < 7 a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at z ~ 6.8, they take-over at stellar masses that are a tenth of the knee mass by z = 4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The predicted properties of helium-enriched globular cluster progenitors at high redshift

Globular cluster progenitors may have been detected by Hubble Space Telescope, and are predicted to be observable with James Webb Space Telescope (JWST) and ground-based extremely large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched (ΔY = +0.12) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first ~30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being ~0.40 mag brighter at $\lambda \approx 2.0\, {\mu \rm m}$, and to be 0.30-mag redder in the JWST–NIRCam colour (F070W - F200W), and to actually be fainter for $\lambda \lesssim 0.50 \, {\mu \rm m}$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim \!5{{\ \rm per\ cent}}$, though we show that the Lyman-α escape fraction could end up higher for helium-enriched stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The SDSS/APOGEE catalogue of HgMn stars

ABSTRACT We report on H-band spectra of chemically peculiar Mercury–Manganese (HgMn) stars obtained via the SDSS/APOGEE survey. As opposed to other varieties of chemically peculiar stars such as classical Ap/Bp stars, HgMn stars lack strong magnetic fields and are defined by extreme overabundances of Mn, Hg, and other heavy elements. A satisfactory explanation for the abundance patterns remains to be determined, but low rotational velocity is a requirement and involvement in binary/multiple systems may be as well. The APOGEE HgMn sample currently consists of 269 stars that were identified among the telluric standard stars as those whose metallic absorption content is limited to or dominated by the H-band Mn ii lines. Due to the fainter magnitudes probed by the APOGEE survey as compared to past studies, only 9/269 stars in the sample were previously known as HgMn types. The 260 newly identified HgMn stars represents a more than doubling of the known sample. At least 32 per cent of the APOGEE sample are found to be binary or multiple systems, and from multi-epoch spectroscopy, we were able to determine orbital solutions for at least one component in 32 binaries. Many of the multilined systems include chemically peculiar companions, with noteworthy examples being the HgMn+Ap/Bp binary HD 5429, the HgMn+HgMn binary HD 298641, and the HgMn+Bp Mn + Am triple system HD 231263. As a further peculiarity, roughly half of the sample produces narrow emission in the C i 16895 Å line, with widths and radial velocities that match those of the Mn ii lines.

Chojnowski, S. Drew↗
DOE OSTI2020

Multiplicity functions of quasars: predictions from the MassiveBlackII simulation

ABSTRACT We examine multiple active galactic nucleus (AGN) systems (triples and quadruples, in particular) in the MassiveBlackII simulation over a redshift range of 0.06 ≲ z ≲ 4. We identify AGN systems (with bolometric luminosity $L_{\mathrm{bol}}\gt 10^{42}~\mathrm{erg\, s}^{-1}$) at different scales (defined by the maximum distance between member AGNs) to determine the AGN multiplicity functions. This is defined as the volume/surface density of AGN systems per unit richness R, the number of AGNs in a system. We find that gravitationally bound multiple AGN systems tend to populate scales of ${\lesssim}0.7~\mathrm{cMpc}\, h^{-1}$; this corresponds to angular separations of ≲100 arcsec and a line-of-sight velocity difference ${\lesssim}200~\mathrm{km\, s}^{-1}$. The simulation contains ∼10 and ∼100 triples/quadruples per deg2 up to depths of DESI (g ≲ 24) and LSST (g ≲ 26) imaging, respectively; at least $20{{\ \rm per\ cent}}$ of these should be detectable in spectroscopic surveys. The simulated quasar ($L_{\mathrm{bol}}\gt 10^{44}~\mathrm{erg\, s}^{-1}$) triples and quadruples predominantly exist at 1.5 ≲ z ≲ 3. Their members have black hole masses $10^{6.5}\lesssim M_{\mathrm{ bh}}\lesssim 10^{9}~\mathrm{M}_{\odot }\, h^{-1}$ and live in separate (one central and multiple satellite) galaxies with stellar masses $10^{10}\lesssim M_{*}\lesssim 10^{12}~\mathrm{M}_{\odot }\, h^{-1}$. They live in the most massive haloes (e.g. ${\sim}10^{13}~\mathrm{M}_{\odot }\, h^{-1}$ at z = 2.5; ${\sim}10^{14}~\mathrm{M}_{\odot }\, h^{-1}$ at z = 1) in the simulation. Their detections provide an exciting prospect for understanding massive black hole growth and their merger rates in galaxies in the era of multimessenger astronomy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A large, deep 3 deg2 survey of H α, [O iii ], and [O ii ] emitters from LAGER: constraining luminosity functions

ABSTRACT We present our measurements of the H α, [O iii], and [O ii] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 z = 0.47 H α-, 3933 z = 0.93 [O iii]-, and 5367 z = 1.59 [O ii]-selected emission line galaxies in a 3 deg2 single, CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5σ depths of 8.2 × 10−18 erg s−1 cm−2 and comoving volumes of (1−7) × 105 Mpc3 making our survey one of the deepest narrow-band surveys. We select our emission line galaxies via spectroscopic confirmation, photometric redshifts, and colour–colour selections. We measure the observed luminosity functions for each sample and find best fits of $\phi ^\star = 10^{-3.16^{+0.09}_{-0.09}}$ Mpc−3 and $L^\star = 10^{41.72^{+0.09}_{-0.09}}$ erg s−1 for H α, $\phi ^\star = 10^{-2.16^{+0.10}_{-0.12}}$ Mpc−3 and $L^\star = 10^{41.38^{+0.07}_{-0.06}}$ erg s−1 for [O iii], and $\phi ^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc−3 and $L^\star = 10^{41.66^{+0.03}_{-0.03}}$ erg s−1 for [O ii], with α fixed to −1.75, −1.6, and −1.3, respectively. An excess of bright >1042 erg s−1 [O iii] emitters is observed and may be due to active galactic nucleus (AGN) contamination. Corrections for dust attenuation are applied assuming AHα = 1 mag. We also design our own empirical rest-frame g − r calibration using SDSS DR12 data, test it against our z = 0.47 H α emitters with zCOSMOS 1D spectra, and calibrate it for (g − r) between −0.8 and 1.3 mag. Dust and AGN-corrected star formation rate densities (SFRDs) are measured as log10ρSFR/(M⊙ yr−1 Mpc−3) = −1.63 ± 0.04, −1.07 ± 0.06, and −0.90 ± 0.10 for H α, [O iii], and [O ii], respectively. We find our [O iii] and [O ii] samples fully trace cosmic star formation activity at their respective redshifts in comparison to multiwavelength SFRDs, while the H α sample traces ∼70 per cent of the total z = 0.47 SFRD.

Khostovan, A. A.↗