Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Massively parallel and universal approximation of nonlinear functions using diffractive processors

Nonlinear computation is essential for a wide range of information processing tasks, yet implementing nonlinear functions using optical systems remains a challenge due to the weak and power-intensive nature of optical nonlinearities. Overcoming this limitation without relying on nonlinear optical materials could unlock unprecedented opportunities for ultrafast and parallel optical computing systems. Here, we demonstrate that large-scale nonlinear computation can be performed using linear optics through optimized diffractive processors composed of passive phase-only surfaces. In this framework, the input variables of nonlinear functions are encoded into the phase of an optical wavefront—e.g., via a spatial light modulator (SLM)—and transformed by an optimized diffractive structure with spatially varying point-spread functions to yield output intensities that approximate a large set of unique nonlinear functions–all in parallel. We provide proof establishing that this architecture serves as a universal function approximator for an arbitrary set of bandlimited nonlinear functions, also covering wavelength-multiplexed nonlinear functions as well as multi-variate and complex-valued functions that are all-optically cascadable. Our analysis also indicates the successful approximation of typical nonlinear activation functions commonly used in neural networks, including the sigmoid, tanh, ReLU (rectified linear unit), and softplus. We numerically demonstrate the parallel computation of one million distinct nonlinear functions, accurately executed at wavelength-scale spatial density at the output of a diffractive optical processor. Furthermore, we experimentally validated this framework using in situ optical learning and approximated 35 unique nonlinear functions in a single shot using a compact setup consisting of an SLM and an image sensor. These results establish diffractive optical processors as a scalable platform for massively parallel universal nonlinear function approximation, paving the way for new capabilities in analog optical computing based on linear materials.

Rahman, Md Sadman Sakib [University of California,↗
DOE OSTI2022

A Massive, Position-Resolving, High-Energy-Resolution Detector for Non-Accelerator Cosmic and Intensity Frontier Particle Physics

We proposed to take the first steps in the development of a detector that promises energy resolution of tens of eV FWHM combined with robust nuclear-recoil discrimination and high fidelity, mm-precision position reconstruction for applications in non-accelerator particle physics at the Cosmic and Intensity Frontiers. The detector would obtain these excellent resolutions by sensing athermal phonons produced by particle interactions in crystalline, dielectric targets using a sensitive, highly multiplexable superconducting phonon sensor, the kinetic inductance detector (KID). This detector would be applicable to: the search for low-mass particle dark matter candidates with masses below 5~GeV via direct detection of scattering of dark matter particles with terrestrial nuclei; detection of coherent elastic neutrino-nucleus scattering to test for new physics such as a non-standard value of the weak nuclear charge, non-standard neutrino interactions (perhaps driven by a neutrino magnetic moment), or the existence of sterile neutrinos; and, searches for neutrinoless double-beta decay. During the funding period, we demonstrated scaling up of the detector concept from a 22-mm by 22-mm by 1-mm, 1-g prototype with 0.9~keV FWHM energy resolution to a 75-mm diameter by 1-mm, 9-g prototype while improving the inferred energy resolution to 0.7~keV~FWHM. In the process, we solved many problems associated with scaling device fabrication to large wafers and vastly reduced the fraction of the surface occupied by inactive but energy-absorbing metal. We also demonstrated a new technique that substantially simplifies the process of characterizing a new detector. These results provide a good foundation for future work scaling up the design to 4-mm thickness substrates and improving the resolution to reach 0.035~keV~FWHM, yielding a detector with compelling potential for the above applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Nematic Fluctuations in the Non-Superconducting Iron Pnictide BaFe 1.9-x Ni 0.1 Cr x As 2

The main driven force of the electronic nematic phase in iron-based superconductors is still under debate. Here, we report a comprehensive study on the nematic fluctuations in a non-superconducting iron pnictide system BaFe 1.9-x Ni 0.1 Cr x As 2 by electronic transport, angle-resolved photoemission spectroscopy (ARPES), and inelastic neutron scattering (INS) measurements. Previous neutron diffraction and transport measurements suggested that the collinear antiferromagnetism persists to x = 0.8, with similar Néel temperature TN and structural transition temperature Ts around 32 K, but the charge carriers change from electron type to hole type around x = 0.5. In this study, we have found that the in-plane resistivity anisotropy also highly depends on the Cr dopings and the type of charge carriers. While ARPES measurements suggest possibly weak orbital anisotropy onset near Ts for both x = 0.05 and x = 0.5 compounds, INS experiments reveal clearly different onset temperatures of low-energy spin excitation anisotropy, which is likely related to the energy scale of spin nematicity. These results suggest that the interplay between the local spins on Fe atoms and the itinerant electrons on Fermi surfaces is crucial to the nematic fluctuations of iron pnictides, where the orbital degree of freedom may behave differently from the spin degree of freedom, and the transport properties are intimately related to the spin dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Environmental sensitivities of shallow-cumulus dilution – Part 1: Selected thermodynamic conditions

Cumulus entrainment, and its consequent dilution of buoyant cloud cores, strongly regulates the life cycle of shallow cumuli yet remains poorly understood. Herein, new insights into this problem are obtained through large-eddy simulations that systematically investigate the sensitivity of shallow-cumulus dilution to cloud-layer relative humidity (RH), cloud- and subcloud-layer depths, and continentality (i.e., the land–ocean contrast). The simulated cloud-core dilution is found to be strongly sensitive to continentality, with fractional dilution rates twice as large over the ocean as over land. Using a similarity theory based on the turbulent-kinetic-energy (TKE) budget, the reduced cloud-core dilution over land is attributed to larger cloud-base mass flux (mb), driven by stronger surface heating and subcloud turbulence. As mb increases, the fractional dilution rate must decrease to maintain energetic equilibrium. A positive sensitivity is also found to cloud-layer RH, with the core dilution increasing by 25 %–50 % for a 10 % enhancement in RH. This sensitivity is interpreted using the buoyancy-sorting hypothesis, in that mixtures of cloud and environmental air are more likely to become negatively buoyant and detrain (rather than diluting the cloud core) in drier cloud layers. By contrast, the sensitivities of (marine) shallow-cumulus dilution to cloud- and subcloud-layer depths are weak, with a 3 % decrease for a doubling for the former and a 4 % reduction in dilution for a 50 % deeper subcloud layer. These surprisingly weak sensitivities are readily explained by offsetting effects in the TKE similarity theory. Altogether, these experimental findings provide useful, though still incomplete, guidance for flow-dependent shallow-cumulus entrainment formulations in large-scale models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

Radio and X-ray observations of the radio halo source in A1367

Radio and X-ray observations of the 'radio halo' in A1367 are presented. The radio spectral index distribution is derived combining data at 610 and 1415 MHz taken with the Westerbork Synthesis Radio Telescope (WSRT). A high-resolution Very Large Array (VLA) map allows us to discriminate between the extended cluster emission and a triple radio source associated with a background stellar object (most probably a QSO) which appears projected onto it. X-ray observations taken with the high-resolution imager on board the Einstein satellite reveal an X-ray source associated with this background object. A second weak X-ray enhancement is found projected onto the extended radio component. Inverse Compton scattering of the synchrotron-emitting particles provides a possible interpretation of such a source. The presence of large-scale magnetic fields and particle acceleration due to the turbulent motion of galaxies through the intergalactic medium is discussed in relationship to the observed peculiarities of the galaxies found in this region.

Gavazzi, G.↗
NASA NTRS1997

An Experimental Investigation of Steady and Unsteady Flow Field in an Axial Flow Turbine

Measurements were made in a large scale single stage turbine facility. Within the nozzle passage measurements were made using a five hole probe, a two-component Laser Doppler Velocimeter (LDV), and a single sensor hot wire probe. These measurements showed weak secondary flows at midchord, and two secondary flow loss cores at the nozzle exit. The casing vortex loss core was the larger of the two. At the exit radial inward flow was found over the entire passage, and was more pronounced in the wake. Nozzle wake decay was found to be more rapid than for an isolated vane row due to the rotor's presence. The midspan rotor flow field was measured using a two-component LDV. Measurements were made from upstream of the rotor to a chord behind the rotor. The distortion of the nozzle wake as it passed through the rotor blade row was determined. The unsteadiness in the rotor flow field was determined. The decay of the rotor wake was also characterized.

Zaccaria, M.↗
NASA NTRS2007

Pixelation Effects in Weak Lensing

Weak gravitational lensing can be used to investigate both dark matter and dark energy but requires accurate measurements of the shapes of faint, distant galaxies. Such measurements are hindered by the finite resolution and pixel scale of digital cameras. We investigate the optimum choice of pixel scale for a space-based mission, using the engineering model and survey strategy of the proposed Supernova Acceleration Probe as a baseline. We do this by simulating realistic astronomical images containing a known input shear signal and then attempting to recover the signal using the Rhodes, Refregier, and Groth algorithm. We find that the quality of shear measurement is always improved by smaller pixels. However, in practice, telescopes are usually limited to a finite number of pixels and operational life span, so the total area of a survey increases with pixel size. We therefore fix the survey lifetime and the number of pixels in the focal plane while varying the pixel scale, thereby effectively varying the survey size. In a pure trade-off for image resolution versus survey area, we find that measurements of the matter power spectrum would have minimum statistical error with a pixel scale of 0.09' for a 0.14' FWHM point-spread function (PSF). The pixel scale could be increased to ~0.16' if images dithered by exactly half-pixel offsets were always available. Some of our results do depend on our adopted shape measurement method and should be regarded as an upper limit: future pipelines may require smaller pixels to overcome systematic floors not yet accessible, and, in certain circumstances, measuring the shape of the PSF might be more difficult than those of galaxies. However, the relative trends in our analysis are robust, especially those of the surface density of resolved galaxies. Our approach thus provides a snapshot of potential in available technology, and a practical counterpart to analytic studies of pixelation, which necessarily assume an idealized shape measurement method.

gravitational lensing↗
NASA NTRS2020

The Prevalence of Precipitation from Polar Supercooled Clouds

Supercooled clouds substantially impact polar surface energy budgets but large-scale models often underestimate their occurrence, which motivates accurately establishing metrics of basic processes. An analysis of long-term measurements at Utqiaġvik, Alaska, and McMurdo Station, Antarctica, combines lidar-validated use of soundings to identify supercooled cloud layers and colocated ground-based profiling radar measurements to quantify cloud base precipitation. We find that more than 85 % (75 %) of sampled supercooled layers are precipitating over the Arctic (Antarctic) site, with more than 75 % (50 %) precipitating continuously to the surface. Such high frequencies can be reconciled with substantially lesser spaceborne estimates by considering differences in radar hydrometeor detection sensitivity. While ice precipitation into supercooled clouds from aloft is common, we also find that the great majority of supercooled cloud layers without ice falling into them are themselves continuously generating precipitation. Such sustained primary ice formation is consistent with continuous activation of immersion-mode ice nucleating particles (INPs), suggesting that supercooled cloud formation is a principal gateway to ice formation at temperatures greater than ~ −38 °C over polar regions. The prevalence of weak precipitation fluxes is also consistent with supercooled cloud longevity, and with well-observed and widely simulated case studies. An analysis of colocated microwave radiometer retrievals suggests that weak precipitation fluxes can be nonetheless consequential to moisture budgets for supercooled clouds owing to small liquid water paths. Finally, we suggest that these ground-based precipitation rate statistics offer valuable guidance for improving the representation of polar cloud processes in large-scale models.

Polar clouds↗
NASA NTRS1982

Numerical prediction of the Mid-Atlantic states cyclone of 18-19 February 1979

A series of forecast experiments was conducted to assess the accuracy of the GLAS model, and to determine the importance of large scale dynamical processes and diabatic heating to the cyclogenesis. The GLAS model correctly predicted intense coastal cyclogenesis and heavy precipitation. Repeated without surface heat and moisture fluxes, the model failed to predict any cyclone development. An extended range forecast, a forecast from the NMC analysis interpolated to the GLAS grid, and a forecast from the GLAS analysis with the surface moisture flux excluded predicted weak coastal low development. Diabatic heating resulting from oceanic fluxes significantly contributed to the generation of low level cyclonic vorticity and the intensification and slow rate of movement of an upper level ridge over the western Atlantic. As an upper level short wave trough approached this ridge, diabatic heating associated with the release of latent heat intensified, and the gradient of vorticity, vorticity advection and upper level divergence in advance of the trough were greatly increased, providing strong large scale forcing for the surface cyclogenesis.

Atlas, R.↗
DOE OSTI2023

Detection of the significant impact of source clustering on higher order statistics with DES Year 3 weak gravitational lensing data

We measure the impact of source galaxy clustering on higher order summary statistics of weak gravitational lensing data. By comparing simulated data with galaxies that either trace or do not trace the underlying density field, we show that this effect can exceed measurement uncertainties for common higher order statistics for certain analysis choices. We evaluate the impact on different weak lensing observables, finding that third moments and wavelet phase harmonics are more affected than peak count statistics. Using Dark Energy Survey (DES) Year 3 (Y3) data, we construct null tests for the source-clustering-free case, finding a p-value of p = 4 × 10 −3 (2.6σ) using third-order map moments and p = 3 × 10 −11 (6.5σ) using wavelet phase harmonics. The impact of source clustering on cosmological inference can be either included in the model or minimized through ad hoc procedures (e.g. scale cuts). We verify that the procedures adopted in existing DES Y3 cosmological analyses were sufficient to render this effect negligible. Failing to account for source clustering can significantly impact cosmological inference from higher order gravitational lensing statistics, e.g. higher order N-point functions, wavelet-moment observables, and deep learning or field-level summary statistics of weak lensing maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Redshift inference from the combination of galaxy colours and clustering in a hierarchical Bayesian model – Application to realistic N -body simulations

ABSTRACT Photometric galaxy surveys constitute a powerful cosmological probe but rely on the accurate characterization of their redshift distributions using only broad-band imaging, and can be very sensitive to incomplete or biased priors used for redshift calibration. A hierarchical Bayesian model has recently been developed to estimate those from the robust combination of prior information, photometry of single galaxies, and the information contained in the galaxy clustering against a well-characterized tracer population. In this work, we extend the method so that it can be applied to real data, developing some necessary new extensions to it, especially in the treatment of galaxy clustering information, and we test it on realistic simulations. After marginalizing over the mapping between the clustering estimator and the actual density distribution of the sample galaxies, and using prior information from a small patch of the survey, we find the incorporation of clustering information with photo-z’s tightens the redshift posteriors and overcomes biases in the prior that mimic those happening in spectroscopic samples. The method presented here uses all the information at hand to reduce prior biases and incompleteness. Even in cases where we artificially bias the spectroscopic sample to induce a shift in mean redshift of $\Delta \bar{z} \approx 0.05,$ the final biases in the posterior are $\Delta \bar{z} \lesssim 0.003.$ This robustness to flaws in the redshift prior or training samples would constitute a milestone for the control of redshift systematic uncertainties in future weak lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Relaxation of the degenerate one-dimensional Fermi gas

In this work, we study how a system of one-dimensional spin-1/2 fermions at temperatures well below the Fermi energy approaches thermal equilibrium. The interactions between fermions are assumed to be weak and are accounted for within the perturbation theory. In the absence of an external magnetic field, spin degeneracy strongly affects relaxation of the Fermi gas. For sufficiently short-range interactions, the rate of relaxation scales linearly with temperature. Focusing on the case of the system near equilibrium, we linearize the collision integral and find the exact solution of the resulting relaxation problem. Furthermore, we discuss the application of our results to the evaluation of the transport coefficients of the one-dimensional Fermi gas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Interaction-driven quantum anomalous Hall insulator in a Dirac semimetal

Here, the interaction-driven quantum anomalous Hall (QAH) insulator has been sought for a long time in a Dirac semimetal with linear band touching points at the Fermi level. By combining exact diagonalization, density matrix renormalization group, and analytical methods, we study a spinless fermion system on the checkerboard lattice with twofold rotational symmetry, which realizes two Dirac band touching points in the absence of interaction. At weak coupling, the Dirac semimetal is stable. At a finite density-density repulsive interaction, we analyze possible symmetry broken states, and find that an QAH state is stabilized when the interaction strength exceeds the energy scale controlling the separation between the Dirac points. Through numerical simulations, we verify the existence of the QAH phase with spontaneous time-reversal symmetry breaking and quantized Chern number C = 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An acoustic investigation of the near-surface turbulence on Mars

The Perseverance rover is carrying out an original acoustic experiment on Mars: the SuperCam microphone records the spherical acoustic waves generated by laser sparks at distances from 2 m to more than 8 m. These N-shaped acoustic waves scatter from the multiple local heterogeneities of the turbulent atmosphere. Therefore, large and random fluctuations of sound travel time and intensity develop as the waves cross the medium. The variances of the travel times and the scintillation index (normalized variance of the sound intensity) are studied within the mathematical formalism of the propagation of spherical acoustic waves through thermal turbulence to infer statistical properties of the Mars atmospheric temperature fluctuation field. The comparison with the theory is made by simplifying assumptions that do not include wind fluctuations and diffraction effects. Two Earth years (about one Martian year) of observations acquired during the maximum convective period (10:00–14:00 Mars local time) show a good agreement between the dataset and the formalism: the travel time variance diverges from the linear Chernov solution exactly where the density of occurrence of the first caustic reaches its maximum. Moreover, on average, waves travel faster than the mean speed of sound due to a fast path effect, which is also observed on Earth. To account for the distribution of turbulent eddies, several power spectra are tested and the best match to observation is obtained with a generalized von Karman spectrum with a shallower slope than the Kolmogorov cascade, φ(k)∝(1 + k 2 L 2 ) -4/3 . It is associated with an outer scale of turbulence, L, of 11 cm at 2 m above the surface and a standard deviation of 6 K over 9 s for the temperature. These near-surface atmospheric properties are consistent with a weak to moderate wave scattering regime around noon with little saturation. Overall, this study presents an innovative and promising methodology to probe the near-surface atmospheric turbulence on Mars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Host Dark Matter Halos of SDSS Red and Blue Quasars: No Significant Difference in Large-scale Environment

The observed optical colors of quasars are generally interpreted in one of two frameworks: unified models that attribute the color to the random orientation of the accretion disk along the line of sight, and evolutionary models that invoke connections between quasar systems and their environments. We test these schemas by probing the dark matter halo environments of optically selected quasars as a function of g – i optical color by measuring the two-point correlation functions of ~0.34 million eBOSS quasars as well as the gravitational deflection of cosmic microwave background photons around ~0.66 million XDQSO photometric quasar candidates. We do not detect a trend of halo bias with optical color through either analysis, finding that optically selected quasars at 0.8 < z < 2.2 occupy halos of characteristic mass M h ~ 3 × 10 12 h –1 M ⊙ regardless of their color. This result implies that a quasar's large-scale halo environment is not strongly connected to its observed optical color. We also confirm the findings of fundamental differences in the radio properties of red and blue quasars by stacking 1.4 GHz FIRST images at their positions, suggesting the observed differences cannot be attributed to orientation. Instead, the differences between red and blue quasars likely arise on nuclear-galactic scales, perhaps owing to reddening by a nuclear dusty wind. Finally, we show that optically selected quasars' halo environments are also independent of their r – W2 optical–infrared colors, while previous work has suggested that mid-infrared-selected obscured quasars occupy more massive halos. We discuss the implications of this result for models of quasar and galaxy coevolution.

79 ASTRONOMY AND ASTROPHYSICS↗