Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency estimation”

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

Improved accuracy and robustness of electron density profiles from JET’s X-mode frequency-modulated continuous-wave reflectometers

JET’s frequency-modulated continuous wave (FMCW) reflectometers have been operating well with the current design since 2005, and density profiles have been automatically calculated intershot since then. However, the calculated profiles had long suffered from several shortcomings: poor agreement with other diagnostics, sometimes inappropriately moving radially by several centimeters, elevated levels of radial jitter, and persistent wriggles (strong unphysical oscillations). Here, in this research, several techniques are applied to the reflectometry data analysis, and the shortcomings are significantly improved. Starting with improving the equilibrium reconstruction that estimates the background magnetic field, adding a ripple correction in the reconstructed magnetic field profile, and adding new inner-wall reflection positions estimated through ray-tracing, these changes not only improve the agreement of reconstructed profiles to other diagnostics but also solve density profile wriggles that were present during band transitions. Other smaller but also persistent wriggles were also suppressed by applying a localized correction to the measured beat frequency where persistent oscillations are present. Finally, the burst analysis method, as introduced by Varela et al., has been implemented to extract the beat frequency from stacked spectrograms. Due to the strong suppression of spurious reflections, the radial jitter that sometimes would span several centimeters has been strongly reduced. The stacking of spectrograms has also been shown to be very useful for stacking recurring events, like small gas puff modulations, and extracting transport coefficients that would otherwise be below the noise level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Design and characterization of a single-channel microwave interferometer for the Helicon Physics Prototype eXperiment

A heterodyne interferometer for measuring the plasma electron density in Helicon Physics Prototype eXperiment (HPPX) has been developed and characterized in the laboratory. The heterodyne-type interferometer uses a probing beam of 85 GHz and a reference beam of 85.7 GHz by using two Gunn oscillators for facilitating frequency modulation. The cutoff density for this frequency is 8.96 x 10 19 m –3 that is much higher than the nominal electron density of the HPPX plasmas. The maximum deflection angle due to the refraction effect is estimated to be 0.63° assuming a parabolic electron density profile, and the impact on the measurement is negligible for the plasmas with a size of 100 mm. Two identical HDPE lenses act as the vacuum windows, and the lenses working together with the transmitter and receiver horns, form a beam with the beam waist of about 20 mm located at the plasma center. The characterization results of the lens and the critical microwave components show good performance. The plasma electron density measurement results are in line with the expectation. Furthermore, the Signal to Noise Ratio (SNR) of the measurement is defined as SNR = 20lg($\bar{n}_{e}$/S), where $\bar{n}_{e}$ and S are respectively the average electron density and the standard deviation of the signal, and the SNR is 80 dB. Finally, this demonstrates the system's capability of providing a clear and detailed density resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High-fidelity qutrit entangling gates for superconducting circuits

Ternary quantum information processing in superconducting devices poses a promising alternative to its more popular binary counterpart through larger, more connected computational spaces and proposed advantages in quantum simulation and error correction. Although generally operated as qubits, transmons have readily addressable higher levels, making them natural candidates for operation as quantum three-level systems (qutrits). Recent works in transmon devices have realized high fidelity single qutrit operation. Nonetheless, effectively engineering a high-fidelity two-qutrit entanglement remains a central challenge for realizing qutrit processing in a transmon device. In this work, we apply the differential AC Stark shift to implement a flexible, microwave-activated, and dynamic cross-Kerr entanglement between two fixed-frequency transmon qutrits, expanding on work performed for the ZZ interaction with transmon qubits. We then use this interaction to engineer efficient, high-fidelity qutrit CZ † and CZ gates, with estimated process fidelities of 97.3(1)% and 95.2(3)% respectively, a significant step forward for operating qutrits on a multi-transmon device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Delensing the CMB with the cosmic infrared background: the impact of foregrounds

ABSTRACT The most promising avenue for detecting primordial gravitational waves from cosmic inflation is through measurements of degree-scale cosmic microwave background (CMB) B-mode polarization. This approach must face the challenge posed by gravitational lensing of the CMB, which obscures the signal of interest. Fortunately, the lensing effects can be partially removed by combining high-resolution E-mode measurements with an estimate of the projected matter distribution. For near-future experiments, the best estimate of the latter will arise from co-adding internal reconstructions (derived from the CMB itself) with external tracers such as the cosmic infrared background (CIB). In this work, we characterize how foregrounds impact the delensing procedure when CIB intensity, I, is used as the matter tracer. We find that higher point functions of the CIB and Galactic dust such as 〈BEI〉c and 〈EIEI〉c can, in principle, bias the power spectrum of delensed B-modes. To quantify these, we first estimate the dust residuals in currently available CIB maps and upcoming, foreground-cleaned Simons Observatory CMB data. Then, using non-Gaussian simulations of Galactic dust – extrapolated to the relevant frequencies, assuming the spectral index of polarized dust emission to be fixed at the value determined by Planck – we show that the bias to any primordial signal is small compared to statistical errors for ground-based experiments, but might be significant for space-based experiments probing very large angular scales. However, mitigation techniques based on multifrequency cleaning appear to be very effective. We also show, by means of an analytical model, that the bias arising from the higher point functions of the CIB itself ought to be negligible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simultaneous Multiwavelength Flare Observations of EV Lacertae

Here we present the first results of our ongoing project conducting simultaneous multiwavelength observations of flares on nearby active M dwarfs. We acquired data of the nearby dM3.5e star EV Lac using five different observatories: NASA’s Transiting Exoplanet Survey Satellite (TESS), NASA’s Neil Gehrels Swift Observatory (Swift), NASA’s Neutron Interior Composition Explorer (NICER), the University of Hawaii 2.2-meter telescope (UH88), and the Las Cumbres Observatory Global Telescope (LCOGT) Network. During the ~25 days of TESS observations, we acquired three simultaneous UV/X-ray observations using Swift that total ~18 ks, 21 simultaneous epochs totaling ~98 ks of X-ray data using NICER, one observation (~3 hr) with UH88, and one observation (~3 hr) with LCOGT. We identified 56 flares in the TESS light curve with estimated energies in the range log E T (erg) = (30.5–33.2), nine flares in the Swift UVM2 light curve with estimated energies in the range log E UV (erg) = (29.3–31.1), 14 flares in the NICER light curve with estimated minimum energies in the range log E N (erg) = (30.5–32.3), and 1 flare in the LCOGT light curve with log E L (erg) = 31.6. We find that the flare frequency distributions (FFDs) of TESS and NICER flares have comparable slopes, β T = -0.67 ± 0.09 and β N = - 0.65 ± 0.19, and the FFD of UVOT flares has a shallower slope (β U = -0.38 ± 0.13). Furthermore, we do not find conclusive evidence for either the first ionization potential (FIP) or the inverse FIP effect during coronal flares on EV Lac.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transverse emittance growth due to rf noise in crab cavities: Theory, measurements, cure, and high luminosity LHC estimates

The High-Luminosity LHC (HL-LHC) upgrade with planned operation from 2029 onward has a goal of achieving a tenfold increase in the integrated number of recorded collisions thanks to a doubling of the intensity per bunch ( 2.2 × 10 11 protons) and a reduction of β * (the β value in the two high luminosity detectors, namely ATLAS and CMS) to 15 cm. Such an increase in recorded collisions would significantly expedite new discoveries and exploration. Crab cavities are an important component of the HL-LHC upgrade and will contribute strongly to achieving an increase in the number of recorded collisions. However, noise injected through the crab cavity radio frequency (rf) system could cause significant transverse emittance growth and limit luminosity lifetime. We presented a theoretical formalism relating transverse emittance growth to rf noise in an earlier work. In this follow-up paper, we summarize measurements in the super-proton synchrotron (SPS) at CERN that validate the theory, we present estimates of the emittance growth rates using state-of-the-art rf and low-level rf (LLRF) technologies, and we set the rf noise specifications to achieve acceptable performance. A novel dedicated feedback system acting through the crab cavities to mitigate emittance growth will be required. In this work, we develop a theoretical formalism to evaluate the performance of such a feedback system in any collider, identify limiting components, present simulation results to validate these studies, and derive key design parameters for an HL-LHC implementation of such a feedback system. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Wind plant wake losses: Disconnect between turbine actuation and control of plant wakes with engineering wake models

Wake losses from neighboring plants may become a major factor in wind plant design and control as additional plants are constructed in areas with high wind resource availability. Because plant wakes span a large range of physical scales, from turbine rotor diameter to tens of kilometers, it is unclear whether conventional wake models or turbine control strategies are effective at the plant scale. Wake steering and axial induction control are evaluated in the current work as means of reducing the impact of neighboring wind plants on power and levelized cost of electricity. FLOw Redirection and Induction in Steady State (FLORIS) simulations were performed with the Gauss–Curl Hybrid and TurbOPark wake models as well as two operation and maintenance models to investigate control setpoint sensitivity to wake representation and economic factors. Both wake models estimate losses across a range of atmospheric conditions, although the wake loss magnitude is dependent on the wake model. Annual energy production and levelized cost of electricity are driven by wind direction frequency, with frequently aligned plants experiencing the greatest losses. However, both wake steering and axial induction are unable to mitigate the impact of upstream plants. Wake steering is constrained by plant geometry, since wake displacement is much less than the plant wake width, while axial induction requires curtailing the majority of turbines in upstream plants. Individual turbine strategies are limited by their effective scale and model representation. New wake models that include plant-scale physics are needed to facilitate the design of effective plant wake control strategies.

Scott, Ryan (ORCID:0000000328107574)↗
DOE OSTI2020

The Potential to Improve the Value of U.S. Geothermal Electricity Generation Through Flexible Operations

Geothermal power plants have typically been operated as baseload plants. The recent expansion of wind and solar power generation creates a potential opportunity to increase the value of geothermal generation through flexible operations. In recent years, California’s wholesale electricity markets have exhibited frequent, short-term periods of negative pricing, indicating that additional flexibility would be valued in the market. Here, we examine local nodal hourly price records at all geothermal plants located in the western United States. We describe how the frequency and temporal characteristics of negative pricing episodes have changed over recent years. Based on these price series, we calculate the value of multiple strategies of flexible operations. Additionally, we use the estimates of future prices, developed through a capacity-expansion model and a dispatch model, to explore how the value of such flexible operations might change along with further penetration of variable renewable power sources. Based on the historical pricing records, we find that simple curtailment of operations during negative pricing episodes could increase the average energy value by 1–2 dollars per MWh and that allowing for increased production during limited high-priced hours, in addition to curtailment during negative priced hours, could potentially double the increase in value (up to 4 dollars per MWh). The forward-looking simulations indicate reduced values of flexibility from geothermal power despite higher penetrations of variable renewable energy. This result highlights the possibility that increasing flexibility options throughout the system may counteract the influence of increased variable renewable energy deployment.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Bi-Level Linear Programming Model for Automatic Load Shedding: A Distributed Wide-Area Measurement System-based Solution

Load shedding is currently implemented as a two-step based approach. In the first step, manual load shedding is taken place, were system operators, using estimates, inform distribution utilities of predicted stressful conditions. Information provided include the potential use of energy reserves, as well as load shedding amount. In a second step, automatic load shedding is done. The latter is realized using protection relays. While considering frequency variation, pre-defined values of load to be shed and correspondent number of stages for such to be realized are transformed into relay settings. Under-frequency protection relays use only local measurements towards decision making, thus operate in a decentralized architecture. Decision making is done in milliseconds plus breaker time. While this approach has provided much system reliability, considering the new smart grid paradigm, where system dynamics are much faster due to increasing renewable resources penetration, in some operating conditions it will generate sub-optimal solutions, such as islanding. Phasor measurement units provide a source of information which can be useful for this problem. Centralized architecture-based solutions for automatic load shedding, as present in the state-of-the-art, require though total processing times which are not acceptable for real-life implementation. In this work, considering the above, a bi-level linear programming model is presented. The model is implemented considering a distributed architecture while leveraging phasor measurement units data. The upper-level model estimates the current system state. Results of this model are embedded in a lower-level model, which decision variables are the location and load value to be shed. Easy-to-implement model, built-on the classic weighted least squares solution, highlight potential aspects towards real-life applications.

Bretas, Arturo Suman↗
DOE OSTI2020

Cost Estimation of a Permanent Magnet Synchronous Machine for Use in Adjustable Speed-Pumped Storage Hydropower

Electric machine design has become an easier process with the availability of computers in recent decades. However, the interest in electric machine design among students has declined. University-level design courses must compete with courses perceived as more desirable by the potential job market. Hence, the power system engineers with this expertise are very few and are nearing retirement. Contrarily, renewable energy has been picking up in the past few decades. The indeterminacy and variation of renewable energy sources require a balancing mechanism to maintain the constant voltage and frequency of the grid. Adjustable-Speed Pumped-Storage Hydropower plants are considered the contender for future energy storage, providing large capacity and fast-acting response. Large-scale hydropower applications depend on the natural characteristics of each project site, making less expensive, off-the-shelf generators unfeasible. This work gives a shortcut to estimate the cost of constructing a permanent magnet generator for low-volume manufacturing.

13 HYDRO ENERGY↗
DOE OSTI2026

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗
DOE OSTI2022

New Lightning-Derived Vertical Total Electron Content Data Provide Unique Global Ionospheric Measurements

A newly released, novel ionospheric data set of global gridded vertical total electron content (VTEC) is introduced in this paper. This VTEC data set, provided by Los Alamos National Laboratory (LANL), is derived from very high frequency (VHF; defined as 30–300 MHz) broadband radio-frequency (RF) measurements of lightning made by U.S. Department of Defense sensing systems on board Global Positioning System (GPS) satellites. This paper presents the new data set (LANL VTEC), discusses the errors inherent in VHF TEC estimation due to ionospheric dispersion, and compares the LANL VTEC to two community standard VTEC gridded products: Jet Propulsion Laboratory's Global Ionospheric Model (JPL GIM) and the CEDAR community's Open Madrigal VTEC gridded measurements of L-band GNSS (global navigation satellite systems) TEC. We find that the LANL VTEC data have an offset of 3 TECU from CEDAR Madrigal GNSS VTEC and a full-width-half-maximum (FWHM) of 6 TECU. In comparison, the offset between LANL VTEC and the JPL GIM model is –3 TECU, but with a FWHM of 5 TECU. We also compare to Jason-3 VTEC measurements over the ocean, finding an offset of less than 0.5 TECU and an FWHM of <5 TECU. Because this technique uses a completely different methodology to determine TEC, the sources of errors are distinct from the typical ground-based GNSS L-band (GHz) TEC measurements. Also, because it is derived from RF lightning signals, this data set provides measurements in regions that are not well covered by ground-based GPS measurements, such as over oceans and over central Africa.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Laser absorption measurements of electron density in nanosecond-scale atmospheric pressure pulsed plasmas

We report on time-resolved measurements of electron number density by continuous-wave laser absorption in a low-energy nanosecond-scale laser-produced spark in atmospheric pressure air. Laser absorption is a result of free-free and bound-free electron excitation, with the absorption coefficient modeled and evaluated using estimates of the time-variation in electron temperature and probe laser absorption path length. Plasma electron number densities are determined to be as high as n e = 7 x 10 19 cm -3 , and decay to 1/e of their peak values over a period of about 50 ns following plasma formation using a 20 mJ, 10 ns pulse width frequency-doubled Nd:YAG laser. The measured plasma densities at later times are shown to be in reasonable agreement with Stark broadening measurements of the 3s[ 5 S o ]-3p[ 5 P] electronic transition in atomic oxygen at 777 nm. This study provides support for the use of such continuous wave laser absorption for time resolved electron density measurements in low energy spark discharges in air, provided that an estimate of the electron temperature and laser path length can be made by accompanying diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Estimating total attenuation using Rayleigh targets at cloud top: applications in multilayer and mixed-phase clouds observed by ground-based multifrequency radars

At millimeter wavelengths, attenuation by hydrometeors, such as liquid droplets or large snowflakes, is generally not negligible. When using multifrequency ground-based radar measurements, it is common practice to use the Rayleigh targets at cloud top as a reference in order to derive attenuation-corrected reflectivities and meaningful dual-frequency ratios (DFRs). By capitalizing on this idea, this study describes a new quality-controlled approach that aims at identifying regions of cloud where particle growth is negligible. The core of the method is the identification of a “Rayleigh plateau”, i.e., a large enough region near cloud top where the vertical gradient of DFR remains small. By analyzing co-located Ka–W band radar and microwave radiometer (MWR) observations taken at two European sites under various meteorological conditions, it is shown how the resulting estimates of differential path-integrated attenuation (ΔPIA) can be used to characterize hydrometeor properties. When the ΔPIA is predominantly produced by cloud liquid droplets, this technique alone can provide accurate estimates of the liquid water path. When combined with MWR observations, this methodology paves the way towards profiling the cloud liquid water, quality-flagging the MWR retrieval for rain and drizzle contamination, and/or estimating the snow differential attenuation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impedance‐matching model‐in‐the‐loop simulation

Abstract Impedance matching is a novel approach to executing model‐in‐the‐loop (MIL) simulations. The concept, key features, and sample applications of this approach have been described in a series of publications. This paper advances theory and implementation for impedance‐matching MIL by: (1) developing a feedback linearization technique for implementing MIL for nonlinear virtual systems, (2) systematically characterizing MIL performance‐stability tradeoffs using three design‐oriented criteria—performance: defined as closeness of the controlled‐actuator impedance to the virtual‐boundary impedance;control effort: low controller response at high frequencies to reduce sensitivity to measurement noise and uncertainties in the actuator model; andstability: pursued using passivity of the controlled actuator system, and (3) assessing, through experimentation, the robustness of two ways to estimate feedback reaction force from the test article. These fundamental advancements to MIL are described using a one‐dimensional configuration consisting of a fluid‐filled cylindrical vessel as the physical test article and a uniaxial hydraulic shake table controlled to imitate different seismic isolation systems at the base of the vessel.

Engineering↗
DOE OSTI2020

Planck intermediate results

In this work, we describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution (SED), which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck’s three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP’s data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. This parameter can be used to extract sub-samples of high-frequency sources with statistically well-understood properties. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1% of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. We describe the characteristics of this specific high-quality subset of PCCS2 and its validation against other data sets, specifically for: the sub-sample of PCCS2 located in low-cirrus areas; the Planck Catalogue of Galactic Cold Clumps; the Herschel GAMA15-field catalogue; and the temperature- and spectral-index-reconstructed dust maps obtained with Planck’s Generalized Needlet Internal Linear Combination method. The results of the BeeP extension of PCCS2, which are made publicly available via the Planck Legacy Archive, will enable the study of the thermal properties of well-defined samples of compact Galactic and extragalactic dusty sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Rippled metamaterials with scale-dependent tailorable elasticity

Thermally induced ripples are intrinsic features of nanometer-thick films, atomically thin materials, and cell membranes, significantly affecting their elastic properties. Despite decades of theoretical studies on the mechanics of suspended thermalized sheets, controversy still exists over the impact of these ripples, with conflicting predictions about whether elasticity is scale-dependent or scale-independent. Experimental progress has been hindered so far by the inability to have a platform capable of fully isolating and characterizing the effects of ripples. This knowledge gap limits the fundamental understanding of thin materials and their practical applications. Here, we show that thermal-like static ripples shape thin films into a class of metamaterials with scale-dependent, customizable elasticity. Utilizing a scalable semiconductor manufacturing process, we engineered nanometer-thick films with precisely controlled frozen random ripples, resembling snapshots of thermally fluctuating membranes. Resonant frequency measurements of rippled cantilevers reveal that random ripples effectively renormalize and enhance the average bending rigidity and sample-to-sample variations in a scale-dependent manner, consistent with recent theoretical estimations. The predictive power of the theoretical model, combined with the scalability of the fabrication process, was further exploited to create kirigami architectures with tailored bending rigidity and mechanical metamaterials with delayed buckling instability.

Applied Physical Sciences↗