Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 901 records · Page 50

A Radiative Transfer Simulator for PACE: Theory and Applications

A radiative transfer simulator was developed to compute the synthetic data of all three instruments onboard NASA’s Plankton Aerosol, Cloud, ocean Ecosystem (PACE) observatory, at the top of the atmosphere (TOA). The instrument suite includes the ocean color instrument (OCI), the HyperAngular Rainbow Polarimeter 2 (HARP2), and the Spectro-Polarimeter for Planetary Exploration 1 (SPEXone). The PACE simulator is wrapped around a monochromatic radiative transfer model based on the successive order of scattering (RTSOS), which accounts for atmosphere and ocean coupling, polarization, and gas absorption. Inelastic scattering, including Raman scattering from pure ocean water, fluorescence due to chlorophyll, and colored dissolved organic matter (CDOM), is also simulated. This PACE simulator can be used to explore the sensitivity of the hyperspectral and polarized reflectance of the Earth system with tunable atmosphere and ocean parameters, which include aerosol and cloud number concentration, refractive indices, and size distribution, ocean particle microphysical parameters, and solar and sensor-viewing geometry. The PACE simulator is used to study two important case studies. One is the impact of the significant uncertainty in pure ocean water absorption coefficient to the radiance field in the ultraviolet (UV) spectral region, which can be as much as 6%. The other is the influence of different amounts of brown carbon aerosols and CDOM on the polarized radiance field at TOA. The percentage variation of the radiance field due to CDOM is mostly for wavelengths smaller than 600 nm, while brown aerosols affect the whole spectrum from 350 to 890 nm, primarily due to covaried soot aerosols. Both case studies are important for aerosol and ocean color remote sensing and have not been previously reported in the literature.

PACE↗

Measurements of Nuclear Effects and the $\bar \nu_\mu +H \to \mu^+ + n$ Cross Section in MINERνA with Neutron Tagging

MINERνA, or Main INjector ExpeRiment for ν-A, at Fermilab, is an experimentdedicated to the study of neutrino-nucleus interactions in the GeV regime. Itsgoal is to illustrate the interplay between hadronic and nuclear physics and measure intranuclear dynamics crucial for the present and future neutrino oscillationmeasurements. We first measure a set of variables sensitive to how Monte Carlo(MC) simulations of neutrino-nucleus interactions implement binding energy andthen move on to measure the antineutrino CCQE cross section on the hydrogentargets in MINERνA’s CH detector. We have developed a method to preferentially select events on the hydrogen by comparing outgoing neutrons’ directions totheoretical neutron directions assuming two-body interactions. We measured thecross section, extracted the axial form factor, and performed a Z-expansion fit. Weobserve larger values in the axial form factor at high Q2than current best fits. Finally, we show a preliminary selection of events with both protons and neutrons toinvestigate nuclear processes responsible for producing these final states.

Cai, Tejin↗

Parameterizing Perennial Bioenergy Crops in Version 5 of the Community Land Model Based on Site‐Level Observations in the Central Midwestern United States

Abstract With projected expansion of biofuel production at a global scale, there is a pressing need to develop adequate representation of bioenergy crops in land surface models to help effectively quantify the biogeophysical and biogeochemical effects of its associated land use changes. This study implements two new perennial bioenergy crops, Miscanthus and switchgrass, into the Community Land Model Version 5 based on site‐level observations from the midwestern United States by modifying parameters associated with photosynthesis, phenology, allocation, decomposition, and carbon cost of nitrogen uptake and integrating concomitantly land management practices. Sensitivity analyses indicate that carbon and energy fluxes of the perennial crops are most sensitive to photosynthesis and phenology parameters. Validation of simulated fluxes against site‐level measurements demonstrates that the model is capable of capturing the overall patterns of energy and carbon fluxes, as well as physiological transitions from leaf emergence to senescence. Compared to annual crops, perennial crops feature longer growing season, greater leaf areas, and higher productivity, leading to increased transpiration, lower annual runoff, and larger carbon uptake. The model simulations suggest that with higher CO 2 assimilation rates and lower demands for nutrients and water, high‐yielding perennial crops are promising alternatives of bioenergy feedstocks compared to traditional annual crops not only for mitigating climate change but also for environmental conservation purposes by reducing fertilizer application and therefore alleviating surface‐ and ground‐water contaminations. Although the local‐scale simulations shed light on potential benefits of using perennial grasses as bioenergy feedstocks, quantifying consequences of their plantations at larger scales warrants additional investigation.

Cheng, Yanyan↗

Discerning Deception: An Empirically-Driven Agent-Based Model of Expert Evaluation of Scientific Content

Both human subject experiments and computational, modeling and simulations have been used to study detection of deception. This work aims to combine these two methods by integrating empirically-derived information (from human subject experiments) into agent-based models to generate novel insights into the complex problems of detection of disinformation content. Computational experiments are used to simulate across multiple scenarios for evaluation and decision-making regarding the validity of potentially deceptive scientific documents. Factors influencing the human agent behaviors in the model were identified through a human subject experiment that was conducted to evaluate and characterize decision making related to disinformation discernment. Correlation and regression analyses were used to translate insights from the human subjects experiment to inform the parameterization of agent features and scenario development. Three scenarios were evaluated with the agent-based models to help evaluate the replicability of the simulations (validation analysis) and assess the influence of human agent and document features (sensitivity analyses). A replication of the human participant experiment demonstrated that the agent-based simulations compare favorably to empirical findings. The agent-based modeling was then used to conduct sensitivity analysis on the accuracy of deception detection as a function of document proportions and human agent features. Results indicate that precision values are adversely impacted when the proportion of deceptive documents is lower in the overall sample, whereas recall values are more sensitive to changes in human agent features. These findings indicate important nuances in accuracy evaluations that should be further considered (including consideration of potential alternate metrics) in future agent-based models of disinformation. Additional areas for future exploration include extension of simulations to consider other ways to align the agent-based model design with psychological theory and inclusion of agent-agent interactions, especially as it pertains to sharing of scientific information within an organizational context.

99 GENERAL AND MISCELLANEOUS↗

Adaptive Planning: Understanding Organizational Workload to Capability/ Capacity through Modeling and Simulation

In August 2003, the Secretary of Defense (SECDEF) established the Adaptive Planning (AP) initiative [1] with an objective of reducing the time necessary to develop and revise Combatant Commander (COCOM) contingency plans and increase SECDEF plan visibility. In addition to reducing the traditional plan development timeline from twenty-four months to less than twelve months (with a goal of six months)[2], AP increased plan visibility to Department of Defense (DoD) leadership through In-Progress Reviews (IPRs). The IPR process, as well as the increased number of campaign and contingency plans COCOMs had to develop, increased the workload while the number of planners remained fixed. Several efforts from collaborative planning tools to streamlined processes were initiated to compensate for the increased workload enabling COCOMS to better meet shorter planning timelines. This paper examines the Joint Strategic Capabilities Plan (JSCP) directed contingency planning and staffing requirements assigned to a combatant commander staff through the lens of modeling and simulation. The dynamics of developing a COCOM plan are captured with an ExtendSim [3] simulation. The resulting analysis provides a quantifiable means by which to measure a combatant commander staffs workload associated with development and staffing JSCP [4] directed contingency plans with COCOM capability/capacity. Modeling and simulation bring significant opportunities in measuring the sensitivity of key variables in the assessment of workload to capability/capacity analysis. Gaining an understanding of the relationship between plan complexity, number of plans, planning processes, and number of planners with time required for plan development provides valuable information to DoD leadership. Through modeling and simulation AP leadership can gain greater insight in making key decisions on knowing where to best allocate scarce resources in an effort to meet DoD planning objectives.

Hase, Chris↗

(abstract) Using TOPEX/Poseidon Sea Level Observations to Test the Sensitivity of an Ocean Model to Wind Forcing

It has been demonstrated that current-generation global ocean general circulation models (OGCM) are able to simulate large-scale sea level variations fairly well. In this study, a GFDL/MOM-based OGCM was used to investigate its sensitivity to different wind forcing. Simulations of global sea level using wind forcing from the ERS-1 Scatterometer and the NMC operational analysis were compared to the observations made by the TOPEX/Poseidon (T/P) radar altimeter for a two-year period. The result of the study has demonstrated the sensitivity of the OGCM to the quality of wind forcing, as well as the synergistic use of two spaceborne sensors in advancing the study of wind-driven ocean dynamics.

forcing radar altimeter global sea level↗

A Verification of Flux Sensitivity Estimates Using the MCNP Tally Perturbation Tool

Nuclear data is commonly used in applications such as nuclear nonproliferation, safeguards, and criticality safety. More specifically, nuclear data is used in predictive simulation codes like the Monte-Carlo N-Particle (MCNP ® ) transport code, Serpent, and similar radiation transport codes. The improvement of nuclear data enables more precise and accurate simulations, which result in higher fidelity designs and reduced operational/procedural costs. Therefore, the improvement of nuclear data is of paramount importance across the nuclear community. Nuclear data is improved and validated through integral benchmark experiments. The design of benchmark experiments is an extensive process; therefore, these experiments are often optimized on multiple characteristics, including sensitivity to the nuclear data, during the design process. Sensitivity is a measure of how much a quantity changes due to changes in independent variables such as experimental configuration. An experimental design that has a larger sensitivity to the nuclear data of interest will have a larger impact on the accuracy and precision of the validated data. Past integral benchmark experiments have primarily used the effective multiplication factor ($k_{eff}$) as the predominant measured quantity; however, experiments designed with other quantities in mind would be able to optimize on validating different areas of the nuclear data. A primary goal of the EUCLID project is to design, constrain, and reduce compensating errors in experiments focused on quantities other than $k_{eff}$ to better validate nuclear data across the board. Currently, there is a capability in MCNP to easily calculate the sensitivity of $k_{eff}$ to specific nuclear data of numerous reactions types and isotopes (KSEN card); however, the sensitivity of other quantities must be estimated in more strenuous manners. For example, the perturbation feature (PERT card) of MCNP can be used to estimate first-order sensitivities of some response in fixed source simulations. A recent announcement revealed that the first- and second-order perturbation features in previous releases of MCNP contained a bug. It was identified that particles were being scored into the wrong energy bin. The bug is in the most recent public release (MCNP6.2); however, a patch has been added to the most up to date version (MCNP6.2.2) that has not been released publicly. A direct comparison of the PERT card results for an F4 (neutron flux averaged over a cell) tally before and after the patch are shown in figure 1. All simulations used in the sensitivity estimates in this report were performed with MCNP6.2.2. This work verifies the patched MCNP perturbation tool by comparing first order sensitivities made using the PERT card to estimates made using manual perturbation of the compact ENDF (ACE) files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Improving Stratocumulus Cloud Amounts in a 200‐m Resolution Multi‐Scale Modeling Framework Through Tuning of Its Interior Physics

Abstract High‐Resolution Multi‐scale Modeling Frameworks (HR)—global climate models that embed separate, convection‐resolving models with high enough resolution to resolve boundary layer eddies—have exciting potential for investigating low cloud feedback dynamics due to reduced parameterization and ability for multidecadal throughput on modern computing hardware. However low clouds in past HR have suffered a stubborn problem of over‐entrainment due to an uncontrolled source of mixing across the marine subtropical inversion manifesting as stratocumulus dim biases in present‐day climate, limiting their scientific utility. We report new results showing that this over‐entrainment can be partly offset by using hyperviscosity and cloud droplet sedimentation. Hyperviscosity damps small‐scale momentum fluctuations associated with the formulation of the momentum solver of the embedded large eddy simulation. By considering the sedimentation process adjacent to default one‐moment microphysics in HR, condensed phase particles can be removed from the entrainment zone, which further reduces entrainment efficiency. The result is an HR that can produce more low clouds with a higher liquid water path and a reduced stratocumulus dim bias. Associated improvements in the explicitly simulated sub‐cloud eddy spectrum are observed. We report these sensitivities in multi‐week tests and then explore their operational potential alongside microphysical retuning in decadal simulations at operational 1.5° exterior resolution. The result is a new HR having desired improvements in the baseline present‐day low cloud climatology, and a reduced global mean bias and root mean squared error of absorbed shortwave radiation. We suggest it should be promising for examining low cloud feedbacks with minimal approximation.

54 ENVIRONMENTAL SCIENCES↗

Perturbation Modeling and Navigation Tuning for a Crewed Station in Near Rectilinear Halo Orbit

NASA's Gateway program will build a crew-tended station in an Earth-Moon Near Rectilinear Halo Orbit (NRHO) to support deep space missions to the lunar surface and heliocentric space. The station in the NRHO will be visited by crewed vehicles including lunar landers that are high mass relative to Gateway. Docking events with large vehicles induce significant thruster plume and docking contact velocity perturbations. Crewed vehicles induce significant venting perturbations that can stress the navigation filter's ability to remain converged about an estimated state. The torques that arise from venting and lunar gravity gradient require reaction control system (RCS) thruster intervention, which themselves induce a velocity perturbation. The best method to model RCS thruster perturbations is currently in development. This paper simulates different RCS perturbation models, schedules, and an increasing magnitude to investigate the navigation sensitivity to different RCS perturbation models. From the RCS perturbation models simulated, a worst-case scenario is chosen for further scrutiny. The navigation filter is further tuned to handle the most disruptive perturbation modeling and the resulting performance is analyzed.

Navigation↗

Dose Rate Analysis of the WCS Consolidated Interim Storage Facility

This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗

Injection of electron beams into two laser wakefields and generation of electron rings

Throughout this work, mutual injection of electron beams into two laser plasma wakefields was observed experimentally when driving laser pulses interfered in plasma at a small crossing angle and were slightly relatively delayed, approximately by one pulse duration. Particle-in-cell simulations revealed that the mutual injection was sensitive to the spatial overlap of the laser pulses, which therefore could be used to control the mutual injection. The dual synchronized, femtosecond electron beams are potentially useful for pump-probe experiments in ultrafast science. In addition, out-of-axis ring-shaped electron beams were detected in both experiments and simulations.

47 OTHER INSTRUMENTATION↗

Behavior and mechanisms of Doppler wind lidar error in complex terrain: stable flow case study at Perdigão

A numerical experiment is carried out investigating the magnitude of biases in ground-based lidar measurements in complex flow conditions. Biases assessed include those arising from flow curvature and from the interaction of turbulence with the wind field reconstruction (WFR) algorithms used by a WindCube lidars and anemometers. RANS-CFD and WRF-LES simulations were performed for the Perdig˜ao Field Experiment site for a range of atmospheric conditions. Virtual anemometer and lidar data were generated for four locations: two near exposed ridge tops and two in low-speed regions in the valley. The LES data at these four locations show that the scalar inflation terms (the relation between scalar and vector averaged wind speed) for virtual lidar and virtual cups agree very well with predictions using perturbation theory. While the lidar errors vary greatly with location and height, the contribution from the flow curvature tends to be larger than the differences arising from scalar inflation. For one lidar/mast pair near the ridge top, comparisons between simulations and measurements are carried out for a resonant mountain wave event on June 14th, 2017, and for the whole duration of the Perdigão campaign for winds perpendicular to the ridges. The lidar error during the mountain wave, a period of strong stability and low inversion height, is significantly larger than the campaign average. The sensitivity of the lidar error to atmospheric stability is confirmed by the RANS simulations, which suggests strong sensitivity of flow curvature error to stability conditions and to the shape of the wind speed profile near the top of the boundary layer.

17 WIND ENERGY↗

Sensitivity of a Cloud-Resolving Model to the Bulk and Explicit Bin Microphysical Schemes

A cloud-resolving model is used to study sensitivities of two different microphysical schemes, one is the bulk type, and the other is an explicit bin scheme, in simulating a mid-latitude squall line case. Simulations using different microphysical schemes are compared with each other and also with the observations. Both the bulk and bin models reproduce the general features during the developing and mature stage of the system. The leading convective zone, the trailing stratiform region, the horizontal wind flow patterns, pressure perturbation associated with the storm dynamics, and the cool pool in front of the system all agree well with the observations. Both the observations and the bulk scheme simulation serve as validations for the newly incorporated bin scheme. However, it is also shown that the bulk and bin simulations have distinct differences, most notably in the stratiform region. Weak convective cells exist in the stratiform region in the bulk simulation, but not in the bin simulation. These weak convective cells in the stratiform region are remnants of the previous stronger convections at the leading edge of the system. The bin simulation, on the other hand, has a horizontally homogeneous stratiform rain structure, which agrees better with the observations. Further analysis of the downdraft core strength, the potential temperature perturbation, and the evaporative cooling rate shows that the differences between the bulk and bin models are due mainly to the stronger low-level evaporative cooling simulated in the bulk model. Further quantitative analysis and sensitivity tests for this case using both the bulk and bin models will be presented in the talk.

Tao, Wei-Kuo↗

Sensitivity of a Cloud-Resolving Model to the Bulk and Explicit Bin Microphysical Schemes: Validations with a PRE-STORM Case - Part 1

A cloud-resolving model is used to study sensitivities of two different microphysical schemes, one is the bulk type, and the other is an explicit bin scheme, in simulating a mid-latitude squall line case (PRE-STORM, June 10-11, 1985). Simulations using different microphysical schemes are compared with each other and also with the observations. Both the bulk and bin models reproduce the general features during the developing and mature stage of the system. The leading convective zone, the trailing stratiform region, the horizontal wind flow patterns, pressure perturbation associated with the storm dynamics, and the cool pool in front of the system all agree well with the observations. Both the observations and the bulk scheme simulation serve as validations for the newly incorporated bin scheme. However, it is also shown that, the bulk and bin simulations have distinct differences, most notably in the stratiform region. Weak convective cells exist in the stratiform region in the bulk simulation, but not in the bin simulation. These weak convective cells in the stratiform region are remnants of the previous stronger convections at the leading edge of the system. The bin simulation, on the other hand, has a horizontally homogeneous stratiform cloud structure, which agrees better with the observations. Preliminary examinations of the downdraft core strength, the potential temperature perturbation, and the evaporative cooling rate show that the differences between the bulk and bin models are due mainly to the stronger low-level evaporative cooling in convective zone simulated in the bulk model. Further quantitative analysis and sensitivity tests for this case using both the bulk and bin models will be presented in a companion paper.

Li, Xiao-Wen↗

Combining Satellite Microwave Radiometer and Radar Observations to Estimate Atmospheric Latent Heating Profiles

In this study, satellite passive microwave sensor observations from the TRMM Microwave Imager (TMI) are utilized to make estimates of latent + eddy sensible heating rates (Q1-QR) in regions of precipitation. The TMI heating algorithm (TRAIN) is calibrated, or "trained" using relatively accurate estimates of heating based upon spaceborne Precipitation Radar (PR) observations collocated with the TMI observations over a one-month period. The heating estimation technique is based upon a previously described Bayesian methodology, but with improvements in supporting cloud-resolving model simulations, an adjustment of precipitation echo tops to compensate for model biases, and a separate scaling of convective and stratiform heating components that leads to an approximate balance between estimated vertically-integrated condensation and surface precipitation. Estimates of Q1-QR from TMI compare favorably with the PR training estimates and show only modest sensitivity to the cloud-resolving model simulations of heating used to construct the training data. Moreover, the net condensation in the corresponding annual mean satellite latent heating profile is within a few percent of the annual mean surface precipitation rate over the tropical and subtropical oceans where the algorithm is applied. Comparisons of Q1 produced by combining TMI Q1-QR with independently derived estimates of QR show reasonable agreement with rawinsonde-based analyses of Q1 from two field campaigns, although the satellite estimates exhibit heating profile structure with sharper and more intense heating peaks than the rawinsonde estimates. 2

Grecu, Mircea↗

Detectability of subsolar mass neutron stars through a template bank search

Here, we study the detectability of gravitational-wave signals from subsolar-mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as 78.4% for an equal mass binary system of chirp mass 0.17M ⊙ , in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for subsolar-mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of subsolar-mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for subsolar-mass binary neutron star signals.

79 ASTRONOMY AND ASTROPHYSICS↗

Cloud-radiative effects on implied oceanic energy transport as simulated by atmospheric general circulation models

This paper summarizes the ocean surface net energy flux simulated by fifteen atmospheric general circulation models constrained by realistically-varying sea surface temperatures and sea ice as part of the Atmospheric Model Intercomparison Project. In general, the simulated energy fluxes are within the very large observational uncertainties. However, the annual mean oceanic meridional heat transport that would be required to balance the simulated surface fluxes is shown to be critically sensitive to the radiative effects of clouds, to the extent that even the sign of the Southern Hemisphere ocean heat transport can be affected by the errors in simulated cloud-radiation interactions. It is suggested that improved treatment of cloud radiative effects should help in the development of coupled atmosphere-ocean general circulation models.

Gleckler, P. J.↗

Imaging Photonic Resonances within an All‐Dielectric Metasurface via Photoelectron Emission Microscopy

Dielectric metasurfaces, through volume‐type photonic resonances, enable precise control of light‐matter interactions for applications including imaging, holography, and sensing. The application space of dielectric metasurfaces has extended from infrared to visible wavelengths by incorporating high refractive index materials, such as titanium dioxide (TiO 2 ). Understanding the fundamental and fabrication limits for these applications requires metrology with nanoscale resolution, sensitivity to electromagnetic fields within the meta‐atom volume, and far‐field excitation. In this work, photoelectron emission microscopy (PEEM) is used to image field distributions of photonic resonances in a TiO 2 metasurface excited with far‐field, visible‐wavelength illumination. The local volumetric field variations within the meta‐atoms are analyzed as a function of illumination angle and polarization by comparing photoelectron images to finite‐difference time‐domain simulations. This study determines the inelastic mean free path of very low‐energy (<1 eV) photoelectrons to be 35 ± 10 nm, which is comparable to the meta‐atom height thereby highlighting PEEM sensitivity to resonances within the volume. Additionally, the simulations reveal high sensitivity of PEEM images to an in‐plane component of the illumination k ‐vector. These results demonstrate that photoelectron imaging with subwavelength resolution offers unique advantages for examining light‐matter interactions in volume‐type (as opposed to surface) photonic modes within dielectric nanophotonic structures.

electron inelastic mean free path↗