Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling bias”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Assessing the Ability of Instantaneous Aircraft and Sonde Measurements to Characterize Climatological Means and Long-Term Trends in Tropospheric Composition

Over four decades of measurements exist that sample the 3-D composition of reactive trace gases in the troposphere from approximately weekly ozone sondes, instrumentation on civil aircraft, and individual comprehensive aircraft field campaigns. An obstacle to using these data to evaluate coupled chemistry-climate models (CCMs)the models used to project future changes in atmospheric composition and climateis that exact space-time matching between model fields and observations cannot be done, as CCMs generate their own meteorology. Evaluation typically involves averaging over large spatiotemporal regions, which may not reflect a true average due to limited or biased sampling. This averaging approach generally loses information regarding specific processes. Here we aim to identify where discrete sampling may be indicative of long-term mean conditions, using the GEOS-Chem global chemical-transport model (CTM) driven by the MERRA reanalysis to reflect historical meteorology from 2003 to 2012 at 2o by 2.5o resolution. The model has been sampled at the time and location of every ozone sonde profile available from the Would Ozone and Ultraviolet Radiation Data Centre (WOUDC), along the flight tracks of the IAGOSMOZAICCARABIC civil aircraft campaigns, as well as those from over 20 individual field campaigns performed by NASA, NOAA, DOE, NSF, NERC (UK), and DLR (Germany) during the simulation period. Focusing on ozone, carbon monoxide and reactive nitrogen species, we assess where aggregates of the in situ data are representative of the decadal mean vertical, spatial and temporal distributions that would be appropriate for evaluating CCMs. Next, we identically sample a series of parallel sensitivity simulations in which individual emission sources (e.g., lightning, biogenic VOCs, wildfires, US anthropogenic) have been removed one by one, to assess where and when the aggregated observations may offer constraints on these processes within CCMs. Lastly, we show results of an additional 31-year simulation from 1980-2010 of GEOS-Chem driven by the MACCity emissions inventory and MERRA reanalysis at 4o by 5o. We sample the model at every WOUDC sonde and flight track from MOZAIC and NASA field campaigns to evaluate which aggregate observations are statistically reflective of long-term trends over the period.

Atmospheric composition↗
DOE OSTI2021

Understanding collective human movement dynamics during large-scale events using big geosocial data analytics

Conventional approaches for modeling human mobility pattern often focus on human activity and movement dynamics in their regular daily lives and cannot capture changes in human movement dynamics in response to large-scale events. With the rapid advancement of information and communication technologies, many researchers have adopted alternative data sources (e.g., cell phone records, GPS trajectory data) from private data vendors to study human movement dynamics in response to large-scale natural or societal events. Big geosocial data such as georeferenced tweets are publicly available and dynamically evolving as real-world events are happening, making it more likely to capture the real-time sentiments and responses of populations. However, precisely-geolocated geosocial data is scarce and biased toward urban population centers. In this research, we developed a big geosocial data analytical framework for extracting human movement dynamics in response to large-scale events from publicly available georeferenced tweets. The framework includes a two-stage data collection module that collects data in a more targeted fashion in order to mitigate the data scarcity issue of georeferenced tweets; in addition, a variable bandwidth kernel density estimation(VB-KDE) approach was adopted to fuse georeference information at different spatial scales, further augmenting the signals of human movement dynamics contained in georeferenced tweets. To correct for the sampling bias of georeferenced tweets, we adjusted the number of tweets for different spatial units (e.g., county, state) by population. To demonstrate the performance of the proposed analytic framework, we chose an astronomical event that occurred nationwide across the United States, i.e., the 2017 Great American Eclipse, as an example event and studied the human movement dynamics in response to this event. Finally, this analytic framework can easily be applied to other types of large-scale events such as hurricanes or earthquakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Bayesian model for multivariate discrete data using spatial and expert information with application to inferring building attributes

When modeling sparsely observed multivariate data, strong prior information elicited from experts can be used to bolster predictive accuracy and counteract sampling bias. Similarly, modeling autocorrelation in space can help make use of co-occurrence patterns present in many types of spatial data. To make use of both expert prior information and spatial structure, we propose a novel graphical model for a spatial Bayesian network developed specifically to address challenges in inferring the attributes of buildings from geographically sparse observational data. This model is implemented as the sum of a spatial multivariate Gaussian random field and a tabular conditional probability function in real-valued space prior to projection onto the probability simplex. This modeling form is especially suitable for the usage of prior information in the form of sets of atomic rules obtained from experts. To perform inference with missing data, we implement a Markov chain Monte Carlo scheme composed of alternating steps of Gibbs sampling of missing entries and Hamiltonian Monte Carlo for model parameters. A case study in building attribution is presented to highlight the advantages and limitations of this approach.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗
NASA NTRS2011

Combining Satellite and in Situ Data with Models to Support Climate Data Records in Ocean Biology

The satellite ocean color data record spans multiple decades and, like most long-term satellite observations of the Earth, comes from many sensors. Unfortunately, global and regional chlorophyll estimates from the overlapping missions show substantial biases, limiting their use in combination to construct consistent data records. SeaWiFS and MODIS-Aqua differed by 13% globally in overlapping time segments, 2003-2007. For perspective, the maximum change in annual means over the entire Sea WiFS mission era was about 3%, and this included an El NinoLa Nina transition. These discrepancies lead to different estimates of trends depending upon whether one uses SeaWiFS alone for the 1998-2007 (no significant change), or whether MODIS is substituted for the 2003-2007 period (18% decline, P less than 0.05). Understanding the effects of climate change on the global oceans is difficult if different satellite data sets cannot be brought into conformity. The differences arise from two causes: 1) different sensors see chlorophyll differently, and 2) different sensors see different chlorophyll. In the first case, differences in sensor band locations, bandwidths, sensitivity, and time of observation lead to different estimates of chlorophyll even from the same location and day. In the second, differences in orbit and sensitivities to aerosols lead to sampling differences. A new approach to ocean color using in situ data from the public archives forces different satellite data to agree to within interannual variability. The global difference between Sea WiFS and MODIS is 0.6% for 2003-2007 using this approach. It also produces a trend using the combination of SeaWiFS and MODIS that agrees with SeaWiFS alone for 1998-2007. This is a major step to reducing errors produced by the first cause, sensor-related discrepancies. For differences that arise from sampling, data assimilation is applied. The underlying geographically complete fields derived from a free-running model is unaffected by solar zenith angle requirements and obscuration from clouds and aerosols. Combined with in situ dataenhanced satellite data, the model is forced into consistency using data assimilation. This approach eliminates sampling discrepancies from satellites. Combining the reduced differences of satellite data sets using in situ data, and the removal of sampling biases using data assimilation, we generate consistent data records of ocean color. These data records can support investigations of long-term effects of climate change on ocean biology over multiple satellites, and can improve the consistency of future satellite data sets.

Gregg, Watson↗
NASA NTRS1994

On the effects of filtering on convective-core statistics

Analysis of convective events north of Australia during the equatorial Monsoon Experiment (EMEX) reinforces accumulating evidence that convective cores over the tropical ocean are weak, with small diameters and vertical velocities. Superposition of smaller-scale turbulence on the convective core signal could yield artificially small cores, and thus produce sampling bias by aircraft. To investigate this potential source of bias, a side-by-side analysis of filtered and unfiltered vertical velocity data was performed while investigating the EMEX convective cores. Data were collected by the NOAA P3 and NCAR Electra using filtered (Graham filter) and non-filtered cores. The effects of filtering vertical velocity data to identify updraft and downdraft cores were noticeable but smaller than expected. Updraft cores had a vertical velocity greater than 1 m/sec for at least 500 m, and downdraft cores were defined analogously so that the filter eliminated events with less than a 500 m apparent diameter. Median and 10 percent core diameters were increased by 20-25 percent at most, with larger increases at higher levels (fewer small events). The maximum vertical velocity was changed only by changing the individual events in the distribution, but changes were detectable. At the lower levels, where many small, weak cores are introduced by filtering, average vertical velocity for the updraft cores was lowered by up to around 10 percent. Mass flux was changed the least, since diameter and mean vertical velocity respond to filtering in opposite ways. The filtering procedure increased total mass flux slightly due to inclusion of more upward moving air into the sample.

Lemone, Margaret A.↗
NASA NTRS2016

Water in the Cratonic Mantle Lithosphere

The fact that Archean and Proterozoic cratons are underlain by the thickest (>200 km) lithosphere on Earth has always puzzled scientists because the dynamic convection of the surrounding asthenosphere would be expected to delaminate and erode these mantle lithospheric "keels" over time. Although density and temperature of the cratonic lithosphere certainly play a role in its strength and longevity, the role of water has only been recently addressed with data on actual mantle samples. Water in mantle lithologies (primarily peridotites and pyroxenites) is mainly stored in nominally anhydrous minerals (olivine, pyroxene, garnet) where it is incorporated as hydrogen bonded to structural oxygen in lattice defects. The property of hydrolytic weakening of olivine [4] has generated the hypothesis that olivine, the main mineral of the upper mantle, may be dehydrated in cratonic mantle lithospheres, contributing to its strength. This presentation will review the distribution of water concentrations in four cratonic lithospheres. The distribution of water contents in olivine from peridotite xenoliths found in kimberlites is different in each craton (Figure 1). The range of water contents of olivine, pyroxene and garnet at each xenolith location appears linked to local metasomatic events, some of which occurred later then the Archean and Proterozoic when these peridotites initially formed via melting. Although the low olivine water contents (<10 ppm wt H2O) at > 6 GPa at the base of the Kaapvaal cratonic lithosphere may contribute to its strength, and prevent its delamination, the wide range of those from Siberian xenoliths is not compatible with providing a high enough viscosity contrast with the asthenophere. The water content in olivine inclusions from Siberian diamonds, on the other hand, have systematically low water contents (<20 ppm wt H2O). The xenoliths may represent a biased sample of the cratonic lithosphere with an over-­abundance of metasomatized peridotites with high water contents. The olivine inclusions, however, may have been protected from metasomatism by their host diamond and record the overall low olivine water content of the cratonic lithosphere. Water may thus still play a role in cratonic keel longevity.

Peslier, A. H.↗
DOE OSTI2024

Adaptive Scalpel Scanning Probe Microscopy for Enhanced Volumetric Sensing in Tomographic Analysis

Controlling nanoscale tip‐induced material removal is crucial for achieving atomic‐level precision in tomographic sensing with atomic force microscopy (AFM). While advances have enabled volumetric probing of conductive features with nanometer accuracy in solid‐state devices, materials, and photovoltaics, limitations in spatial resolution and volumetric sensitivity persist. This work identifies and addresses in‐plane and vertical tip‐sample junction leakage as sources of parasitic contrast in tomographic AFM, hindering real‐space 3D reconstructions. Novel strategies are proposed to overcome these limitations. First, the contrast mechanisms analyzing nanosized conductive features are explored when confining current collection purely to in‐plane transport, thus allowing reconstruction with a reduction in the overestimation of the lateral dimensions. Furthermore, an adaptive tip‐sample biasing scheme is demonstrated for the mitigation of a class of artefacts induced by the high electric field inside the thin oxide when volumetrically reduced. This significantly enhances vertical sensitivity by approaching the intrinsic limits set by quantum tunneling processes, allowing detailed depth analysis in thin dielectrics. The effectiveness of these methods is showcased in tomographic reconstructions of conductive filaments in valence change memory, highlighting the potential for application in nanoelectronics devices and bulk materials and unlocking new limits for tomographic AFM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2022

Molecular Driving Force for Facet Selectivity of Sequence-Defined Amphiphilic Peptoids at Au–Water Interfaces

Shape-controlled colloidal nanocrystal syntheses often require facet-selective solution-phase chemical additives to regulate surface free energy, atom addition/migration fluxes or particle attachment rates. Because of their highly tunable properties and robustness to a wide range of experimental conditions, peptoids represent a very promising class of next generation functional additives for control over nanocrystal growth. However, understanding the origin of facet selectivity at the molecular level is critical to generalizing their design. In this work, we employ molecular dynamics simulations and biased sampling methods and report stronger selectivity to Au(111) than to Au(100) for Nce3Ncp6, a peptoid that has been shown to assist the formation of five-fold twinned Au nanostars. We find that facet-selectivity is achieved through synergistic effects of both peptoid-surface and solvent-surface interactions. Moreover, the amphiphilic nature of Nce3Ncp6 together with the order of peptoid-peptoid and peptoid-surface binding energies, i.e., peptoid-Au(100) < peptoid-peptoid < peptoid-Au(111), further amplifies its distinct collective behavior on different Au surfaces. Our studies provide a fundamental understanding of the molecular origin of facet-selective adsorption and highlight the possibility of future designs and uses of sequence-defined peptoids for predictive syntheses of nanocrystals with designed shapes and properties.

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

Defects in Highly Anisotropic Transition-Metal Dichalcogenide PdSe 2

The atomic and electronic structures of the pristine PdSe 2 as well as various intrinsic vacancy defects in PdSe 2 are studied comprehensively by combining scanning tunneling microscopy, spectroscopy and density functional theory calculations. Other than the topmost Se atoms, sublayer Pd atoms, the intrinsic Pd and Se vacancy-defects are identified. Both V Se and V Pd defects induce defect states near Fermi level. As a result, the vacancy defects can be negatively charged by tip gating effect. At negative sample bias, the screened Coulomb interaction between STM tip and the charged vacancies creates disk-like protrusion around V Pd and crater-like features around V Se . The magnification effect of the long-range charge localization at the charged defect site makes sublayer defects as deep as 1nm visible even in STM images. In conclusion, this result proves that by gating the probe, scanning probe microscopy can be used as an easy tool for characterizing sublayer defects in a nondestructive way.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Avoided Level Crossing and Entangled States of Interacting Hydrogen Molecules Detected by the Quantum Superposition Microscope

Pump–probe measurements by ultrashort THz pulses can be used to excite and follow the coherence dynamics in the time domain of single hydrogen molecules (H 2 ) in the junction of a scanning tunneling microscope (STM). By tailoring the resonance frequency through the sample bias, we identified two spectral signatures of the interactions among multiple H 2 molecules. First, the avoided level crossing featured by energy gaps ranging from 20 to 80 GHz was observed because of the level repulsion between two H 2 molecules. Second, the tip can sense the signal of H 2 outside the junction through the projective measurement on the H 2 inside the junction, owing to the entangled states created through the interactions. A dipolar-type interaction was integrated into the tunneling two-level system model of H 2 , enabling accurate reproduction of the observed behaviors. Furthermore, our results obtained by the quantum superposition microscope reveal the intricate quantum mechanical interplay among H 2 molecules and additionally provide a 2D platform to investigate unresolved questions of amorphous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Automated Construction of Artificial Lattice Structures with Designer Electronic States

Manipulating matter with a scanning tunneling microscope (STM) enables the creation of atomically defined artificial structures that host designer quantum states. However, the time-consuming nature of the manipulation process, coupled with the sensitivity of the STM tip, constrains the exploration of diverse configurations and limits the size of the designed features. In this study, we present a reinforcement learning (RL)-based framework for creating artificial structures by spatially manipulating carbon monoxide (CO) molecules on a copper substrate by using the STM tip. The automated workflow combines molecule detection and manipulation, employing deep-learning-based object detection to locate CO molecules and linear assignment algorithms to allocate these molecules to designated target sites. We initially perform molecule maneuvering based on randomized parameter sampling for sample bias, tunneling current set point, and manipulation speed. This data set is then structured into an action trajectory used to train an RL agent. The model is subsequently deployed on the STM for real-time fine-tuning of the manipulation parameters during structure construction. Our approach incorporates path-planning protocols coupled with active drift compensation to enable atomically precise fabrication of structures with significantly reduced human input while realizing larger-scale artificial lattices with the desired electronic properties. Furthermore, using our approach, we demonstrate the automated construction of an extended artificial graphene lattice and confirm the existence of a characteristic Dirac point in its electronic structure. Further challenges regarding the RL-based structural assembly scalability are discussed.

Algorithms↗
DOE OSTI2020

Why Are There so Few Reports of High-Energy Electron Drift Resonances? Role of Radial Phase Space Density Gradients

Models of monochromatic Pc5 (2–7 mHz) ultralow frequency (ULF) wave interactions with high energy (greater than ~1 MeV) electrons predict drift resonant interactions that can cause rapid radial transport and acceleration. There are few reports of electron drift resonance at energies greater than ~1 MeV, in contrast to lower energies; moreover, all previous reports occur in the aftermath of interplanetary shocks. These two facts are difficult to reconcile with theory and numerical simulations predicting that greater than ~1 MeV drift resonances should occur more often and in a wider variety of driving conditions. In this study, we show that a combination of observational sampling biases and nominal radial phase space density gradients is one explanation for this discrepancy between theory and observations. In particular, we examine electron dynamics in two case studies with very similar satellite coverage, solar wind conditions, and Pc5 wave properties, yet with different radial phase space density profiles. Using global wave and particle observations, we show that the events have vastly different particle responses despite having similar wave properties. Placing these results in context with past studies, we further show here that nominal radial PSD gradients near geostationary orbit can mask the expected drift resonance particle response and explain (1) the small number of past greater than ~1 MeV drift resonance reports and (2) the restriction of these reports to interplanetary shock events. We argue that future observational studies characterizing radial transport via drift resonance should examine global particle dynamics, including observations of the radial phase space density profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Global terrestrial nitrogen fixation and its modification by agriculture

Biological nitrogen fixation (BNF) is the largest natural source of new nitrogen (N) that supports terrestrial productivity, yet estimates of global terrestrial BNF remain highly uncertain. Here, in this study, we show that this uncertainty is partly because of sampling bias, as field BNF measurements in natural terrestrial ecosystems occur where N fixers are 17 times more prevalent than their mean abundances worldwide. To correct this bias, we develop new estimates of global terrestrial BNF by upscaling field BNF measurements using spatially explicit abundances of all major biogeochemical N-fixing niches. We find that natural biomes sustain lower BNF, 65 (52–77) Tg N yr −1 , than previous empirical bottom-up estimates, with most BNF occurring in tropical forests and drylands. We also find high agricultural BNF in croplands and cultivated pastures, 56 (54–58) Tg N yr −1 . Agricultural BNF has increased terrestrial BNF by 64% and total terrestrial N inputs from all sources by 60% over pre-industrial levels. Our results indicate that BNF may impose stronger constraints on the carbon sink in natural terrestrial biomes and represent a larger source of agricultural N than is generally considered in analyses of the global N cycle, with implications for proposed safe operating limits for N use.

Reis Ely, Carla R. [Oregon State Univ., Corvallis,↗
DOE OSTI2023

Identifying the secondary electron cutoff in ultraviolet photoemission spectra for work function measurements of non-ideal surfaces

Abstract Absolute values of work functions can be determined in ultraviolet photoemission spectroscopy (UPS) by measuring the minimum kinetic energy of secondary electrons generated by a known photon energy. However, some samples can produce spectra that are difficult to interpret due to additional intensity below the true secondary electron cutoff. Disordered absorbates on elemental metals add small intensity below the onset for the transition metal surfaces studied, which can be attributed to energy losses after photoelectrons are generated. In contrast, spectra from WO 3−x films can produce multiple onsets with comparable intensity which do not fit this model. False onsets (in the context of work function measurements) can be minimized by optimizing experimental detection parameters including limiting analyzer acceptance angles and pass energy. True work functions can be identified by examining the onsets as the sample bias is varied.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗