Lattice Distortion and Low-Frequency Anharmonic Phonons Suppress Charge Recombination in Lead Halide Perovskites upon Pseudohalide Doping: Time-Domain Ab Initio Analysis
Not Available
SEARCH · Engineering Papers
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.
Not Available
In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.
A Programmable Logic Controller (PLC) is a microprocessor-based controller that is used to automate physical processes in critical infrastructure and various other industries and manufacturing sectors. Initially, PLCs were completely isolated from the Internet, and cyber security was not incorporated at the time of development. The introduction of industry 4.0 and the evolution of ICS systems to communicate over public IP addresses from the Internet enhanced productivity and efficiency, but Internet connectivity exposed the systems and their vulnerabilities, which led to an increase in cyber attacks. When a system is sabotaged/compromised, security analysts need to get to the root cause of the attack as quickly as possible to recover the system. To do so, memory forensic analysis is critical to provide a unique insight into the run-time memory activities and extract a reliable source of evidence. In this paper, we analyze the memory structure of the Schneider Electric Modicon M221 PLC. To build a memory profile, we reverse engineer the communication protocol and conduct differential analysis to gain knowledge about the structure of the memory and the low-level representation of control logic instructions. We then identify dynamic and static memory regions by modifying different project fields and conducting differential analysis, which allows us to identify boundaries of critical memory structures and extract important forensic artifacts that can be found in the memory. The Python implementation of the memory profile can help reduce the time and effort required for manual analysis in case of cyber incident or system failure.
This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.
The importance of high intensity few- to single-cycle laser pulses for applications such as intense isolated attosecond pulse generation is constantly growing, and with the breakdown of the monochromatic approximation in field ionization models, the few-cycle pulse (FCP) interaction with solids near the damage threshold has ushered a new paradigm of nonperturbative light–matter interaction. In this Letter, we systematically study and contrast how femtosecond laser-induced damage and ablation behaviors of SiO 2 /HfO 2 -based reflective multilayer dielectric thin film systems vary between FCP and 110 fs pulses. With time-resolved surface microscopy and ex situ analysis, we show that there are distinct differences in the interaction depending on the pulse duration, specifically in the “blister” morphology formation at lower fluences (damage) as well as in the dynamics of debris formation at higher fluences (ablation).
Abstract X-ray free electron lasers (XFELs) have emerged as powerful sources of short and intense x-ray pulses. We propose a simple and robust procedure which takes advantage of the inherent stochasticity of self-amplified stimulated emission (SASE) pulses to enhance the time-resolution and signal strength of the recorded data. Notably, the proposed method is able to enhance the average signal without knowledge of the signal strength of individual shots. Simple metrics for the probe pulses are introduced, such as an effective pulse duration applicable to SASE pulses characterised in the time domain using e.g. an X-band transverse cavity. The approach is evaluated using simulated and real pulse data in the context of ultrafast electron dynamics in a molecule. Utilising H 2 as a model system, we demonstrate the efficacy of the method theoretically, successfully enhancing the predicted nonresonant ultrafast x-ray scattering signal associated with electron dynamics. The method presented is broadly applicable and offers a general strategy for enhancing time-dependent observables at XFELs.
Polyurethane (PU) thermosets are usually landfilled at the end of their service lifetimes and cannot be recycled through conventional means. Although dibutyltin dilaurate (DBTDL) is an effective catalyst for reprocessing PU thermosets, organotins are immunotoxic and teratogenic, which is concerning for catalysts embedded within covalent adaptable networks (CANs). Here, we identify Zr(acac) 4 and Zr(tmhd) 4 as catalysts to reprocess PU CANs more sustainably. Thermoset PU foams containing Zr-based catalysts exhibit characteristic stress relaxation times between 6 and 200 s, similar or faster than DBTDL. These catalysts are capable of reprocessing a PU foam for four or five cycles. Dynamic mechanical thermal analysis indicates that zirconium-based catalysts preserve glass transition temperature and crosslink density between cycles and are promising catalysts for bulk reprocessing of thermoset PUs. Here, a solvent-free reprocessing method was developed by co-milling the PU foam and the Zr catalyst. Ultimately, these reprocessing methods using Zr-based catalysts will impart greater sustainability and circularity to PU plastics.
Increasing interest in sustainable chemistry coupled with the quest to explore new reactivity has spurred research on first-row transition metal complexes for potential applications in a variety of settings. One of the more active areas of research is photoredox catalysis, where the synthetically tunable nature of their electronic structures provides a rich palette of options for tailoring their reactivity to a desired chemical transformation. Understanding the mechanism of excited-state reactivity is critical for the informed development of next-generation catalysts, which in turn requires information concerning the propensity of their electronic excited states to engage in the desired electron or energy transfer processes. Herein we provide direct evidence of the highly oxidizing nature of the lowest-energy ligand-field (LF) excited state of a first-row d 6 -low-spin Co(III) photosensitizer [Co(4,4'-Br 2 bpy) 3 ] 3+ (where 4,4'-Br 2 bpy is 4,4'-dibromo-2,2'-bipyridine). The redox potential associated with the LF excited state of the Co(III) complex was bracketed by performing bimolecular quenching studies by using a series of simple organic electron donors. Time-resolved absorption spectroscopy confirmed a dynamic quenching process attributed to reductive quenching of the lowest-energy ligand-field excited state of the Co(III) chromophore. Analysis of the Stern–Volmer plots for each chromophore-quencher pair revealed a limiting value of E red * ~ 1.25 V vs Fc/Fc + for the metal-centered excited state, which is significantly stronger than that of more commonly employed transition metal-based photoredox agents such as [Ru(bpy) 3 ] 2+ (E red * = 0.32 V vs Fc/Fc + ) and [Ir(ppy) 2 (bpy)] + (E red * = 0.27 V vs Fc/Fc + ). Finally, these results suggest that this class of chromophores could find utility in applications requiring the activation of oxidatively resistant organic substrates for photoredox catalysis.
Abstract Due to its high energy density, lithium metal is a promising electrode for future energy storage. However, its practical capacity, cyclability and safety heavily depend on controlling its reactivity in contact with liquid electrolytes, which leads to the formation of a solid electrolyte interphase (SEI). In particular, there is a lack of fundamental mechanistic understanding of how the electrolyte composition impacts the SEI formation and its governing processes. Here, we present an in-depth model-based analysis of the initial SEI formation on lithium metal in a carbonate-based electrolyte. Thereby we reach for significantly larger length and time scales than comparable molecular dynamic studies. Our multiscale kinetic Monte Carlo/continuum model shows a layered, mostly inorganic SEI consisting of LiF on top of Li 2 CO 3 and Li after 1 µs. Its formation is traced back to a complex interplay of various electrolyte and salt decomposition processes. We further reveal that low local Li + concentrations result in a more mosaic-like, partly organic SEI and that a faster passivation of the lithium metal surface can be achieved by increasing the salt concentration. Based on this we suggest design strategies for SEI on lithium metal and make an important step towards knowledge-driven SEI engineering.
Income and its distribution profile are important determinants of residential energy demand and carry direct implications for human well-being and climate. We explore the sensitivity of residential energy systems to income growth and distribution across shared socioeconomic pathway-representative concentration pathways scenarios using a global, integrated, multisector dynamics model, Global Change Analysis Model, which tracks national/regional household energy services and fuel choice by income decile. Nation/region energy use patterns across deciles tend to converge over time with aggregate income growth, as higher-income consumers approach satiation levels in floorspace and energy services. However, in some regions, existing within-region inequalities in energy consumption persist over time due to slow income growth in lower income groups. Due to continued differences in fuel types, lower income groups will have higher exposure to household air pollution, despite lower contributions to greenhouse gas emissions. We also find that the share of income dedicated to energy is higher for lower deciles, with strong regional differences.
Microbial communities assemble through a complex set of interactions between microbes and their environment, and the resulting metabolic impact on the host ecosystem can be profound. Microbial activity is known to impact human health, plant growth, water quality, and soil carbon storage which has lead to the development of many approaches and products meant to manipulate the microbiome. In order to understand, predict, and improve microbial community engineering, genome-scale modeling techniques have been developed to translate genomic data into inferred microbial dynamics. However, these techniques rely heavily on simulation to draw conclusions which may vary with unknown parameters or initial conditions, rather than more robust qualitative analysis. To better understand microbial community dynamics using genome-scale modeling, we provide a tool to investigate the network of interactions between microbes and environmental metabolites over time. Using our previously developed algorithm for simulating microbial communities from genome-scale metabolic models (GSMs), we infer the set of microbe-metabolite interactions within a microbial community in a particular environment. Because these interactions depend on the available environmental metabolites, we refer to the networks that we infer as metabolically contextualized , and so name our tool MetConSIN: Met abolically Con textualized S pecies I nteraction N etworks.
Identification of atmospheric conditions within a multivariable atmospheric data set is a necessary step in the validation of emerging and existing high-fidelity models used to simulate wind plant flows and operation.Atmospheric conditions relevant for wind energy research include stationary conditions, given the need for well-converged statistics for model validation, as well as conditions observed less frequently, such as extreme atmospheric events, which are used in wind turbine and wind plant design.Aggregation of observations without regard to covariance between time series discounts the dynamical nature of the atmosphere and is not sufficiently representative of atmospheric conditions.Identification and characterization of continuous time periods with atmospheric conditions that have a high value for analysis or simulation set the stage for more advanced model validation and the development of real-time control and operational strategies.The current work explores a single metric for variation in a multivariate data sample that quantifies variability within each channel as well as covariance between channels.The total variation is used to identify conditions of interest that conform to desired objective functions, such as stationary conditions, ramps or waves of wind speed, and changes in wind direction.Total variation is somewhat sensitive to the presence of outliers in the input data, and the method is best complemented by quality-control procedures to ensure reliable results.The direct detection and classification of events or conditions of interest within atmospheric data sets is vital to developing our understanding of wind plant response and to the formulation of forecasting and control models.
Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.
Here, the exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.
In fusion neutral gas analysis, such as with the Diagnostic Residual Gas Analyzer (DRGA) for ITER, the primary measurement range of interest comprises the low-amu species (1 to 6), especially deuterium and helium. The challenge in successfully obtaining accurate measurements is two-fold. First, the sensitivity of the method must be sufficient to resolve trace amounts accurately; typically, one percent or less. Second, the gas signal from the fusion processes must be free of bias caused by the latent presence (from system outgassing and/or vacuum backstreaming) of these gases to enable accurate interpretation of the measured signal. This latter criterion can be problematic for the lightest gases since there is a propensity for some fraction of the pumped gas load to undergo a phenomenon known as backstreaming. This behavior is manifested in pumping systems for gas properties related to relative atomic weight (lightest) and size (smallest). Backstreaming results in a significant amount of the pumped gas undertaking a reverse flow and re-entering the measurement region; thus, contaminating the forward, real-time measurement. To fully eliminate this adverse effect, a conductance-limiting device – or orifice – has been installed in the high-vacuum pumping system of the present ITER DRGA prototype. The system was already equipped with a secondary turbomolecular pump (TMP), but with limited effectiveness against backstreaming in the inter-pump volume (IPV). This orifice is placed within the suction inlet coupling of the secondary TMP, which is downstream of the IPV. Its objective is to eliminate the backstreaming phenomenon by increasing the back pressure in the IPV. However, the orifice sizing must take into consideration other factors, such as the diagnostic measurement objectives. For example, in the ITER DRGA, one of the measurement requirements is a dynamic response time of ~1s. Fortunately, an added benefit of the pumping restriction created by the orifice is that the upstream pressure increase is beneficial for the DRGA’s optical gas analysis (OGA) sensors. These sensors are attached to the IPV in the present design. The glow discharges, when used as an OGA light source, will typically have a brighter light emission with increasing plasma cell pressure. In addition to the fusion machine research sector, there are other potential applications of this pumping technique where the monitoring of lighter gas concentrations is essential, such as the photolithography process for the semiconductor fabrication of integrated circuits. This presentation will describe the vacuum system used to demonstrate a process to eliminate backstreaming as well as show test results to verify the accomplishment of this critical objective.
We numerically investigate an adaptive version of the parareal algorithm in the context of molecular dynamics. This adaptive variant has been originally introduced in [1]. We focus here on test cases of physical interest where the dynamics of the system is modelled by the Langevin equation and is simulated using the molecular dynamics software LAMMPS. In this work, the parareal algorithm uses a family of machine-learning spectral neighbor analysis potentials (SNAP) as fine, reference, potentials and embedded-atom method potentials (EAM) as coarse potentials. We consider a self-interstitial atom in a tungsten lattice and compute the average residence time of the system in metastable states. Our numerical results demonstrate significant computational gains using the adaptive parareal algorithm in comparison to a sequential integration of the Langevin dynamics. We also identify a large regime of numerical parameters for which statistical accuracy is reached without being a consequence of trajectorial accuracy.
Mathematical models have shed light on the dynamics of HIV- 1 infection in vivo. In this paper, we generalize continuous mathematical models of drug therapy for HIV-1 by Perelson et al. (Science 271:1582–1586, 1996) and Perelson and Nelson (SIAM Rev 41:3–44, 1999) on time scales, i.e., a nonempty closed subset of real numbers in order to derive new discrete models that predict the total concentration of plasma virus as a function of time. Here, one of our main goals is to compare discrete mathematical models with the continuous model in Perelson et al. (1996) where HIV infected patients were given protease inhibitors and sampled frequently thereafter. For the comparison, we use experimental data collected in Perelson et al. (1996) and estimate the parameters such as the virion clearance rate and the rate of loss of infected cells by fitting the total concentration of plasma virus to this data set. Our results show that discrete systems describe the best fit. In the previous models of this study, the efficacy of protease inhibitor is assumed to be perfect. Motivated by Perelson and Nelson (1999), we end the paper with a mathematical model of imperfect protease inhibitor and reverse transcriptase (RT) inhibitor combination therapy of HIV-1 infection on time scales with its stability analysis.
Exciton dynamics of perovskite nanoclusters has been investigated for the first time using femtosecond transient absorption (TA) and time-resolved photoluminescence (TRPL) spectroscopy. The TA results show two photoinduced absorption signals at 420 and 461 nm and a photoinduced bleach (PB) signal at 448 nm. The analysis of the PB recovery kinetic decay and kinetic model uncovered multiple processes contributing to electron-hole recombination. The fast component (~8 ps) is attributed to vibrational relaxation within the initial excited state, and the medium component (~60 ps) is attributed to shallow carrier trapping. The slow component is attributed to deep carrier trapping from the initial conduction band edge (~666 ps) and the shallow trap state (~40 ps). The TRPL reveals longer time dynamics, with modeled lifetimes of 6.6 and 93 ns attributed to recombination through the deep trap state and direct band edge recombination, respectively. In conclusion, the significant role of exciton trapping processes in the dynamics indicates that these highly confined nanoclusters have defect-rich surfaces.