Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dependability”

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 1,459 records · Page 81

Thickness- and Twist-Angle-Dependent Interlayer Excitons in Metal Monochalcogenide Heterostructures

Interlayer excitons, or bound electron-hole pairs whose constituent quasiparticles are located in distinct stacked semiconducting layers, are being intensively studied in heterobilayers of two dimensional semiconductors. They owe their existence to an intrinsic type-II band alignment between both layers that convert these into p-n junctions. Here, we unveil a pronounced interlayer exciton (IX) in heterobilayers of metal monochalcogenides, namely gamma-InSe on epsilon-GaSe, whose pronounced emission is adjustable just by varying their thicknesses given their number of layers dependent direct bandgaps. Time-dependent photoluminescense spectroscopy unveils considerably longer interlayer exciton lifetimes with respect to intralayer ones, thus confirming their nature. The linear Stark effect yields a bound electron-hole pair whose separation d is just (3.6 ± 0.1) Å with d being very close to d Se = 3.4 Å which is the calculated interfacial Se separation. The envelope of IX is twist angle dependent and describable by superimposed emissions that are nearly equally spaced in energy, as if quantized due to localization induced by the small moiré periodicity. Furthermore, these heterostacks are characterized by extremely flat interfacial valence bands making them prime candidates for the observation of magnetism or other correlated electronic phases upon carrier doping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Size-Dependent Biexciton Quantum Yields and Carrier Dynamics of Quasi-Two-Dimensional Core/Shell Nanoplatelets

Quasi-two-dimensional nanoplatelets (NPLs) possess fundamentally different excitonic properties from zero-dimensional quantum dots. We study lateral size-dependent photon emission statistics and carrier dynamics of individual NPLs using second-order photon correlation (g( 2 )(τ)) spectroscopy and photoluminescence (PL) intensity-dependent lifetime analysis. Room-temperature radiative lifetimes of NPLs can be derived from maximum PL intensity periods in PL time traces. It first decreases with NPL lateral size and then stays constant, deviating from the electric dipole approximation. Analysis of the PL time traces further reveals that the single exciton quantum yield in NPLs decreases with NPL lateral size and increases with protecting shell thickness, indicating the importance of surface passivation on NPL emission quality. Second-order photon correlation (g( 2 )(τ)) studies of single NPLs show that the biexciton quantum yield is strongly dependent on the lateral size and single exciton quantum yield of the NPLs. In large NPLs with unity single exciton quantum yield, the corresponding biexciton quantum yield can reach unity. In conclusion, these findings reveal that by careful growth control and core–shell material engineering, NPLs can be of great potential for light amplification and integrated quantum photonic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvent- and Wavelength-Dependent Photoluminescence Relaxation Dynamics of Carbon Nanotube sp 3 Defect States

Photoluminescent sp 3 defect states introduced to single wall carbon nanotubes (SWCNTs) through low-level covalent functionalization create new photophysical behaviors and functionality as a result of defect sites acting as exciton traps. Evaluation of relaxation dynamics in varying dielectric environments can aid in advancing a more complete description of defect-state relaxation pathways and electronic structure. For this study, we exploit helical wrapping polymers as a route to suspending (6,5) SWCNTs covalently functionalized with 4-methoxybenzene in solvent systems including H 2 O, D 2 O, methanol, dimethylformamide, tetrahydrofuran, and toluene, spanning a range of dielectric constants from 80 to 3. Defect-state photoluminescence decays were measured as a function of emission wavelength and solvent environment. Emission decays are biexponential, with short lifetime components on the order of 65 ps and long components ranging from around 100 to 350 ps. Both short and long decay components increase as emission wavelength increases, while only the long lifetime component shows a solvent dependence. We demonstrate that the wavelength dependence is a consequence of thermal detrapping of defect-state excitons to produce mobile E 11 excitons, providing an important mechanism for loss of defect-state population. Deeper trap states (i.e., those emitting at longer wavelengths) result in a decreased rate for thermal loss. The solvent-independent behavior of the short lifetime component is consistent with its assignment as the characteristic time for redistribution of exciton population between bright and dark defect states. Furthermore, the solvent dependence of the long lifetime component is shown to be consistent with relaxation via an electronic to vibrational energy transfer mechanism, in which energy is resonantly lost to solvent vibrations in a complementary mechanism to multiphonon decay processes.

36 MATERIALS SCIENCE↗

Quantifying Wavelength-Dependent Plasmonic Hot Carrier Energy Distributions at Metal/Semiconductor Interfaces

Hot carriers generated from the nonradiative decay of localized surface plasmons are capable of driving charge-transfer reactions at the surfaces of metal nanostructures. Photocatalytic devices utilizing plasmonic hot carriers are often based on metal nanoparticle/semiconductor heterostructures owing to their efficient electron–hole separation ability. The rapid thermalization of hot carriers generated at the metal nanoparticles yields a distribution of carrier energies that determines the capability of the photocatalytic device to drive redox reactions. Here in this paper, we quantify the thermalized hot carrier energy distribution generated at Au/TiO 2 nanostructures using wavelength-dependent scanning electrochemical microscopy and a series of molecular probes with different redox potentials. We determine the quantum efficiencies and oxidizing power of the hot carriers from wavelength-dependent reaction rates and photocurrent across the metal/semiconductor interface. The wavelength-dependent reaction efficiency tracks the surface plasmon resonance spectrum of the Au nanoparticles, showing that the reaction is facilitated by plasmon excitation, while the responses from molecules with different redox potentials shed light on the energy distribution of the hot holes generated at metal nanoparticle/semiconductor heterostructures. The results provide important insight into the energies of the plasmon-generated hot carriers and quantum efficiencies of plasmonic photocatalytic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure–Transport Properties Governing the Interplay in Humidity-Dependent Mixed Ionic and Electronic Conduction of Conjugated Polyelectrolytes

Polymeric mixed ionic-electronic conductors (MIECs) are of broad interest in the field of energy storage and conversion, opto-electronics, and bioelectronics. A class of polymeric MIECs are conjugated polyelectrolytes (CPEs), which possess a p-conjugated backbone imparting electronic transport characteristics along with side chains comprised of a pendant ionic group to allow for ionic transport. Here, our study focuses on the humidity-dependent structure-transport properties of poly[3-(potassium-n-alkanoate) thiophene-2,5-diyl] (P3KnT) CPEs with varied side-chain lengths of n = 4, 5, 6, and 7. UV-Vis spectroscopy along with electronic paramagnetic resonance (EPR) spectroscopy reveal the infiltration of water leads to a hydrated, self-doped state that allows for electronic transport. The resulting humidity-dependent ionic conductivity (σ i ) of the thin films shows a monotonic increase with relative humidity (RH) while electronic conductivity (σ e ) follows a nonmonotonic profile. The values of σ e continue to rise with increasing RH reaching a local maximum after which σ e begins to decrease. P3KnTs with higher n values demonstrate greater resiliency to increasing RH before suffering decrease in σ e . This drop in σ e is attributed to two factors. First, disruption of the locally-ordered π-stacked domains observed through in situ humidity-dependent grazing incidence wide angle X-ray scattering (GIWAXS) experiments can account for some of the decrease in σ e . A second and more dominant factor is attributed to the swelling of the amorphous domains where electronic transport pathways connecting ordered domains are impeded. P3K7T is most resilient to swelling (based on ellipsometry and water uptake measurements) where sufficient hydration allows for high σ i (1.0 × 10 -1 S/cm at 95% RH) while not substantially disrupting σ e (1.7 × 10 -2 S/cm at 85% RH and 8.0 × 10 -3 S/cm at 95% RH). Overall, our study highlights the complexity of balancing electronic and ionic transport in hydrated CPEs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multivalent Macromolecules Redirect Nucleation-Dependent Fibrillar Assembly into Discrete Nanostructures

Manipulating the size and shape of noncovalent multivalent assemblies is an ongoing challenge in the field of supramolecular polymers. Following a mechanistic approach, we reasoned that nucleation–elongation kinetics presents unique opportunities for controlled growth since the final outcome is likely to depend on the structure and dynamics of critical-nucleus formation. Taking fibrillar assembly of amyloid β (Aβ) peptide as the model system of nucleation-dependent supramolecular polymerization, here we report multivalent polymer–peptide conjugates (mPPCs) that redirect fibrillar assembly of Aβ to form discrete nanostructures. The mPPCs were rationally designed to target Aβ intermediates formed prior to critical nucleation. Atomic force microscopy and transmission electron microscopy studies show that in the presence of mPPCs, Aβ self-assembles into zero-dimensional discrete nanostructures with lateral dimensions approximately in 5–35 nm, while Aβ alone self-assembles into one-dimensional fibrils in micrometer. Thioflavin T kinetics fluorescence assays demonstrate that mPPCs suppress Aβ fibrillogenesis. As a result, the mPPCs may thus represent a prototypical molecular design of multivalent macromolecules able to control the final shape of supramolecular polymers assembled via a nucleation-dependent mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mo Atom Rearrangement Drives Layer-Dependent Reactivity in Two-Dimensional MoS 2

Two-dimensional (2D) materials offer a valuable platform for manipulating and studying chemical reactions at the atomic level, owing to the ease of controlling their microscopic structure at the nanometer scale. While extensive research has been conducted on the structure-dependent chemical activity of 2D materials, the influence of structural transformation during the reaction has remained largely unexplored. In this work, we report the layer-dependent chemical reactivity of MoS 2 during a nitridation atomic substitution reaction and attribute it to the rearrangement of Mo atoms. Our results show that the chemical reactivity of MoS 2 decreases as the number of layers is reduced in the few-layer regime. In particular, monolayer MoS 2 exhibits significantly lower reactivity compared with its few-layer and multilayer counterparts. Atomic-resolution transmission electron microscopy (TEM) reveals that MoN nanonetworks form as reaction products from monolayer and bilayer MoS 2 , with the continuity of the MoN crystals increasing with layer number, consistent with the local conductivity mapping data. The layer-dependent reactivity is attributed to the relative stability of the hypothetically formed MoN phase, which retains the number of Mo atomic layers present in the precursor. Specifically, the low chemical reactivity of monolayer MoS 2 is attributed to the high energy cost associated with Mo atom diffusion and migration necessary to form multilayer Mo lattices in the thermodynamically stable MoN phase. In conclusion, this study underscores the critical role of lattice rearrangement in governing chemical reactivity and highlights the potential of 2D materials as versatile platforms for advancing the understanding of materials chemistry at the atomic scale.

Chemical reactivity↗

Driving Force Dependent Photoinduced Charge Transfer Dynamics in Polymer-Wrapped Semiconducting Single-Walled Carbon Nanotubes

Here, we investigate the thermodynamic driving-force dependences of photoinduced charge separation (CS) and subsequent charge transfer dynamics in single-walled carbon nanotube (SWNT)–perylenediimide (PDI) donor–acceptor (D–A) superstructures. Pump–probe spectroscopy reveals that [SWNT (•+)n ]-(PDI –• ) n CS states form on an ∼100 fs time scale following photoexcitation; these dynamics are invariant across an ∼400 mV driving force range, indicating that SWNT hole polaron formation time scales are determined by nanotube lattice and solvent relaxation. These CS states feature SWNT hole polarons adjacent to (geminate) and nearby (nongeminate) PDI radical anions. Analysis of the free energy dependence for charge recombination (CR) of [SWNT •+ ] geminate -(PDI –• ) CS states highlights an ∼2 meV value for D–A electronic coupling (HAB) and ∼0.93 eV for the total reorganization energy (λ T ). A corresponding driving force dependence of the CR dynamics for [SWNT •+ ] nongeminate -(PDI –• ) CS states indicates a diminished H AB value (∼0.6 meV) and a larger λ T (∼1.1 eV), consistent with larger transfer distances. SWNT excitons that persist following photoinduced CS drive photooxidation of PDI –• components of [SWNT (•+)n ]-(PDI –• ) n CS states ( 1 SWNT* + PDI •– → PDI + SWNT •– ); this reaction manifests a significantly reduced λ T value (∼0.67 eV) as the initially prepared SWNT reduced state bears the character of a conduction band injected electron ([SWNT •– ] CB ). This reaction thus gives rise to relaxed, nongeminate SWNT electron and hole polarons on the same nanotube; these polarons react on a 10 2 ps time scale independent of the electronic structure of these SWNT-PDI superstructures.

anions↗

Explanation of Dramatic pH-Dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High pH

Hydrogen oxidation reaction (HOR) and hydrogen evolution reaction (HER) are both 2 orders slower in alkaline electrolyte than in acidic electrolyte, but no explanation has been provided. The first step toward understanding this dramatic pH-dependent HOR/HER performance is to explain the pH-dependent hydrogen binding to the electrode, a perplexing behavior observed experimentally. In this work, we carried out Quantum Mechanics Molecular Dynamics (QMMD) with explicit considerations of solvent and applied voltage ( U ) to in situ simulate water/Pt(100) interface in the condition of under-potential adsorption of hydrogen ( H UPD ). We found that as U is made more negative, the electrode tends to repel water, which in turn increases the hydrogen binding. We predicted a 0.13 eV increase in hydrogen binding from pH = 0.2 to pH = 12.8 with a slope of 10 meV/pH, which is close to the experimental observation of 8 to 12 meV/pH. Thus, we conclude that the changes in water adsorption are the major causes of pH-dependent hydrogen binding on a noble metal. Furthermore, the new insight of critical role of surface water in modifying electrochemical reactions provides a guideline in designing HER/HOR catalyst targeting for the alkaline electrolyte.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Observation of Symmetry-Dependent Electron–Phonon Coupling in Black Phosphorus

Electron–phonon coupling in two-dimensional nanomaterials plays a fundamental role in determining their physical properties. Such interplay is particularly intriguing in semiconducting black phosphorus (BP) due to the highly anisotropic nature of its electronic structure and phonon dispersions. Here we report the direct observation of symmetry-dependent electron–phonon coupling in BP by performing the polarization-selective resonance Raman measurement in the visible and ultraviolet regimes, focusing on the out-of-plane A g 1 and in-plane A g 2 phonon modes. Their intrinsic resonance Raman excitation profiles (REPs) were extracted and quantitatively compared. The in-plane A g 2 mode exhibits remarkably strong resonance enhancement across the excitation wavelengths when the excitation polarization is parallel to the armchair (A g 2 //AC) direction. In contrast, a dramatically weak resonance effect was observed for the same mode with the polarization parallel to zigzag (A g 2 //ZZ) direction and for the out-of-plane A g 1 mode (A g 1 //AC and A g 1 //ZZ). Analysis on quantum perturbation theory and first-principles calculations on the anisotropic electron distributions in BP demonstrated that electron–phonon coupling considering the symmetry of the involved excited states and phonon vibration patterns is responsible for this phenomenon. Further analysis of the polarization-dependent REPs for A g phonons allows us to resolve the existing controversies on the physical origin of Raman anomaly in BP and its dependence on excitation energy, sample thickness, phonon modes, and crystalline orientation. Finally, our study gives deep insights into the underlying interplay between electrons and phonons in BP and paves the way for manipulating the electron–phonon coupling in anisotropic nanomaterials for future device applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Size- and Composition-Dependent Radio Frequency Magnetic Permeability of Iron Oxide Nanocrystals

We investigate the size- and composition-dependent ac magnetic permeability of superparamagnetic iron oxide nanocrystals for radio frequency (RF) applications. The nanocrystals are obtained through high-temperature decomposition synthesis, and their stoichiometry is determined by Mossbauer spectroscopy. Two sets of oxides are studied: (a) as-synthesized magnetite-rich and (b) aged maghemite nanocrystals. All nanocrystalline samples are confirmed to be in the superparamagnetic state at room temperature by SQUID magnetometry. Through the one-turn inductor method, the ac magnetic properties of the nanocrystalline oxides are characterized. In magnetite-rich iron oxide nanocrystals, size-dependent magnetic permeability is not observed, while maghemite iron oxide nanocrystals show clear size dependence. The inductance, resistance, and quality factor of hand-wound inductors with a superparamagnetic composite core are measured. The superparamagnetic nanocrystals are successfully embedded into hand-wound inductors to function as inductor cores.

Yun, H↗

Simultaneous Dependence of the Earthquake-Size Distribution on Faulting Style and Depth

We analyze two high-quality Southern Californian earthquake catalogues, one with focal mechanisms, to statistically model and test for dependencies of the earthquake-size distribution, the b values, on both faulting style and depth. In our null hypothesis, b is assumed constant. We then develop and calibrate one model based only on faulting style, another based only on depth dependence and two models that assume a simultaneous dependence on both parameters. We develop a new maximum-likelihood estimator corrected for the degrees of freedom to assess models' performances. Our results show that all models significantly reject the null hypothesis. The best performing is the one that simultaneously takes account of depth and faulting style. Our results suggest that differential stress variations in the Earth's crust systematically influence b values and that this variability should be considered for contemporary seismic hazard studies.

54 ENVIRONMENTAL SCIENCES↗

Scale-Dependent Bedform Migration and Deformation in the Physical and Spectral Domains

Multi-scale bathymetries observed in laboratory channels and natural rivers have posed several challenges to the description of bedform geometry and kinematics. To better quantify scale-dependent migration velocity, high resolution spatio-temporal bed evolution data are analyzed using three independent methods: (i) a bedform tracking method that identifies individual bedforms in longitudinal bed elevation profiles, (ii) a cross-correlation based, image matching technique in consecutive bathymetry scans similar to Particle Image Velocimetry, and (iii) two-dimensional frequency-wavenumber spectra marking a dispersion relation between bedform period and length, allowing estimation of Fourier-based, scale-dependent migration velocities. The comparative studies show that the spectral and bedform tracking methods present small differences in the migration velocity of medium-large bedforms. However, significant deviations occur in the small bedforms, mostly due to bedform sheltering or amplified exposure to the flow, which also contribute to their enhanced scale-dependent deformation. The analysis also shows that some form of scale separation is needed to extract the large and slowly migrating bedforms as they are masked by secondary features when several bedform orders coexist. Bathymetric Image Velocimetry reasonably estimates the averaged migration velocity and exhibits potential to study relations between spanwise bed surface velocity and local bedform kinematics.

58 GEOSCIENCES↗

Explaining Forcing Efficacy With Pattern Effect and State Dependence

The magnitude of global surface temperature change in response to unit radiative forcing depends on the type and magnitude of forcing agent—a concept known as a “forcing efficacy.” However, the mechanisms behind the forcing efficacy are still unclear. In this study, we perform a set of simulations using CESM1 to calculate the efficacy of 10 different forcing agents defined in terms of fixed-SST effective radiative forcing, and then use a Green's function approach to show that each forcing efficacy can be largely understood in terms of the radiative feedbacks associated with the different surface temperature patterns induced by the forcing agents (a pattern effect). We also quantify how the state dependence of feedbacks on global mean surface temperature anomalies impacts forcing efficacies. The results show that the forcing efficacy can be well reconstructed with a combination of pattern effect and state dependence.

58 GEOSCIENCES↗

An Open‐Source, Physics‐Based, Tropical Cyclone Downscaling Model With Intensity‐Dependent Steering

Abstract An open‐source, physics‐based tropical cyclone (TC) downscaling model is developed, in order to generate a large climatology of TCs. The model is composed of three primary components: (a) a random seeding process that determines genesis, (b) an intensity‐dependent beta‐advection model that determines the track, and (c) a non‐linear differential equation set that determines the intensification rate. The model is entirely forced by the large‐scale environment. Downscaling ERA5 reanalysis data shows that the model is generally able to reproduce observed TC climatology, such as the global seasonal cycle, genesis locations, track density, and lifetime maximum intensity distributions. Inter‐annual variability in TC count and power‐dissipation is also well captured, on both basin‐wide and global scales. Regional TC hazard estimated by this model is also analyzed using return period maps and curves. In particular, the model is able to reasonably capture the observed return period curves of landfall intensity in various sub‐basins around the globe. The incorporation of an intensity‐dependent steering flow is shown to lead to regionally dependent changes in power dissipation and return periods. Advantages and disadvantages of this model, compared to other downscaling models, are also discussed.

Meteorology & Atmospheric Sciences↗

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗

Quantifying the Effect of Pore‐Size Dependent Wettability on Relative Permeability Using Capillary Bundle Model

Abstract Relative permeability is a key parameter for characterizing the multiphase flow dynamics in porous media at macroscopic scale while it can be significantly impacted by wettability. Recently, it has been reported in microfluidic experiments that wettability is dependent on the pore size (Van Rooijen et al., 2022). To investigate the effect of pore‐size‐dependent wettability on relative permeability, we propose a theoretical framework informed by digital core samples to quantify the deviation of relative permeability curves due to wettability change. We find that the significance of impact is highly dependent on two factors: (i) the function between contact angle and pore size (ii) overall pore size distribution. Under linear function, this impact can be significant for tight porous media with a maximum deviation of 1,000%.

Yu, Siqin↗

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗