Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “importance sampling”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Supported molybdenum oxides for the aldol condensation reaction of acetaldehyde

The (retro-)aldol condensation reaction is an important chemical transformation in the upgrading of biomass-derived compounds into fuels and valuable specialty chemicals. In this study, we found that supported molybdenum oxide (MoO x ) catalysts were active and selective for the aldol condensation of acetaldehyde to crotonaldehyde under steady-state reactor conditions. Through a combination of transmission electron microscopy (TEM), ultraviolet–visible (UV–VIS) diffuse reflectance spectroscopy, Fourier transform infrared (FTIR) spectroscopy of adsorbed pyridine, and steady-state reactor testing, we determined that highly dispersed MoO x has a strong interaction with a γ-Al 2 O 3 support resulting in optimal catalyst performance at low weight loadings. In contrast, MoO x particles supported on SiO 2 have a weaker interaction with the support, resulting in a monotonic relationship between Mo loading and aldol condensation activity. The Lewis acid site density and strength are important parameters for predicting aldol condensation activity across all samples. The concentration of weak acid sites had a poor correlation with aldol condensation activity, most likely because these sites are too weak to activate acetaldehyde for the reaction. Medium and strong acid sites both had good correlations to aldol condensation activity. Results from X-ray absorption near edge structure (XANES) and acetaldehyde temperature programmed desorption (TPD) indicated that partially reduced MoO x was more active for aldol condensation, but pretreatment in reducing or oxidizing environments had no significant effect on steady-state catalytic activity. Finally, characterization of spent catalyst samples through temperature programmed oxidation (TPO) and thermogravimetric analysis (TGA) revealed that catalysts with high densities of strong acid sites tended to form more carbonaceous deposits on the surface over the course of the reaction.

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

DNA viruses from different stages of a wastewater treatment plant in southwest Ohio

Wastewater treatment plants (WWTPs) house diverse populations of microbial communities that are dynamic across treatment stages. Although viruses, especially non-pathogenic ones are relatively less studied compared to their bacterial counterparts, they are equally important players in WWTP microbiomes. In this study, we sampled influent, activated sludge, and effluent stages of a WWTP in southwest Ohio to sequence their viral population using multiple displacement amplification (MDA) and metagenomic approaches. We recovered 1003 viral genomes, the majority of which were ssDNA (single-stranded DNA) viruses that formed distinct clusters representing novel species and phylogenetic groups not closely related to known human gut or terrestrial viruses. Additionally, tracking viruses across treatment stages showed several persisted across multiple treatment stages. These results suggest that further studies are needed to understand the persistence of viral populations, particularly non-pathogenic ones, across different stages of WWTPs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Unequal {110} Facets: The Potential Role of Intraparticle Heterogeneity and Facet Termination in Photoelectrochemical Activity of Single BiVO 4 Particles

BiVO 4 photoanodes are promising for solar water splitting, with photogenerated electrons and holes preferentially reacting at top {010} and lateral {110} facets, respectively. However, the mechanisms driving this facet-dependent reactivity remain unclear. Here, in this study, we investigate facet-dependent photocurrent and material heterogeneity using correlative scanning photoelectrochemical microscopy (SPCM), electron beam induced current (EBIC) mapping, and mid-IR scattering scanning near-field optical microscopy (s-SNOM). SPCM measurements of 62 BiVO 4 particles confirmed higher photocurrents at lateral {110} facets compared to top {010} facets, but unexpectedly revealed variations in photocurrent among lateral facets within the same particle. Variations in lateral facet surface termination could explain the intraparticle-level reactivity heterogeneity, consistent with theoretical predictions. Nano-FTIR spectroscopy and Raman microspectroscopy indicated significant materials chemistry heterogeneity within individual particles and facets that could be attributed to variations in lattice vibration distortions that enhance the overlap between Bi 6s and O 2p orbitals. The increased orbital overlap is significant as it potentially increases hole mobility in the valence band and potentially explains the lateral facet-dependent charge separation efficiency observed in photocurrent maps. Facet-dependent electrical and EBIC measurements showed no space charge regions at interfacet junctions or metal-BiVO 4 contacts under vacuum, suggesting that photogenerated holes beneath top {010} facets are unlikely to transport to lateral {110} facets to drive water/sulfite oxidation. These findings indicate the potential influence of distinct bulk properties and surface termination chemistries across different particles and facets, highlighting the importance of carefully controlling defects and surface chemistry during sample growth to optimize photocatalytic performance.

BiVO4↗
DOE OSTI2024

Use of Ionic Liquids for Hemp Fiber Degumming

Degumming or retting to remove pectin, lignin, and hemicellulose is a critical step in the production of high-quality hemp fibers. Current technologies, including traditional water or dew retting, as well as chemical, enzymatic, and physical degumming processes, have inherent environmental risks and produce fibers of varying qualities. Herein we explore the use of an ionic liquid (IL) capable of dissolving pectin, lignin, and hemicellulose, but not cellulose, in combination with different oxidating agents to degum hemp. In this study, we observed that treating raw hemp fibers with the IL choline acetate ([Cho][OAc]), alone or in combination with aqueous H 2 O 2 , led to a reduction in the gummy materials and lignin content of the resulting samples of up to 63 and 78%, respectively. More importantly, the high crystallinity and low reduction in tensile strength of the resulting fibers make these novel and relatively simple IL-based processes promising alternatives toward the sustainable preparation of hemp-based materials suitable for use in the textile industry.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Influence of growth temperature of KTiOAsO{sub 4} single crystals on their physicochemical parameters and formation of domain structures

A potassium titanyl arsenate (KTiOAsO{sub 4}, KTA) crystal 50 × 80 × 60 mm in size has been grown by upgraded Czochralski method from flux (TGGS) with a decrease in temperature from 900 to 770 °C during pulling. It is shown that the spectroscopic properties of the parts of KTA crystals grown at 900 and 770 °C are close, whereas the electrical conductivity of the low-temperature (770 °C) KTA part turned out to be an order of magnitude lower than that of the high-temperature part. Visualisation of the domain structure by second-harmonic generation microscopy revealed a more efficient domain intergrowth (throughout the sample) in the low-temperature KTA, which is important for forming a regular domain structure (RDS) in a KTA-based nonlinear optical element. It is established that the quantum efficiency of parametric generation of light in the RDS formed in low-temperature KTA is several times higher than in the case of high-temperature KTA. The results obtained are important for optimising RDS parameters. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A tale of two planets: Disparate evolutionary models for Mars inferred from radiogenic isotope compositions of Martian meteorites

The radiogenic isotopic compositions of basaltic Martian meteorites (shergottites) and clinopyroxene/olivine cumulate meteorites (nakhlite/chassignites) are used to define the global evolution of Mars. However, the two main groups of meteorites demonstrate that their sources underwent divergent styles of magmatic evolution. The shergottites portray a planet that differentiated ~4.52 billion years ago via solidification of a magma ocean, producing incompatible element-depleted and -enriched reservoirs that remained isolated until melt production. In contrast, the reservoir from which the nakhlite/chassignites derive may have formed earlier, produced melts that fractionated Sm/Nd and Hf/W differently, was compositionally less variable, and experienced a significantly more complex history following primordial differentiation than the shergottite sources. The disparate histories recorded by these two groups of meteorites elucidate important questions that could be addressed by acquiring additional samples. Obtaining samples that shared the isotopic systematics of the shergottites would provide confidence that extrapolating the primordial differentiation history of Mars from shergottite radiogenic isotope systematics is reasonable. Returned samples from Mars will also constrain the physical locations of the meteorite source regions, providing insights into the general structure of the Martian mantle. In addition, they will help constrain the phases present in the martian mantle during melting and the conditions under which they are stable. Finally, identifying an evolved lithology that satisfies the geochemical and isotopic constraints placed on the incompatible element-enriched endmember observed in the shergottites would define the nature of magmatic evolution on Mars and whether it is more akin to processes on the Earth or the Moon.

Mars↗
DOE OSTI2023

Micro- and Nanoscale Surface Analysis of Late Iron Age Glass from Broborg, a Vitrified Swedish Hillfort

Abstract Archaeological glasses with prolonged exposure to biogeochemical processes in the environment can be used to understand glass alteration, which is important for the safe disposal of vitrified nuclear waste. Samples of mafic and felsic glasses with different chemistries, formed from melting amphibolitic and granitoid rocks, were obtained from Broborg, a Swedish Iron Age hillfort. Glasses were excavated from the top of the hillfort wall and from the wall interior. A detailed microscopic, spectroscopic, and diffraction study of surficial textures and chemistries were conducted on these glasses. Felsic glass chemistry was uniform, with a smooth surface showing limited chemical alteration (<150 nm), irrespective of the position in the wall. Mafic glass was heterogeneous, with pyroxene, spinel, feldspar, and quartz crystals in the glassy matrix. Mafic glass surfaces in contact with topsoil were rougher than those within the wall and had carbon-rich material consistent with microbial colonization. Limited evidence for chemical or physical alteration of mafic glass was found; the thin melt film that coated all exposed surfaces remained intact, despite exposure to hydraulically unsaturated conditions, topsoil, and associated microbiome for over 1,500 years. This supports the assumption that aluminosilicate nuclear waste glasses will have a high chemical durability in near-surface disposal facilities.

Matthews, Bethany E. (ORCID:0000000171932583)↗
DOE OSTI2020

Accurate and efficient description of interacting carriers in quantum nanostructures by selected configuration interaction and perturbation theory

In this paper, we present a method to accurately and efficiently calculate many-body states of interacting carriers in quantum nanostructures based on a combination of iterative selection of configurations and perturbation theory. This method enables investigations of large excitonic complexes and multielectron systems with near full configuration interaction accuracy, even though only a small subspace of the full many-body Hilbert space is sampled, thus saving orders of magnitudes in computational resources. Important advantages of this method are that the convergence is controlled by a single parameter, the threshold, and that ground and excited states can be treated on an equal footing. On the example of InAsP nanowire quantum dots described using a million-atom tight-binding approach, we demonstrate the extreme efficiency of the method by numerical studies of large excitonic complexes filling up to the fourth electronic shell. We find that the method generally converges fast as a function of the threshold, profiting from a significant acceleration due to the perturbative corrections. The role of the choice of single-particle basis states is discussed. It is found that the algorithm converges faster in the Hartree-Fock basis only for highly charged systems, where Coulomb repulsion dominates. Finally, based on the observation that second-order perturbative energy corrections only depend on off-diagonal elements of the many-body Hamiltonian, we present a way to accurately calculate many-body states that requires only a relatively small number of Coulomb matrix elements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fire, insect and disease‐caused tree mortalities increased in forests of greater structural diversity during drought

Abstract Structural diversity is an emerging dimension of biodiversity that accounts for size variations in organs among individuals in a community. Previous studies show significant effects of structural diversity on forest growth, but its effects on forest mortality are not known, particularly at a large scale. To address this knowledge gap, we quantified structural diversity using stem structural diversity (SSD) based on both tree diameter and height. We obtained U.S. Forest Service Forest Inventory and Analysis (FIA) data from over 2400 plots across southcentral U.S. forests that have suffered a recent drought. Using data from multiple sampling times, we calculated SSD and compared the relative importance of SSD, species diversity, functional diversity and other stand attributes in determining tree mortalities caused by fire, insects and diseases. We also used FIRETEC, a physics‐based fire model, to test the effect of SSD on canopy consumption by fire. Our results showed that (1) SSD was positively associated with tree mortalities caused by all three disturbances; (2) species richness was negatively associated with insect‐ and disease‐caused mortalities; (3) functional diversity was negatively associated with fire‐ and disease‐caused mortalities and (4) more phylogenetically related species had more similar mortality rates by insect and disease but not fire. Moreover, the FIRETEC model showed increasing canopy consumption by fire in stands with greater SSD. Together, the different tree mortalities during drought associated with SSD more consistently than the other biodiversity metrics were evaluated. Synthesis . Our results suggest that SSD could be considered in modelling forest dynamics and planning management to sustain forest health under disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Power Production, Inter- and Intra-Array Wake Losses from the U.S. East Coast Offshore Wind Energy Lease Areas

There is an urgent need to develop accurate predictions of power production, wake losses and array–array interactions from multi-GW offshore wind farms in order to enable developments that maximize power benefits, minimize levelized cost of energy and reduce investment uncertainty. New, climatologically representative simulations with the Weather Research and Forecasting (WRF) model are presented and analyzed to address these research needs with a specific focus on offshore wind energy lease areas along the U.S. east coast. These, uniquely detailed, simulations are designed to quantify important sources of wake-loss projection uncertainty. They sample across different wind turbine deployment scenarios and thus span the range of plausible installed capacity densities (ICDs) and also include two wind farm parameterizations (WFPs; Fitch and explicit wake parameterization (EWP)) and consider the precise WRF model release used. System-wide mean capacity factors for ICDs of 3.5 to 6.0 MWkm−2 range from 39 to 45% based on output from Fitch and 50 to 55% from EWP. Wake losses are 27–37% (Fitch) and 11–19% (EWP). The discrepancy in CF and wake losses from the two WFPs derives from two linked effects. First, EWP generates a weaker ‘deep array effect’ within the largest wind farm cluster (area of 3675 km2), though both parameterizations indicate substantial within-array wake losses. If 15 MW wind turbines are deployed at an ICD of 6 MWkm−2 the most heavily waked wind turbines generate an average of only 32–35% of the power of those that experience the freestream (undisturbed) flow. Nevertheless, there is no evidence for saturation of the resource. The wind power density (electrical power generation per unit of surface area) increases with ICD and lies between 2 and 3 Wm−2. Second, EWP also systematically generates smaller whole wind farm wakes. Sampling across all offshore wind energy lease areas and the range of ICD considered, the whole wind farm wake extent for a velocity deficit of 5% is 1.18 to 1.38 times larger in simulations with Fitch. Over three-quarters of the variability in normalized wake extents is attributable to variations in freestream wind speeds, turbulent kinetic energy and boundary layer depth. These dependencies on meteorological parameters allow for the development of computationally efficient emulators of wake extents from Fitch and EWP.

17 WIND ENERGY↗
DOE OSTI2023

Reconstructing Ly$α$ Fields from Low-resolution Hydrodynamical Simulations with Deep Learning

Hydrodynamical cosmological simulations are a powerful tool for accurately predicting the properties of the intergalactic medium (IGM) and for producing mock skies that can be compared against observational data. However, the need to resolve density fluctuation in the IGM puts a stringent requirement on the resolution of such simulations, which in turn limits the volumes that can be modeled, even on the most powerful supercomputers. In this work, we present a novel modeling method that combines physics-driven simulations with data-driven generative neural networks to produce outputs that are qualitatively and statistically close to the outputs of hydrodynamical simulations employing eight times higher resolution. We show that the Ly$α$ flux field, as well as the underlying hydrodynamic fields, have greatly improved statistical fidelity over a low-resolution simulation. Importantly, the design of our neural network allows for sampling multiple realizations from a given input, enabling us to quantify the model uncertainty. Using test data, we demonstrate that this model uncertainty correlates well with the true error of the Ly$α$ flux prediction. Ultimately, our approach allows for training on small simulation volumes and applying it to much larger ones, opening the door to producing accurate Ly$α$ mock skies in volumes of Hubble size, as will be probed with DESI and future spectroscopic sky surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A High-Granularity Approach to Modeling Energy Consumption and Savings Potential in the U.S. Residential Building Stock: Preprint

Building simulations are increasingly used in various applications related to energy efficient buildings. For individual buildings, applications include: design of new buildings, prediction of retrofit savings, ratings, performance path code compliance and qualification for incentives. Beyond individual building applications, larger scale applications (across the stock of buildings at various scales: national, regional and state) include: codes and standards development, utility program design, regional/state planning, and technology assessments. For these sorts of applications, a set of representative buildings are typically simulated to predict performance of the entire population of buildings. Focusing on the U.S. single-family residential building stock, this paper will describe how multiple data sources for building characteristics are combined into a highly-granular database that preserves the important interdependencies of the characteristics. We will present the sampling technique used to generate a representative set of thousands (up to hundreds of thousands) of building models. We will also present results of detailed calibrations against building stock consumption data.

building stock↗
DOE OSTI2025

Field expedient stool collection methods for gut microbiome analysis in deployed military environments

ABSTRACT Field expedient devices and protocols for the collection, storage, and shipment of stool samples in deployed settings are needed for the advancement of microbiome research in military health. Relevant assessments include the evaluation of microbiome signatures associated with susceptibility to travelers’ diarrhea and recovery of gut function following infection. However, inherent biases in microbial measurements due to preservatives and sampling methods are unclear and should be assessed for an accurate evaluation of the microbiome. We performed shotgun metagenomic sequencing and compared the microbiome composition in paired fecal samples collected using Flinters Technology Associates (FTA) cards and OMNIgene (OG) Gut tubes, prior to and during international travel, from 49 adult participants, 39 of whom remained asymptomatic and 10 experienced travelers’ diarrhea. Higher concentrations of nucleic acid and sequencing libraries were observed in OG samples. A majority of genera (82.9%) were detected with both methods, and detections of genera limited to one collection method were not highly prevalent across samples and were present in extremely low relative abundances (<0.01%). Differences in beta diversity were largely explained by inter-individuality of microbiome composition, followed by the effect of collection method and timepoint-disease states. Differential abundance analysis indicated that Corynebacterium and Blautia were consistently higher in abundance across all groups with FTA and OG collection, respectively. The observed differences in microbiome composition between methods suggest the need for consistent and standardized protocols within a study. Overall, the data presented here could help guide the future design of fecal microbiome study protocols in field and military deployment settings. IMPORTANCE The assessment of field-deployable methods for fecal sample collection and storage is required to reliably capture samples collected in remote and austere locations. This study describes a comparative metagenomics analysis between samples collected by two different commercially available methods in a military-deployed setting. The results presented here are foundational for the future design of fecal microbiome study protocols in an operational context.

field study↗
DOE OSTI2023

Effects of finite sampling on fatigue damage estimation of wind turbine components: A statistical study

The variability of the wind turbine loads complicates fatigue assessment in the design phase, as performing simulations covering the entire lifetime is computationally expensive. The current work provides important information for assessing the uncertainty in fatigue damage estimation due to finite data. Here we study the sample size effect on mean, variance, and skewness of damage in each wind bin, identify the important wind bins, and study the uncertainty propagation from each wind bin to the lifetime damage using 3600 aeroelastic simulations and bootstrapping. To achieve less than 1% error in the damage estimation across all load channels in the current case study, at least 100 turbulence seeds are needed. Damage in different wind bins follows a lognormal distribution when using the conventional approach of six seeds. The provided insights and information allow the designer to achieve a specific level of accuracy for a given computational cost using strategic bin sampling.

17 WIND ENERGY↗
DOE OSTI2022

Illustrating galaxy–halo connection in the DESI era with illustrisTNG

ABSTRACT We employ the hydrodynamical simulation illustrisTNG to inform the galaxy–halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ∼ 0.8. Specifically, we model the galaxy colours of illustrisTNG and apply sliding DESI colour–magnitude cuts, matching the DESI target densities. We study the halo occupation distribution (HOD) model of the selected samples by matching them to their corresponding dark matter haloes in the illustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as ‘galaxy assembly bias’. The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modelling galaxy–halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10–15 ${{\ \rm per\ cent}}$ excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Complete resolution across the neodymium/samarium isotopic envelope with a liquid sampling‐atmospheric pressure glow discharge — Orbitrap mass spectrometer

Rationale Nd and Sm isotope ratios play an important role in geological dating and as nuclear forensic signatures; however, the overlap of the respective 144, 148, 150 Nd/Sm isobars requires prior separations to be performed before analysis on typical MS platforms. The work presented here overcomes these isobaric interferences using ultrahigh‐mass resolution to alleviate interference without prior chemical separations. Methods A liquid sampling‐atmospheric pressure glow discharge ion source was coupled to a standard, QExactive Focus Orbitrap mass spectrometer, providing a mass resolution of ~80 k. A Spectroswiss FTMS booster X2 data acquisition package was used to collect extended transients, providing much higher mass resolution; ~230 k and ~600 k are employed here for Nd and Sm isotopes. Results While the standard Orbitrap resolution is far greater than typical “atomic” MS platforms, it was insufficient to alleviate all isobars. The use of a resolution of ~230 k resulted in baseline separation across the entire isotopic envelope for both Nd and Sm. Isotope ratios obtained from Nd:Sm mixtures using high‐resolution were equivalent to those found for individual‐element solutions, while isotope ratios obtained at a resolution of ~80 k (standard for the OEM data system) showed large deviations. Conclusions Use of ultrahigh‐resolution is an attractive alternative to extensive chemical separations to alleviate severe isobaric interferences. Sufficient mass resolution greatly reduces/eliminates the need for sample manipulations (separations) before analysis while reducing costs and total analysis times.

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

Grover-QAOA for 3-SAT: quadratic speedup, fair-sampling, and parameter clustering

Abstract The SAT problem is a prototypical NP-complete problem of fundamental importance in computational complexity theory with many applications in science and engineering; as such, it has long served as an essential benchmark for classical and quantum algorithms. This study shows numerical evidence for a quadratic speedup of the Grover Quantum Approximate Optimization Algorithm (G-QAOA) over random sampling for finding all solutions to 3-SAT (All-SAT) and Max-SAT problems. G-QAOA is less resource-intensive and more adaptable for these problems than Grover’s algorithm, and it surpasses conventional QAOA in its ability to sample all solutions. We show these benefits by classical simulations of many-round G-QAOA on thousands of random 3-SAT instances. We also observe G-QAOA advantages on the IonQ Aria quantum computer for small instances, finding that current hardware suffices to determine and sample all solutions. Interestingly, a single-angle-pair constraint that uses the same pair of angles at each G-QAOA round greatly reduces the classical computational overhead of optimizing the G-QAOA angles while preserving its quadratic speedup. We also find parameter clustering of the angles. The single-angle-pair protocol and parameter clustering significantly reduce obstacles to classical optimization of the G-QAOA angles.

Zhang, Zewen (ORCID:000000032258613X)↗
DOE OSTI2021

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces

The scattering-type scanning near-field optical microscope (s-SNOM) has emerged as a powerful tool for resolving nanoscale inhomogeneities in laterally heterogeneous samples. However, most analytical models used to predict the scattering near-field signals are assuming homogenous landscapes (bulk materials), resulting in inconsistencies when applied to samples with more complex configurations. In this work, we combine the point-dipole model (PDM) to the finite-element method (FEM) to account for the lateral and vertical heterogeneities while keeping the computation time manageable. Full images, spectra, or hyperspectral line profiles can be simulated by calculating the self-consistent dipole radiation demodulated at higher harmonics of the tip oscillation, mimicking real experimental procedures. Using this formalism, we clarify several important yet puzzling experimental observations in near-field images on samples with rich typography and complex material compositions, heterostructures of two-dimensional material flakes, and plasmonic antennas. The developed method serves as a basis for future investigations of nano-systems with nontrivial topography.

36 MATERIALS SCIENCE↗