Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level Set Method”

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 397 records · Page 22

Data-Driven Computation of Probabilistic Marching Cubes for Efficient Visualization of Level-Set Uncertainty

Uncertainty visualization is an important emerging research area. Being able to visualize data uncertainty can help scientists improve trust in analysis and decision-making. However, visualizing uncertainty can add computational overhead, which can hinder the efficiency of analysis. In this paper, we propose novel data-driven techniques to reduce the computational requirements of the probabilistic marching cubes (PMC) algorithm. PMC is an uncertainty visualization technique that studies how uncertainty in data affects level-set positions. However, the algorithm relies on expensive Monte Carlo (MC) sampling for the multivariate Gaussian uncertainty model because no closed-form solution exists for the integration of multivariate Gaussian. In this work, we propose the eigenvalue decomposition and adaptive probability model techniques that reduce the amount of MC sampling in the original PMC algorithm and hence speed up the computations. Our proposed methods produce results that show negligible differences compared with the original PMC algorithm demonstrated through metrics, including root mean squared error, maximum error, and difference images. We demonstrate the performance and accuracy evaluations of our data-driven methods through experiments on synthetic and real datasets.

Athawale, Tushar↗

Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: IV. Vibrationally-resolved ionization of the ground and excited electronic states

Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact ionization of molecular hydrogen and its isotopologues (H 2 , D 2 , T 2 , HD, HT, and DT) in both the ground and excited electronic states. We apply the adiabatic-nuclei molecular convergent close-coupling (MCCC) method to calculate cross sections from threshold to 1000 eV for ionization of the ground and excited vibrational levels of the X 1 Σ$^{+}_{g}$, B 1 Σ$^{+}_{u}$, C 1 π u , EF 1 Σ$^{+}_{g}$, a 3 Σ$^{+}_{g}$, and c 3 π u electronic states, representing all states with united-atoms-limit principle quantum number n=1–2. The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files. The data can also be downloaded from the MCCC database at mccc-db.org.

74 ATOMIC AND MOLECULAR PHYSICS↗

Transmission Interface Limits for High-Spatial Resolution Capacity Expansion Modeling

Large-scale capacity expansion models typically rely on estimates of the power transfer limits between modeled zones. Accurate estimation of these interface transfer limits (ITLs) requires modeling the underlying transmission network. Here we expand on a maximum flow optimization method that uses linearized power flow to estimate transfer limits. We apply this method to a data set of the U.S. transmission network to estimate ITLs between U.S. counties. By calculating ITLs using different subsets of the network, we evaluate how the size of the network used in the estimation affects the results. The results show diminishing returns to ITL accuracy after six hops, suggesting that a network subset can reasonably be used to approximate ITLs. The county-level estimates produced in this study will support more spatially resolved capacity expansion modeling and will help inform policy making at local and national levels.

capacity planning↗

A Bilevel Voltage Regulation Operation for Distribution Systems With Self-Operated Microgrids

The emerging of microgrids in distribution systems has significantly enhanced the resilience of power grids. However, the operators of a distribution system and microgrids therein can be different and have accessibility to different devices. To model the operation of such a grid, this work proposes a bilevel formulation and probes into the voltage regulation operation, considering the interaction between different systems. The proposed bilevel formulation considers the cooperation of active energy resources (AER), transformer tap-changers, and capacitor banks that are controlled by different operators. To facilitate the solution time of the target bilevel optimization, the lower-level problems with different objectives are modeled using deep neural networks (DNNs) which are then converted into a set of constraints. Hence, the bilevel problem can be reformed to a single-level problem. Lastly, the proposed solution procedures are validated using a customized joint system constructed by the IEEE 123-bus system and a real distribution system in Iowa. According to the numerical validation results, the solution time of the proposed nonlinear activation function based DNN model is 69 times faster than other methods in solving voltage regulation with a bilevel structure.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Data‐driven identification of environmental variables influencing phenotypic plasticity to facilitate breeding for future climates

Summary Phenotypic plasticity describes a genotype's ability to produce different phenotypes in response to different environments. Breeding crops that exhibit appropriate levels of plasticity for future climates will be crucial to meeting global demand, but knowledge of the critical environmental factors is limited to a handful of well‐studied major crops. Using 727 maize ( Zea mays L.) hybrids phenotyped for grain yield in 45 environments, we investigated the ability of a genetic algorithm and two other methods to identify environmental determinants of grain yield from a large set of candidate environmental variables constructed using minimal assumptions. The genetic algorithm identified pre‐ and postanthesis maximum temperature, mid‐season solar radiation, and whole season net evapotranspiration as the four most important variables from a candidate set of 9150. Importantly, these four variables are supported by previous literature. After calculating reaction norms for each environmental variable, candidate genes were identified and gene annotations investigated to demonstrate how this method can generate insights into phenotypic plasticity. The genetic algorithm successfully identified known environmental determinants of hybrid maize grain yield. This demonstrates that the methodology could be applied to other less well‐studied phenotypes and crops to improve understanding of phenotypic plasticity and facilitate breeding crops for future climates.

Kusmec, Aaron↗

Study on Electrostatic Separation of Quinoline Insolubles from Coal Tar Pitch

The feasibility of electrical separation in the removal of quinoline insoluble (QI) particles from coal tar pitch (CTP) was experimentally investigated. QI particle involvement prohibits the effective fabrication of high-quality value-added products, such as carbon fibers, from CTP. A substantial and sustainable CTP market exists around the world; therefore, a strong incentive exists to develop a viable technical approach to effectively and economically remove QI particles from CTP. The electrical separation method shows promise to achieve this technical goal (Cao et al., 2012). Even with the given setup (wire- cylinder adapted), critical issues remain to be addressed for this method to be applicable to the CTP: (1) identifying the wash oil used in the original QI separation, (2) understanding the mechanism of QI separation, (3) characterizing the deposit, and (4) identifying the QI. This study integrated a set of experimental and analysis techniques into the electrical separation tests to address the aforementioned issues. Key findings are as follows: 1. Electric current responses: • The electric current level of the CTP–wash oil solution reported by Cao et al. (2012) can be attained by using a mixture of 25% quinoline and 75% toluene for the CTPs examined in this study under the same electrical load condition. • The electric current tends to decrease during the test period because of the decreasing number of charged particles. • The reversed field corresponded to the configuration of electrostatic precipitation for positive corona discharge. The high electric field can result in dielectric breakdown and lead to an abrupt current surge. 2 Deposit response and solvent candidates: • Deposit of particles in CTP mixture can be effectively implemented in a wire-cylinder configuration as proposed and examined in this study. • Deposit depends on the solvents. Among the solvents tested, two- and three-part solvents that included quinoline and ethanol (i.e., 25% quinoline and 75% toluene; 50% quinoline and 50% toluene; and 33% wash oil, 33% BTX, and 34% ethanol) produced the highest deposit weight. • The deposit process examined in this study is derived from the charged particles. An appreciable relation exists between deposit weight and electric charge. 3. QI removal efficiency (RE): • The RE of the electrostatic separation can reach as high as 76.2% for Carbores, and 61.6% for Koppers. • The RE can be further enhanced if the field level increases from 0.16 kV/mm that was used in the current study to 0.23 kV/mm, according to the relation established between the RE and the electric field. 4. Testing of the modeling system: • The motion of particles is originally driven by the charge-based electric force in the cases tested. The solid particle separation mechanism is similar to that of QI deposition in CTP. • The mechanical movement of solid particles can be strongly affected by the gravitational and viscous forces in the electric field. 5. EDS analysis and chemical compositions: • The main chemical composition of deposit QI matches that of as-prepared CTP QI. The deposit QI is derived from the same group of the CTP QI. • In addition to carbon, the QI contains oxygen, sodium, aluminum, silicon, iron, and sulfur. The work for the near future is also discussed.

01 COAL, LIGNITE, AND PEAT↗

Transmission Interface Limits for High-Spatial Resolution Capacity Expansion Modeling

Large-scale capacity expansion models typically rely on estimates of the power transfer limits between modeled zones. Accurate estimation of these interface transfer limits (ITLs) requires modeling the underlying transmission network. Here we expand on a maximum flow optimization method that uses linearized power flow to estimate transfer limits. We apply this method to a data set of the U.S. transmission network to estimate ITLs between U.S. counties. By calculating ITLs using different subsets of the network, we evaluate how the size of the network used in the estimation affects the results. The results show diminishing returns to ITL accuracy after six hops, suggesting that a network subset can reasonably be used to approximate ITLs. The county-level estimates produced in this study will support more spatially resolved capacity expansion modeling and will help inform policy making at local and national levels.

capacity planning↗

Multi-Pressure Rail System Design with Variable Pressure Control Strategy

This paper presents the design and the control method of a Multiple Pressure Rail system (MPR) with a variable pressure control logic that allows minimization of throttling loss, therefore, increasing the overall fluid power actuation efficiency. An MPR system allows controlling hydraulic actuators with a pressure control logic, as opposed to traditional flow control logics which are well established in mobile applications. A proper setting of the number of pressure rails and the pressure level in each rail permits optimizing the system efficiency of an MPR solution. MPR systems have been recently proposed with different purposes, but seldom adopt variable pressure control strategies to manage the pressure variations in the rails. Moreover, the related past work attempts to maximize the energy recovery during overrunning load conditions. This paper, instead, addresses the case of an MPR system that handles a multitude of hydraulic functions without significant overrunning loads, as it occurs in many mobile applications such as in agriculture. The paper first presents the design alternatives considered for the MPR system, including the supply, the control valves, and the controller. Second, a dedicated test rig is developed to support experimental activities and validations of a lumped parameter model developed within this research to develop the MPR control strategy. The results describe the operating features of the proposed MPR solution. By considering the case of a traditional LS system, the paper results also show how system efficiency can be highly improved, with energy savings in the order of 49%

24 POWER TRANSMISSION AND DISTRIBUTION↗

Integrating Multi-Source Data for Bi-Level Traffic Simulator Calibration: A Literature Review and Highway Case Study

Traffic simulation serves as a powerful tool for pre-evaluating policies and technologies. In this context, simulation-based Dynamic traffic assignment (DTA) models are capable of capturing traffic dynamics. They are well-known as critical tools in controlling and predicting traffic situations. The reliability of simulation results heavily depends on the calibration process. Most studies in the literature formulate and calibrate simulators based on a single source of collected data or multiple data sets with the same spatiotemporal characteristics. However, in practice, traffic data is collected by various tools with usually different spatial and temporal resolutions. This study introduces a novel approach to taking into account diverse input data from a variety of sources. An iterative bi-level solution is proposed. to equally treat traffic flow and speed data. The upper level solves flow calibration with the exact solution method, and the lower level calibrates the speed with the simultaneous perturbation stochastic approximation (SPSA) algorithm. Subsequently, the effectiveness of the proposed model is investigated using data from a six-mile section of Nashville's I-24 highway in Tennessee. The results demonstrate that our proposed model creates an effective feedback loop between the optimizer and the simulator for calibrating flow and speed to reduce the error between simulated and real data.

42 ENGINEERING↗

A Short-Term Solar Forecasting Platform Using a Physics-Based Smart Persistence Model and Data Imputation Method

Electrical energy plays vital role in our socio-economic activity and therefore ensuring the reliability of the electric grid, from the generation, transmission and distribution level is critical. In order to maintain the power system parameter viz., frequency, voltage, etc., optimally, balancing of generation and consumption is very much essential. However, solar energy is infirm power by nature this is due to cloud cover / other local phenomena. Hence, Photovoltaic (PV) power generation brings a significant challenge to the grid operator due to the variability of the solar energy. The complexity of this challenge in terms of planning and dispatch ability of PV resources, aggravates with the high penetration of solar energy into the electric grid. In this setting, reliable solar radiation forecasting models based on accurate and quality input data become essential. In order to develop a suitable model for predicting solar radiation, quality historical / real time measurement is also needed. Under this study NIWE and NREL jointly developed / tested short-term solar forecasting frameworks using a smart persistence and physics-based smart persistence models for intra-hour forecasting of solar radiation (PSPI) and benchmarked 9 different data imputation techniques in 15 Solar Radiation Resource Assessment (SRRA) stations, located at different parts of India. During any measurement campaign, due to various technical reasons, we may miss few observations. However, the missing observation often reduce the performance of any forecasting model. Therefore, suitable data imputation method would assist us to obtain continuous observation of solar radiation. A station-by-station and method-by-method analysis was carried out to understand the performance of each model. Based on our analysis, among all the data imputation methods, the Kalman data imputation method is better for Indian Weather condition. In addition, Kalman StructTS, Linear, Stine and Arima methods yield slightly inferior accuracy compared to Kalman, but outperform the other methods. The extended solar radiation data are used by solar forecasting models to provide the prediction of solar radiation at 15 SRRA stations. As far as short term forecasting model is concerned, the PSPI model outperforms the Smart Persistence model. However, the forecast error is increases with the forecasting horizon.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOME: Directional medical embedding vectors from Electronic Health Records

Motivation: The increasing availability of Electronic Health Record (EHR) systems has created enormous potential for translational research. Recent developments in representation learning techniques have led to effective large-scale representations of EHR concepts along with knowledge graphs that empower downstream EHR studies. However, most existing methods require training with patient-level data, limiting their abilities to expand the training with multi-institutional EHR data. On the other hand, scalable approaches that only require summary-level data do not incorporate temporal dependencies between concepts. Methods: We introduce a DirectiOnal Medical Embedding (DOME) algorithm to encode temporally directional relationships between medical concepts, using summary-level EHR data. Specifically, DOME first aggregates patient-level EHR data into an asymmetric co-occurrence matrix. Then it computes two Positive Pointwise Mutual Information (PPMI) matrices to correspondingly encode the pairwise prior and posterior dependencies between medical concepts. Following that, a joint matrix factorization is performed on the two PPMI matrices, which results in three vectors for each concept: a semantic embedding and two directional context embeddings. They collectively provide a comprehensive depiction of the temporal relationship between EHR concepts. Results: We highlight the advantages and translational potential of DOME through three sets of validation studies. First, DOME consistently improves existing direction-agnostic embedding vectors for disease risk prediction in several diseases, for example achieving a relative gain of 5.5% in the area under the receiver operating characteristic (AUROC) for lung cancer. Second, DOME excels in directional drug-disease relationship inference by successfully differentiating between drug side effects and indications, correspondingly achieving relative AUROC gain over the state-of-the-art methods by 10.8% and 6.6%. Finally, DOME effectively constructs directional knowledge graphs, which distinguish disease risk factors from comorbidities, thereby revealing disease progression trajectories. The source codes are provided at https://github.com/celehs/Directional-EHRembedding.

60 APPLIED LIFE SCIENCES↗

OzDES reverberation mapping program: Lag recovery reliability for 6-yr C iv analysis

ABSTRACT We present the statistical methods that have been developed to analyse the OzDES reverberation mapping sample. To perform this statistical analysis we have created a suite of customizable simulations that mimic the characteristics of each source in the OzDES sample. These characteristics include: the variability in the photometric and spectroscopic light curves, the measurement uncertainties, and the observational cadence. By simulating the sources in the OzDES sample that contain the C iv emission line, we developed a set of criteria that rank the reliability of a recovered time-lag depending on the agreement between different recovery methods, the magnitude of the uncertainties, and the rate at which false positives were found in the simulations. These criteria were applied to simulated light curves and these results used to estimate the quality of the resulting Radius–Luminosity relation. We grade the results using three quality levels (gold, silver, and bronze). The input slope of the R–L relation was recovered within 1σ for each of the three quality samples, with the gold standard having the lowest dispersion with a recovered a R–L relation slope of 0.454 ± 0.016 with an input slope of 0.47. Future work will apply these methods to the entire OzDES sample of 771 AGN.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of tomographic redshift bin width errors on cosmological probes

ABSTRACT Systematic errors in the galaxy redshift distribution n(z) can propagate to systematic errors in the derived cosmology. We characterize how the degenerate effects in tomographic bin widths and galaxy bias impart systematic errors on cosmology inference using observational data from the Deep Lens Survey. For this we use a combination of galaxy clustering and galaxy–galaxy lensing. We present two end-to-end analyses from the catalogue level to parameter estimation. We produce an initial cosmological inference using fiducial tomographic redshift bins derived from photometric redshifts, then compare this with a result where the redshift bins are empirically corrected using a set of spectroscopic redshifts. We find that the derived parameter S8 ≡ σ8(Ωm/.3)1/2 decreases from $0.841^{+0.062}_{-0.061}$ to $0.781^{0.061}_{0.054}$ upon correcting the n(z) errors in the second method.

79 ASTRONOMY AND ASTROPHYSICS↗

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗

Generative adversarial networks for scintillation signal simulation in EXO-200

Generative Adversarial Networks trained on samples of simulated or actual events have been proposed as a way of generating large simulated datasets at a reduced computational cost. In this work, a novel approach to perform the simulation of photodetector signals from the time projection chamber of the EXO-200 experiment is demonstrated. The method is based on a Wasserstein Generative Adversarial Network — a deep learning technique allowing for implicit non-parametric estimation of the population distribution for a given set of objects. Our network is trained on real calibration data using raw scintillation waveforms as input. We find that it is able to produce high-quality simulated waveforms an order of magnitude faster than the traditional simulation approach and, importantly, generalize from the training sample and discern salient high-level features of the data. In particular, the network correctly deduces position dependency of scintillation light response in the detector and correctly recognizes dead photodetector channels. Furthermore, the network output is then integrated into the EXO-200 analysis framework to show that the standard EXO-200 reconstruction routine processes the simulated waveforms to produce energy distributions comparable to that of real waveforms. Finally, the remaining discrepancies and potential ways to improve the approach further are highlighted.

47 OTHER INSTRUMENTATION↗

RT-EOM-CCSD Calculations of Inner and Outer Valence Ionization Energies and Spectral Functions

Photoelectron spectroscopy (PES) is a standard experimental method for material characterization, but its interpretation can be hampered by its reliance on standard materials. To facilitate the study of unknown systems, theoretical methods are desirable. Here we present a real-time equation-of-motion coupled cluster (RT-EOM-CC) approach for valence PES, extending our core-level development. Here, we demonstrate that RT-EOM-CC yields ionization energies and spectral functions in good agreement with experiment and CI-based methods, even for some more correlated cases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Residential Building Energy Efficiency Field Studies: Low-Rise Multifamily

In recent years, the U.S. Department of Energy (DOE) has conducted a series of research studies to validate energy efficient building technologies in the field. Much of the work has focused on single-family construction, and some has also addressed commercial energy codes. The work detailed in this DOE-funded study (EE0007616) focuses on low-rise multifamily buildings (three stories or fewer above grade) in various regions of the United States, and reports on how state-level building codes are being implemented, both in terms of observed characteristics and also in terms of estimated energy impacts. Nearly 100 buildings across four states—Illinois, Minnesota, Oregon, and Washington—were sampled, which represent a range of climate types from mild temperature to very cold continental. Both common entry and outdoor entry buildings were included, and a parallel research project evaluated envelope air tightness and current still-evolving air tightness testing methods. Finally, a set of structured interviews of building designers and other relevant professionals was carried to out to gain more insight into this market. To the greatest extent possible, the methodology developed under the project for low-rise multifamily buildings mirrored the approach established by Pacific Northwest National Laboratory (PNNL) for single-family residential buildings (https://www.energy.gov/eere/buildings/downloads/residential-building-energy-code-field-study). This included the general approach to sampling, recruitment, and data collection, as well as data analysis and presentation. The range of permitting dates for the sites encompassed two energy code cycles in most regions. All states in the study had adopted a variation of the International Energy Conservation Code (IECC) for the structure of their state code. The low-rise multifamily occupancy presents a hybrid building type: most of the building’s conditioned floor area was covered by the residential chapter of the code while portions of the building (such as corridors and common spaces) fell under the commercial code chapter. The key items assessed in this work were: Building Shell—exterior wall insulation, ceiling insulation, foundation insulation, windows. Common Areas—HVAC and lighting. Living Units—lighting, ventilation. A few items were not assessed in detail, given their relative paucity in this occupancy type; these included duct leakage, pipe insulation, and hot water circulation controls. Building characteristics were collected via a combination of architectural, mechanical, electrical, and plumbing plan reviews and field inspections, and entered into a spreadsheet-based tool that was later queried to build a database. Data went through quality control both upon arrival and via a later semi-automated review and assurance process. Most of the data are presented graphically so that the reader can quickly assess compliance with the applicable energy codes (both by state and by code year). As a final step, EnergyPlus™ simulations were created for all buildings in the study to estimate both the as-found energy use intensity (EUI) and the energy and CO 2 that could be saved if features that were found to not meet code minimums were brought up to code. The savings estimates were tabulated for each of the four states in the study. The research team found that the single-family approach was largely applicable to low-rise multifamily buildings. This applies to both the data collection and the prototype EUI analysis. Most of the occupied space is living units and falls under residential energy codes, and many characteristics use similar envelope construction and relatively straightforward mechanical systems and lighting. One of the most challenging aspects of this work was to build an effective spreadsheet-based data collection instrument that could allow efficient collection of both building plan and field data. The research team is of the view that other methods could be equally effective if the work is done carefully with diligent quality control. The primary findings for the work center around the thermal envelope and mechanical systems and lighting at the sites: For thermal envelope components, the majority of buildings met or were better than the prescriptive code.This suggests that building designers and builders are aware of code requirements. In some cases, surveyed buildings were designed to qualify for energy efficiency certification programs. These buildings made up at least 20% of sampled buildings in each state. Almost all buildings met mechanical system efficiency requirements (for both living units and common areas). In some cases, sites employed systems that were considerably more efficient than required by the applicable energy code. Dwelling units had a majority of high-efficacy lighting, often in excess of the state’s residential code requirements. While high-efficacy fixtures were also typical in common areas (corridors and stairwells), lighting power densities (LPDs) in these areas were sometimes higher than levels dictated by the applicable part of the state commercial energy code. The simulation models run on a series of low-rise multifamily prototypes, informed by a composite of the field data collected, calculated annual EUIs of between 20 and 50 kBtu/ft2-yr, with the range representing the effects of both building characteristics and building location (climate zone). A detailed process (based on simulations of prototype buildings) was used to estimate the amount of avoided energy use that would occur if 100% adherence to energy codes were attained. The results indicated modest savings are attainable for items such as window thermal performance and common area lighting. The result is overall only a modest potential for additional energy savings, averaging about 10% of EUI.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Size Effect of Local Current-Voltage Characteristics of MX 2 Nanoflakes: Local Density of States Reconstruction from Scanning Tunneling Microscopy Experiments

Local current-voltage characteristics for low-dimensional transition-metal dichalcogenides (LDTMD), as well as the reconstruction of their local density of states (LDOS) from scanning tunneling microscopy (STM) experiments, are of fundamental interest and can be useful for advanced applications. Most of the existing models either have limited applicability for complex-shaped LDTMDs (e.g., those based on the Simmons approach) or require solving of an ill-defined integral equation to deconvolute the unknown LDOS (e.g., those based on the Tersoff approach). Using a serial expansion of the Tersoff formulas, we propose a flexible method to reconstruct the LDOS from local current-voltage characteristics measured in STM experiments. We establish a set of key physical parameters, which characterize the tunneling current of a STM-probe–sample contact and the sample LDOS expanded in Gaussian functions. Using a direct variational method coupled with probabilistic analysis, we determine these parameters from the STM experiments for MoS2 nanoflakes with different numbers of layers. The main result is the reconstruction of the LDOS in a relatively wide energy range around the Fermi level, which allows us to gain insight into the local band structure of LDTMDs. The reconstructed LDOS reveal pronounced size effects for the single-layer, two-layer, and three-layer MoS 2 nanoflakes, which we relate to the low dimensionality and strong bending or corrugation of the nanoflakes. We hope that the proposed elaboration of the Tersoff approach, allowing LDOS reconstruction, will be of critical interest for the quantitative description of STM experiments and also of use to better understand the microscopic physical aspects of the surface, strain, and bending contributions to the LDTMDs’ electronic properties.

42 ENGINEERING↗