Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constants”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Self-diffusion of liquid deuterium hydride and liquid tritium

Here, we present a quasi-elastic neutron scattering study of liquid deuterium hydride carried out using the Disk Chopper Spectrometer at the National Institute of Standards and Technology. Under saturated vapor pressure, the self-diffusion constant of deuterium hydride obeys an Arrhenius law D = D 0 exp(-E A /k B T), where the prefactor D 0 is given by D 0 = 9.5 ± 1.2 Å 2 /ps and the activation energy is given by E A = 58 ± 2 K. We apply the quantum law of corresponding states to the known diffusion constants of the hydrogen isotopologues. From this application, we estimate that D 0 ≈ 9.1 Å 2 /ps and E A ≈ 75 K in liquid tritium. Young’s theory of quantum-mechanical effects in van der Waals fluids is shown to apply to the diffusion constants of the liquid hydrogens. Our results underscore the importance of nuclear quantum effects in shaping the properties and behavior of the hydrogen isotopologues.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗
DOE OSTI2026

Synthesis and composition control of epitaxial FeWO4 thin films using pulsed laser deposition

We report the growth of epitaxial FeWO4 thin films on c-plane sapphire via pulsed laser deposition using a Fe 2 O 3 /WO 3 target ablated in an O 2 atmosphere. At a constant O 2 pressure of 1 mTorr, x-ray diffraction (XRD) confirms FeWO 4 as the major crystalline phase for substrate temperatures from 500 to 800 °C. As temperature increases, x-ray fluorescence (XRF) shows the Fe/W ratio remains nearly constant at 0.90 ± 0.02, while x-ray photoelectron spectroscopy (XPS) shows Fe 3+ to Fe 2+ conversion increases from 20% to 35%. Morphological analysis reveals phase separation, likely of amorphous Fe 3+ oxide from crystalline FeWO 4 , increasing with Fe 3+ conversion. This is attributed to an O-rich laser ablated flux, where conversion is driven by the Arrhenius temperature dependence of Fe 3+ → Fe 2+ reduction on the film surface. At a constant substrate temperature of 750 °C, XRD confirms FeWO 4 formation for O 2 pressures from 0.5 to 10 mTorr. As pressure increases, XRF shows the Fe/W ratio decreases from 0.98 to 0.70, while XPS shows Fe 3+ conversion rises from 15% to 70%. Morphology shows phase separation decreasing with increasing Fe 3+ conversion. This is attributed to scattering, where higher O 2 pressure makes the laser ablated flux O-deficient relative to Fe and W, facilitating Fe 2+ formation. Films with Fe 3+ conversion above ∼30% and Fe/W ratios from 0.86 to 0.96 exhibit FeWO 4 optical and electronic properties suitable for photoanode applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction ( 1 ) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k 1 = 7.5 × 10 −12 cm 3 s −1 . High-level accuracy coupled cluster calculations in combination with E,J -resolved two-dimensional master equation analyses yield k 1 = (2.4 ± 0.5) × 10 −12 cm 3 s −1 for relevant atmospheric conditions ( T = 260–310 K and P = 0–1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction HCHO + H 2 O ⇌ HOCH 2 OH (2) with a HOCH 2 OH to HCHO ratio of only ca . 2%. Although HOCH 2 OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5–5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2–8.5 Tg/y for the global HCOOH production from cloud-derived HOCH 2 OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant K eq (2) determined in this work also allows us to estimate the Henry’s law constant of methanediol (8.1 × 10 5 M atm −1 at 280 K).

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

Effect of dielectric target properties on plasma surface ionization wave propagation

Surface ionization waves (SIWs) propagating along dielectric covered, grounded surfaces have been studied for various dielectric bulk and surface conditions; a dependence on the propagation velocity with respect to dielectric electrical thickness and near surface permittivity profiles are observed. SIWs generated by an atmospheric pressure plasma source are imaged interacting with planar dielectric surface. Surface wave velocity is obtained by tracking emission intensity as a function of time. Target dielectric thickness is varied from $d = 0.15-10$ mm and dielectric constant is varied from $\epsilon_r = 6.21 - 9.4$. The propagation of SIWs can be generally predicted by relating their velocity to the RC time constant of the circuit generated between the plasma and the dielectric surface, but it is found that this approximation breaks down for dielectric substrates of sufficient thickness and wave velocity becomes constant. The results show that wave velocity is stable and predictable for target thicknesses beyond a certain point determined by the permittivity of the target material. It is also shown that SIW propagation is strongly driven by the dielectric material near to the surface of the target in addition to the bulk material. The possible mechanisms driving these thickness dependent behaviors is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fiber-optic bolometer with low detection limit fabricated using thermal release tape and precise laser heating

Fiber-optic bolometers (FOBs) based on a fiber-tipped silicon Fabry–Perot (FP) interferometric temperature sensor and a gold disk absorber have been shown to be an attractive alternative to conventional resistive bolometers for plasma radiation measurement in fusion devices. Either a high-finesse FP or a low-finesse FP can be used, each with trade-offs between noise performance and fabrication complexity. In this paper, we present an FOB design that overcomes these limitations by combining a low-finesse long silicon FP cavity with a large gold disk absorber to achieve enhanced sensitivity and noise performance without increasing the fabrication complexity and the time constant. We also demonstrated a fabrication method for the sensor head facilitated by thermal release tape and precise laser heating. Our FOB demonstrates a temperature resolution of 0.08 mK, a cooling time constant of 230 ms, and a noise equivalent power density of 0.015 W m −2 . This represents an eightfold improvement over previous high-finesse FOBs and 26-fold improvement over previous low-finesse FOBs with similar demodulation bandwidths and similar cooling time constants.

Fabry–Perot interferometer↗
DOE OSTI2022

Vacuum energy density measured from cosmological data

Within the ΛCDM cosmological model, the absolute value of Einstein's cosmological constant Λ, sometimes expressed as the gravitating mass-energy density ρΛ of the physical vacuum, is a fundamental constant of nature, whose accurate measurement plays a central role in testing some proposed theories of quantum gravity. Several combinations of currently public cosmological data and an assumed flat ΛCDM cosmological model are used here to make a joint Bayesian inference on the combination of conventional parameters Ω Λ h 2 that corresponds to the absolute physical density ρΛ. In physical units, we obtain ρ Λ = (60.3 ± 1.3)× 10 -31 g/cm 3 , the most accurate constraint to date, with an absolute calibration of cosmological measurements based on CMB temperature. Significantly different values are obtained with calibrations that use a local distance scale, mainly connected to systematic differences in the value of the Hubble constant. Furthermore, it is suggested that future comprehensive cosmological parameter studies assuming the ΛCDM model include constraints on the vacuum density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A framework for testing soil carbon dynamics post land-use transition in a multisector dynamics model

Soil carbon plays a crucial role in the global carbon cycle. Changes in land use can determine whether carbon is stored or is emitted into the atmosphere as carbon dioxide, which has broad implications for the human and Earth systems. These feedbacks to the carbon cycle and their socio-economic drivers are modelled by many global multisector dynamics models to project future possibilities for the human-Earth system. One notable model of this class is the Global Change Analysis Model (GCAM), which uses a simplified process to model soil organic carbon (SOC) content after land-use transition across 384 land units. While the current GCAM soil carbon framework is based on scientific principles, it has not been tested against experimental data. This work examines rates of SOC change from GCAM input data. Specifically, first order rate constants derived from model inputs were compared to values from two syntheses to assess GCAM’s accuracy. Welch’s t-tests and linear models were used to determine if rate constants were consistent across all tested geographical areas and land-use transition types. While we found that there was general agreement on the direction and magnitude (i.e., rate) of SOC change, the rate constant derived from GCAM and empirical values differed strongly in a subset of specific instances. These results indicate that GCAM’s current SOC dynamics during land use transition successfully capture broad patterns of change in this critical carbon pool, but should be interpreted with caution at finer spatial scales. One potential cause of these discrepancies is our highly aggregated variable, soil timescale, which could be made more granular to improve accuracy. When using economically rooted multisector dynamics models, such as GCAM, it is critical to understand such model limitations for representing specific Earth system processes.

carbon↗
DOE OSTI2020

A viscoelastic model for seismic attenuation using fractal mechanical networks

SUMMARY Seismic attenuation (quantified by the quality factor Q) has a significant impact on the seismic waveforms, especially in the fluid-saturated rocks. This dissipative process can be phenomenologically represented by viscoelastic models. Previous seismological studies show that the Q value of Earth media exhibits a nearly frequency-independent behaviour (often referred to as constant-Q in literature) in the seismic frequency range. Such attenuation can be described by the mathematical Kjartansson constant-Q model, which lacks of a physical representation in the viscoelastic sense. Inspired by the fractal nature of the pore fluid distribution in patchy-saturated rocks, here we propose two fractal mechanical network (FMN) models, that is, a fractal tree model and a quasi-fractal ladder model, to phenomenologically represent the frequency-independent Q behaviour. As with the classic viscoelastic models, the FMN models are composed of mechanical elements (spring and dashpots) arranged in different hierarchical patterns. A particular parametrization of each model can produce the same complex modulus as in the Kjartansson model, which leads to the constant-Q. Applying the theory to several typical rock samples, we find that the seismic attenuation signature of these rocks can be accurately represented by either one of the FMN models. Besides, we demonstrate that the ladder model in particular exhibits the realistic multiscale fractal structure of the saturated rocks. Therefore, the FMN models as a proxy could provide a new way to estimate the microscopic rock structure property from macroscopic seismic attenuation observation.

Geochemistry & Geophysics↗
DOE OSTI2023

Moderate high temperature is beneficial or detrimental depending on carbon availability in the green alga Chlamydomonas reinhardtii

High temperatures impair plant growth and reduce agricultural yields, but the underlying mechanisms remain elusive. The unicellular green alga Chlamydomonas reinhardtii is an excellent model to study heat responses in photosynthetic cells due to its fast growth rate, many similarities in cellular processes to land plants, simple and sequenced genome, and ample genetic and genomics resources. Chlamydomonas grows in light by photosynthesis and with externally supplied acetate as an organic carbon source. Understanding how organic carbon sources affect heat responses is important for the algal industry but remains understudied. We cultivated wild-type Chlamydomonas under highly controlled conditions in photobioreactors at 25 °C (control), 35 °C (moderate high temperature), or 40 °C (acute high temperature) with or without constant acetate supply for 1 or 4 day. Treatment at 35 °C increased algal growth with constant acetate supply but reduced algal growth without sufficient acetate. The overlooked and dynamic effects of 35 °C could be explained by induced acetate uptake and metabolism. Heat treatment at 40 °C for more than 2 day was lethal to algal cultures with or without constant acetate supply. Furthermore, our findings provide insights to understand algal heat responses and help improve thermotolerance in photosynthetic cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

An outflow powers the optical rise of the nearby, fast-evolving tidal disruption event AT2019qiz

ABSTRACT At 66 Mpc, AT2019qiz is the closest optical tidal disruption event (TDE) to date, with a luminosity intermediate between the bulk of the population and the faint-and-fast event iPTF16fnl. Its proximity allowed a very early detection and triggering of multiwavelength and spectroscopic follow-up well before maximum light. The velocity dispersion of the host galaxy and fits to the TDE light curve indicate a black hole mass ≈106 M⊙, disrupting a star of ≈1 M⊙. By analysing our comprehensive UV, optical, and X-ray data, we show that the early optical emission is dominated by an outflow, with a luminosity evolution L ∝ t2, consistent with a photosphere expanding at constant velocity (≳2000 km s−1), and a line-forming region producing initially blueshifted H and He ii profiles with v = 3000–10 000 km s−1. The fastest optical ejecta approach the velocity inferred from radio detections (modelled in a forthcoming companion paper from K. D. Alexander et al.), thus the same outflow may be responsible for both the fast optical rise and the radio emission – the first time this connection has been observed in a TDE. The light-curve rise begins 29 ± 2 d before maximum light, peaking when the photosphere reaches the radius where optical photons can escape. The photosphere then undergoes a sudden transition, first cooling at constant radius then contracting at constant temperature. At the same time, the blueshifts disappear from the spectrum and Bowen fluorescence lines (N iii) become prominent, implying a source of far-UV photons, while the X-ray light curve peaks at ≈1041 erg s−1. Assuming that these X-rays are from prompt accretion, the size and mass of the outflow are consistent with the reprocessing layer needed to explain the large optical to X-ray ratio in this and other optical TDEs, possibly favouring accretion-powered over collision-powered outflow models.

Nicholl, M.↗
DOE OSTI2025

New scaling and nuclear structure aspects in heavy-ion fusion reactions

Three new behaviors have been found in comparisons of fusion cross sections for different collision systems. root (1) Replacing the energy E with a scaling one, E scal = (E-V g )/($\sqrt{2}$W g ), is successful for washing out the Coulomb interaction in the spectra of fusion cross sections, where V g and W g are barrier height and width of the single-Gaussian barrier distribution model. (2) In a representation of σE vs the scaling energy, E scal , all data sets display in parallel. Here, the ratio for sigma E from any two fusion systems over the whole range is a constant value. That behavior is also studied in another representation, in which the data sets display as parallel horizontal lines for any heavy-ion fusion system. (3) The constant ratio value is the ratio of parameter products, $R^2_gW_g$, of the two systems; where R g is the barrier radius obtained in the single-Gaussian barrier distribution model. Moreover, when comparing neighboring collision systems at the same E scal , the ratio of sigma is near a constant value within a few percent over the whole range. Thus a quantitative comparison for the fusion enhancement for neighboring systems is developed. The present finding could be beneficial for predicting unmeasured fusion cross sections.

Jiang, C. L. [Argonne National Laboratory (ANL), A↗
DOE OSTI2026

FLAG review 2024

We review lattice results related to pion, kaon, 𝐷-meson, 𝐵-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor 𝑓+⁡(0) arising in the semileptonic 𝐾 → 𝜋 transition at zero momentum transfer, as well as the decay-constant ratio 𝑓 𝐾 ⁡/𝑓 𝜋 and its consequences for the Cabibbo–Kobayashi–Maskawa (CKM) matrix elements 𝑉 𝑢⁢𝑠 and 𝑉 𝑢⁢𝑑 . We review the determination of the 𝐵 𝐾 parameter of neutral kaon mixing as well as the additional four 𝐵 parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for 𝑚 𝑐 and 𝑚 𝑏 as well as those for the decay constants, form factors, and mixing parameters of charmed and bottom mesons and baryons. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant 𝛼 𝑠 . We review the determinations of nucleon charges from the matrix elements of both isovector and flavor-diagonal axial, scalar and tensor local quark bilinears, and momentum fraction, helicity moment and the transversity moment from one-link quark bilinears. We also review determinations of scale-setting quantities. Finally, in this review we have added a new section on the general definition of the low-energy limit of the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Scrambling versus relaxation in Fermi and non-Fermi liquids

In this work, we compute the Lyapunov exponent characterizing quantum scrambling in a family of generalized Sachdev- Ye-Kitaev models, which can be tuned between different low temperature states from Fermi liquids, through non-Fermi liquids, to fast scramblers. The analytic calculation, controlled by a small coupling constant and large $\textit{N}$, allows us to clarify the relations between the quasiparticle relaxation rate $\textit{1/τ}$ and the Lyapunov exponent $λ_L$ characterizing scrambling. In the Fermi liquid states we find that the quasiparticle relaxation rate dictates the Lyapunov exponent. In non-Fermi liquids, where $1/τ \gg T$, we find that $λ_L$ is always $\textit{T}$ linear with a prefactor that is independent of the coupling constant in the limit of weak coupling. Instead it is determined by a scaling exponent that characterizes the relaxation rate. $λ_L$ approaches the general upper bound 2$\textit{πT}$ at the transition to a fast scrambling state. Finally, in a marginal Fermi liquid state the exponent is linear in temperature with a prefactor that vanishes as a nonanalytic function ~$g\text{ln}(1/g)$ of the coupling constant $\textit{g}$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Condensation of collective polar vortex modes

The dynamics of extended objects such as domain walls, domain bubbles, vortex structures, etc., can be described by their equations of motion associated with their effective mass and spring constant. Here we analytically derive the equations of motion for the polarization dynamics and elastodynamics for the structural responses of ferroelectric polar vortices, and theoretically extract their effective mass, spring constant, and mode frequencies. We demonstrate two subterahertz phonon modes and predicted their frequencies, both consistent with our recent experimental measurements and phase-field simulations. We show that elastic modulation of the energy function and spring constants leads to a condensation of a collective mode upon a second-order structural transition from symmetric to asymmetric vortices at a critical strain, analogous to the ferroelectric soft phonon mode at a ferroelectric transition. Furthermore, the present work offers a theoretical framework for predicting and manipulating the ultrafast collective dynamics of polar nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Quasi-deuteron model at low renormalization group resolution

The quasi-deuteron model introduced by Levinger is used to explain cross sections for knocking out high-momentum protons in photoabsorption on nuclei. This is within a framework we characterize as exhibiting high renormalization group (RG) resolution. Assuming a one-body reaction operator, the nuclear wave function must include two-body short-range correlations (SRCs) with deuteronlike quantum numbers. In Phys. Rev. C 104, 034311 (2021), we showed that SRC physics can be naturally accounted for at low RG resolution. We describe the quasi-deuteron model at low RG resolution and determine the Levinger constant, which is proportional to the ratio of nuclear photoabsorption to that for photodisintegration of a deuteron. We extract the Levinger constant based on the ratio of momentum distributions at high relative momentum. We compute momentum distributions evolved under similarity RG (SRG) transformations where the SRC physics is shifted into the operator as a universal two-body term. The short-range nature of this operator motivates using local-density approximations with uncorrelated wave functions in evaluating nuclear matrix elements, which greatly simplifies the analysis. The operator must be consistently matched to the RG scale and scheme of the interaction for a reliable extraction. We apply SRG transformations to different nucleon-nucleon (NN) interactions and use the deuteron wave functions and Weinberg eigenvalues to determine approximate matching scales. We predict the Levinger constant for several NN interactions and a wide range of nuclei comparing to experimental extractions. The predictions at low RG resolution are in good agreement with experiment when starting with a hard NN interaction and the initial operator. Similar agreement is found using soft NN interactions when the additional two-body operator induced by evolution from hard to soft is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

BICEP/ K e c k XIV: Improved constraints on axionlike polarization oscillations in the cosmic microwave background

We present an improved search for axionlike polarization oscillations in the cosmic microwave background (CMB) with observations from the Keck Array . An all-sky, temporally sinusoidal rotation of CMB polarization, equivalent to a time-variable cosmic birefringence, is an observable manifestation of a local axion field and potentially allows a CMB polarimeter to detect axionlike dark matter directly. We describe improvements to the method presented in previous work, and we demonstrate the updated method with an expanded dataset consisting of the 2012–2015 observing seasons. We set limits on the axion-photon coupling constant for mass m in the range 10 - 23 – 10 - 18 eV , which corresponds to oscillation periods on the order of hours to years. Our results are consistent with the background model. For periods between 1 and 30 d ( 1.6 × 10 - 21 ≤ m ≤ 4.8 × 10 - 20 eV ), the 95%-confidence upper limits on rotation amplitude are approximately constant with a median of 0.27°, which constrains the axion-photon coupling constant to g φ γ < ( 4.5 × 10 - 12 GeV - 1 ) m / ( 10 - 21 eV ) , if axionlike particles constitute all of the dark matter. More than half of the collected BICEP dataset has yet to be analyzed, and several current and future CMB polarimetry experiments can apply the methods presented here to achieve comparable or superior constraints. In the coming years, oscillation measurements can achieve the sensitivity to rule out unexplored regions of the axion parameter space.

79 ASTRONOMY AND ASTROPHYSICS↗