Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ionic diffusivity”

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 343 records · Page 19

Understanding interfacial chemistry of positive bias high-voltage degradation in photovoltaic modules

Photovoltaic module degradation from a high system voltage is a prevalent degradation mode in the field, where the enabling degradation mechanisms are inherently dependent on the voltage bias polarity of the installed system. In this study, the effects of positive bias on module performance are confirmed and the underlying chemical degradation processes are more thoroughly investigated to reveal different degradation pathways from those previously reported in negative bias studies. When cells are under +1000 V stress, crystalline silicon mini-modules with poly(ethylene-co-vinyl acetate) (EVA) encapsulant demonstrated a significant photocurrent loss due to EVA discoloration and delamination from increased chemical reactivity at the front-side EVA/cell metallization interface. Brown discoloration of the EVA encapsulant near the cell gridlines is linked to an electrochemical reaction at the Ag gridlines under hot and humid conditions (85 °C, 85% relative humidity). Chemical compositional analysis using X-ray photoelectron spectroscopy (XPS) confirmed that the discoloration is attributed to the formation of silver sulfide (Ag 2 S) and/or silver oxide (Ag 2 O) species at the EVA/Ag gridline interface. The subsequent migration of Ag ions from the cell gridlines into the bulk of the EVA was evident from XPS depth profiling and optical microscopy. However, the Ag signal was not detected at the EVA/glass interface, inferring limited ionic transport through the nominally 0.45 mm thick encapsulant. For the samples studied herein, the sulfur is believed by the process of elimination to come from the ambient air, diffusing into the module through the permeable polymer backsheet.

14 SOLAR ENERGY↗

Framework Short-Range Order Observed in a Spinel-Type Li Superionic Conductor

Solid-state superionic conductors are characterized by rich structural disorders. Though structural complexities are central to their functionalities, they often give rise to short-range order that eludes detection by conventional diffraction-based techniques and is thus overlooked in establishing precise structure-property relationships. In this work, we synthesized single crystals of a recently discovered lithium (Li) superionic conductor Li16.2(1)In9.00(2)Sn1.10(1)O23.8 (LISO) for in-depth characterizations of structural subtleties. LISO exhibits an unusual spinel-like phase with significant Li overstoichiometry and a face-sharing Li network. Single-crystal neutron diffraction confirms significant Li disorder, as manifested in Li site splitting and partial occupancy. More importantly, synchrotron diffuse scattering combined with 3D-ΔPDF analysis and Monte Carlo simulations reveal short-range order in the nonalkali framework that might contribute to the phase stability and ionic conductivity. This work showcases an example in which subtle local energetics can be directly visualized in structurally disordered ionic conductors.

Chen, Yu↗

Salt Transport in Crosslinked Hydrogel Membranes Containing Zwitterionic Sulfobetaine Methacrylate and Hydrophobic Phenyl Acrylate

Produced water is a by-product of industrial operations, such as hydraulic fracturing for increased oil recovery, that causes environmental issues since it includes different metal ions (e.g., Li+, K+, Ni2+, Mg2+, etc.) that need to be extracted or collected before disposal. To remove these substances using either selective transport behavior or absorption-swing processes employing membrane-bound ligands, membrane separation procedures are promising unit operations. This study investigates the transport of a series of salts in crosslinked polymer membranes synthesized using a hydrophobic monomer (phenyl acrylate, PA), a zwitterionic hydrophilic monomer (sulfobetaine methacrylate, SBMA), and a crosslinker (methylenebisacrylamide, MBAA). Membranes are characterized according to their thermomechanical properties, where an increased SBMA content leads to decreased water uptake due to structural differences within the films and to more ionic interactions between the ammonium and sulfonate moieties, resulting in a decreased water volume fraction, and Young’s modulus increases with increasing MBAA or PA content. Permeabilities, solubilities, and diffusivities of membranes to LiCl, NaCl, KCl, CaCl2, MgCl2, and NiCl2 are determined by diffusion cell experiments, sorption-desorption experiments, and the solution-diffusion relationship, respectively. Permeability to these metal ions generally decreases with an increasing SBMA content or MBAA content due to the corresponding decreasing water volume fraction, and the permeabilities are in the order of K+ > Na+ > Li+ > Ni2+ > Ca2+ > Mg2+ presumably due to the differences in the hydration diameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adapting confocal Raman microscopy for in situ studies of redox transformations at electrode-electrolyte interfaces

Confocal Raman microscopy was applied to quantify redox species present within the diffusion layer adjacent to an electrode surface under potentiostatic control. A glass microscope coverslip with a thin indium tin oxide (ITO) coating served as both the working electrode and optical window for a microscope-stage mountable spectroelectrochemical cell. A high numerical aperture objective mounted in an inverted microscope frame just below the stage brought excitation radiation through the coverslip window and to a tight focus a few micrometers above the ITO film surface. Species diffusing into the confocal probe volume defined by the excitation beam focus and the collected light region were detected, identified and quantified based on their Raman scattering frequencies and intensities. Here, in measurements that interrogated the interconversion of ferrocyanide and ferricyanide ions as a function of applied voltage, least-squares regression analysis of spectral datasets predicted the formal potential and relative surface concentrations of the ions in good agreement with the expected Nernstian response. Preliminary studies of methyl viologen reduction at an ITO film/Nafion membrane interface were conducted and show the possibility for estimation of mass transport coefficients of redox species within ionic polymer materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗

Nonlinear ion transport mediated by induced charge in ultrathin nanoporous membranes

Ultrathin membranes with nanoporous conduits show promise for ionic separations and desalination applications, but the mechanisms underlying the nonlinear ionic transport observed in these systems are not well understood. Here, in this work, we demonstrate how induced charge at membrane interfaces can lead to nonlinear ionic transport and voltage-dependent conductance through such channels. The application of an electric field on a polarizable membrane leads to induced charges at the membrane interfaces. The induced charges in turn are screened by diffuse charges in the electrolyte, which are acted upon by the electric field. For extremely thin membranes, the induced charge effect can be significant even for moderate applied voltages commonly used in experiments. We apply a continuum Poisson-Nernst-Planck model to characterize the current-voltage behavior of ultrathin membranes over a wide parameter space. The predictions of the model are compared to recent experiments on graphene and MoS 2 membranes in an electric field. We expect the role of induced charge to be especially pronounced in the limit of atomically thin membranes.

36 MATERIALS SCIENCE↗

Interfacial solvation and slow transport of hydrated excess protons in non-ionic reverse micelles

This work employs molecular dynamics simulations to investigate the solvation and transport properties of hydrated excess protons (with a hydronium-like core structure) in non-ionic Igepal CO-520 reverse micelles of various sizes in a non-polar solvent. Multiscale Reactive Molecular Dynamics (MS-RMD) simulations were used to describe vehicular and hopping diffusion during the proton transport process. As detailed herein, an excess proton shows a marked tendency to localize in the interfacial region of micellar water pools. Slow proton transport was observed which becomes faster with increasing micellar size. Further analysis reveals that the slow diffusion of an excess proton is a combined result of slow water diffusion and the low proton hopping rate. This study also confirms that a low proton hopping rate in reverse micelles stems from the interfacial solvation of hydrated excess protons and the immobilization of interfacial water. The low water density in the interfacial region makes it difficult to form a complete hydrogen bond network near the hydrated excess proton, and therefore locks in the orientation of hydrated proton cations. The immobilization of the interfacial water also slows the relaxation of the overall hydrogen bond network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Membrane rejection of nitrogen compounds

Rejection characteristics of nitrogen compounds were examined for reverse osmosis, nanofiltration, and low-pressure reverse osmosis membranes. The rejection of nitrogen compounds is explained by integrating experimental results with calculations using the extended Nernst-Planck model coupled with a steric hindrance model. The molecular weight and chemical structure of nitrogen compounds appear to be less important in determining rejection than electrostatic properties. The rejection is greatest when the Donnan potential exceeds 0.05 V or when the ratio of the solute radius to the pore radius is greater than 0.8. The transport of solute in the pore is dominated by diffusion, although convective transport is significant for organic nitrogen compounds. Electromigration contributes negligibly to the overall solute transport in the membrane. Urea, a small organic compound, has lower rejection than ionic compounds such as ammonium, nitrate, and nitrite, indicating the critical role of electrostatic interaction in rejection. This suggests that better treatment efficiency for organic nitrogen compounds can be obtained after ammonification of urea.

NASA Discipline Life Support Systems↗

Methods—Meso-Scale Electrodes for Characterizing Diffusion-Reaction Properties of Redox-Active Organics in Viscous Electrolytes

Deep Eutectic Solvents (DESs) have recently gained interest as flow battery electrolytes. Their advantages include a wider electrochemical stability window compared to aqueous electrolytes, higher solubility for redox-active species, and negligible volatility. However, DESs are often highly viscous, and suffer from low ionic conductivities. This can make assessing redox kinetics difficult when attempting to determine their viability for energy storage. In classical voltammetric measurements, low ionic conductivity manifests as high solution resistance, thereby requiring "live" compensation of the electrolyte ohmic drop when performing fast-scan voltammetry. An uncompensated or inadequately-compensated ohmic drop leads to misinterpretation of the voltammetric behavior, e.g., assessing reversibility vs. irreversibility of a redox reaction. Here, we present micro-fabricated electrodes as facile "meso-scale" electrodes, which overcome these issues by nearly eliminating the ohmic drop while retaining uniformity of the current distribution over the electrode surface. Their use in precise transport-kinetics measurements is demonstrated using a redox-active organic, i.e., 4-Hydroxy-TEMPO in an aqueous medium and in ethaline, which is a viscous DES. This study provides a methodical approach to design and to implement voltammetry experiments using meso-scale electrodes leading to reliable measurements of diffusion-reaction properties of 4-Hydroxy-TEMPO.

25 ENERGY STORAGE↗

On-Demand Designing of Cathode Internal Surface Architecture for Dramatic Enhancement of SOFC Performance and Durability

This project is aimed to design and modify the internal surfaces of porous composite cathode from currently commercially viable Solid Oxide Fuel Cells (SOFCs), using additive manufacturing process of Atomic Layer Deposition (ALD). The material systems being investigated are commercial composite electrodes complex three-dimensional topographies. In term of the chemistry of the ALD layer applied on the internal surface of the porous cathode, this project has employed commercially relevant electrolyte, electrocatalyst and noble metal materials set. Such materials are fully compatible with the commercial fuel cells, and this project has developed special nanostructure on the surface of the commercial composite cathodes. The formation of the designed nanoarchitecture on the surface of SOFC cathode has been achieved through precise control of ALD parameters and their effect on overall cell performance and resultant electrochemical reaction mechanism of cathodes has been investigated through full cell electrochemical performance testing and nanostructure characterization by transmission electron microscopy (TEM). Under the support of this award, following has been achieved: (1). For cathode materials in solid oxide fuel cells (SOFCs), such as perovskite mixed conductor La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-x (LSCF), cation surface segregation and consequently losing conductivity and active sites for the oxygen reduction reaction (ORR) are problematic. To mitigate the cation segregation and enhance SOFC durability, further decorating the internal backbone surface using the desired electrocatalysts could be a promoting approach. Commonly, the cation segregation such as Sr is very volatile, so the effective surface decoration is ideally conformal. Nevertheless, the conformal surface coating would inevitably alter the ORR pathways that initially take place on the surface of the backbone. To reveal the impact of the conformal coating on both the catalytic activity and the conductivity of the cathode, the unary electrocatalyst of Pt or CoO x , was applied to the LSCF/SDC composite electrode of inherently functional SOFCs, respectively. Both ALD coating layers evolve strong interaction with the LSCF composite cathode. Upon operations, the Pt coating layer remains conformal on LSCF grain surfaces but turns into discrete particles on SDC grain surfaces. Meanwhile, CoO x conformal coating grows to be the discrete nanograins on both the LSCF and SDC grains. ALD coating of the cathode alone reduces the ohmic resistance up to 28 % for the entire cells. The increased conductivity induced by the ALD coating of Pt or CoO x is ascribed to different mechanisms. For the inherent functional SOFCs, the present study presents a novel and feasible approach to apply a conformal, dense coating layer on the surface of a mixed conductor, simultaneously increasing the conductivity and durability of the SOFC cathode. (2). High resistance of the oxygen electrode still significantly hinders the state-of-the-art Solid Oxide Fuel Cells (SOFCs). In particular, for an oxygen electrode consisting of mixed electronic and ionic conductors, such as perovskite lanthanum strontium cobalt ferrite (LSCF), it deteriorates due to its low chemical stability of the grain surface. Such degradation is often associated with the segregation of cations. To prevent the cation surface segregation and its resultant perovskite phase decomposition, we demonstrate a conformal ultra-thin (7-10 nm) surface heterogeneous coating layer consisting of subjacent discrete Pt nanoparticles capped with a superjacent fully dense conformal CoO x layer. The performance studies indicate the ALD coating reduces the cell series resistance by up to 40 %. The conformal CoO x layer consists of randomly orientated but single-layered nanograins, with high-density intergranular and surface grain boundaries serving as the electrochemical reaction sites and facilitating mass transport. The conformal coating layer appears to have successfully suppressed the Sr outward diffusion and confined the Sr enriched layer to a ~ 2 nm interface perovskite phase between the coating layer and the LSCF grain surface. Moreover, this ultra-thin Sr enriched perovskite layer presumably possesses high oxygen vacancy and high ionic conductivity and further imposes tensile strain to the LSCF grain surfaces. With the combination of a conformal CoO x nanoionics, Sr enriched layer, and its strained interface, the ALD coating induced surface layer is estimated to have a conductivity of ~ 1.27x10 4 S/cm, which is over two orders magnitude of that from LSCF at 750 ºC.

30 DIRECT ENERGY CONVERSION↗

High-speed mapping of surface charge dynamics using sparse scanning Kelvin probe force microscopy

Unraveling local dynamic charge processes is vital for progress in diverse fields, from microelectronics to energy storage. This relies on the ability to map charge carrier motion across multiple length- and timescales and understanding how these processes interact with the inherent material heterogeneities. Towards addressing this challenge, we introduce high-speed sparse scanning Kelvin probe force microscopy, which combines sparse scanning and image reconstruction. This approach is shown to enable sub-second imaging (>3 frames per second) of nanoscale charge dynamics, representing several orders of magnitude improvement over traditional Kelvin probe force microscopy imaging rates. Bridging this improved spatiotemporal resolution with macroscale device measurements, we successfully visualize electrochemically mediated diffusion of mobile surface ions on a LaAlO 3 /SrTiO 3 planar device. Such processes are known to impact band-alignment and charge-transfer dynamics at these heterointerfaces. Furthermore, we monitor the diffusion of oxygen vacancies at the single grain level in polycrystalline TiO 2 . Through temperature-dependent measurements, we identify a charge diffusion activation energy of 0.18 eV, in good agreement with previously reported values and confirmed by DFT calculations. Together, these findings highlight the effectiveness and versatility of our method in understanding ionic charge carrier motion in microelectronics or nanoscale material systems.

25 ENERGY STORAGE↗

Nanoscale structural correlations in a model cuprate superconductor

Understanding the extent and role of inhomogeneity is a pivotal challenge in the physics of cuprate superconductors. While it is known that structural and electronic inhomogeneity is prevalent in the cuprates, it has proven difficult to disentangle compound-specific features from universally relevant effects. Here, in this study, we combine advanced neutron and x-ray diffuse scattering with numerical modeling to obtain insight into bulk structural correlations in HgBa 2 ⁢ CuO 4+δ . This cuprate exhibits a high optimal transition temperature of nearly 100 K, pristine charge-transport behavior, and a simple average crystal structure without long-range structural instabilities, and is therefore uniquely suited for investigations of intrinsic inhomogeneity. We uncover diffuse reciprocal-space patterns that correspond to prominent nanoscale correlations of atomic displacements perpendicular to the CuO 2 planes. The real-space nature of the correlations is revealed through three-dimensional pair distribution function analysis and complementary numerical refinement. We find that relative displacements of ionic and CuO 2 layers play a crucial role, and that the structural inhomogeneity is not directly caused by the presence of conventional point defects. The observed correlations are therefore intrinsic to HgBa 2 ⁢ CuO 4+δ , and thus likely important for the physics of cuprates more broadly. It is possible that the structural correlations are closely related to the unusual superconducting fluctuations and Mott-localization in these complex oxides. As advances in scattering techniques yield increasingly comprehensive data, the experimental and analysis tools developed here for large volumes of diffuse scattering data can be expected to aid future investigations of a wide range of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Local stress within a granular molecular solvent matrix, a mechanism for individual ion hydration

The activity of water is important in many natural systems and directed processes. To better understand water activity, the hydration of individual ions is approximated for a recently developed speciation-based solution model. Here a matrix of values for hydration of complete salts (anion and cation) was used as source data to calculate individual ion hydration values. Some electrolytes were excluded from calculations of individual ion hydration due to the prevalence of ion pairing and resulting liberation of water upon ion pairing. Evaluation of the periodic trends emerging from individual ion hydration values suggests that ionic volume, relative to the volume of molecular water, dictates both modeled hydration requirements (related to vapor-liquid-equilibrium-derived activity) as well as properties historically described by hydrodynamic radius (ion diffusion and solution viscosity). When water is conceptualized as a granular three-dimensional matrix, solutes may function as substitutional defects (dilation or constriction) inducing local compressive or tensile stress. These results have implications on the mechanism of the Hofmeister series and series reversals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DFT Accurate Interatomic Potential for Molten NaCl from Machine Learning

Molten alkali chloride salts are a critical component in concentrated solar power and nuclear applications. Despite their ubiquity, the extreme chemical reactivity of molten alkali chlorides at high temperatures has presented a significant challenge in characterizing atomic structures and dynamic properties experimentally. In this work, we we investigate molten NaCl by performing high temperature molecular dynamics simulations using a Gaussian Approximation Potential (GAP) trained on Density Functional Theory (DFT) datasets. Our GAP model, trained with a meager 1000 atomic configurations, arrives at near ab initio accuracy with a mean absolute error of 1.5 meV/atom thus enabling fast analysis of high temperature salt properties on large length (5000 ion pairs) and time (> 1ns) scales currently inaccessible to ab initio simulations. Calculated structure factors and diffusion constants from our GAP model simulations show excellent agreement with experiments. Our results indicate that GAP models are able to capture the many-body interactions required to accurately model ionic-systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ion Channels – From Structure to Electrophysiology and Back

A reliable way to establish whether our understanding of a channel is satisfactory is to reproduce its measured ionic conductance over a broad range of applied voltages in computer simulations. In molecular dynamics (MD), this can be done by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. Since this approach is computationally very expensive, we have developed a markedly more efficient alternative in which MD is combined with the electrodiffusion (ED) equation. In this approach, the assumptions of the ED equation can be rigorously tested, and the precision and consistency of the calculated conductance can be determined. We have demonstrated that the full current/voltage dependence and the underlying free energy profile for a simple channel can be reliably calculated from equilibrium or non-equilibrium MD simulations at a single voltage. Free energy profiles can be obtained from non-equilibrium simulations without a loss of accuracy even without the knowledge of diffusion coefficient. To carry out MD simulations, a structural model of a channel has to be assumed, which is an important constraint, considering that high-resolution structures are available for only very few simple channels. If the comparison of calculated ionic conductance with electrophysiological data is satisfactory, it greatly increases our confidence that the structure and the function are described sufficiently accurately. We examined the validity of the ED for several channels embedded in phospholipid membranes - four naturally occurring channels: trichotoxin, alamethicin, p7 from hepatitis C virus (HCV) and Vpu from the HIV-1 virus, a synthetic, hexameric channel, formed by a 21-residue peptide that contains only leucine and serine and a bacterial pentameric ligand-gated ion channel, GLIC. All these channels mediate transport of potassium and chloride ions. It was found that the ED equation is satisfactory for these systems. In some of them, experimental and calculated electrophysiological properties are in good agreement, whereas in others there are strong indications that the structural models are incorrect.

Pohorille, Andrew↗

A simulation-based model studying monoethanolamine and aprotic heterocyclic anion ionic liquid (AHA-IL) mixtures for carbon capture

Capturing of waste CO 2 , particularly from large industrial point sources, is necessary for either permanent storage or CO 2 utilization. Amine solvents used for carbon capture have certain drawbacks, particularly the high energy required by the reboiler for solvent regeneration and degradation of the amine with absorption/regeneration cycling. These challenges may be overcome using ionic liquids (ILs) that can be tailored for operating under typical flue gas stream conditions (8–10% CO 2 , 18–20% H 2 O, 2–3% O 2 , and 67–72% N 2 ). ILs, however, oftentimes suffer from high viscosities, leading to diffusion limitations when compared to amine solvents. This work models the combinations of aprotic heterocyclic anion ILs (AHA-IL) with monoethanolamine (MEA) and studies the carbon capture effectiveness of these mixtures. Furthermore, this study found that all AHA-IL/MEA mixtures exhibited 32 % lower enthalpy of reaction (compared to MEA) and 64 % lower viscosity (compared to AHA-IL), which are desirable characteristics for CO 2 capture solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase space densities in the Uranian plasma sheet

Based on its ionic composition, the magnetosphere of Uranus appears to contain little plasma of solar wind origin. If the plasma source is the atmosphere of the planet, both adiabatic convection with uniform losses and radial diffusion would imply negative radial gradients of the distribution function at fixed first and second adiabatic invariants. The analysis of LECP data in the outer magnetosphere reveals different gradients. Evidence is presented that ions may have been non-adiabaticaly heated during the fourth plasma sheet immersion although this is not unambiguously distinguishable from a source at large down tail distance.

Hammond, C. M.↗

Structure–transport correlations in Na 11 Sn 2 SbSe 12 and its sulfide solid solutions

We report a new ion conductor—Na 11 Sn 2 SbSe 12 —as a possible candidate for a solid catholyte in composite cathodes for all-solid state Na-batteries, which exhibits a room temperature ionic conductivity of 0.15 ± 0.03 mS cm −1 and an activation energy of 0.39 ± 0.02 eV. The sulfide solid solutions of Na 11 Sn 2 SbSe 12 , namely, Na 11 Sn 2 SbS 12−x Se x (x = 1 and 6), were also investigated through a combination of Rietveld refinement against powder x-ray diffraction data and electrochemical impedance spectroscopy to reveal the complex structure–property relationships governing ion transport in this class of materials. Meanwhile, broadening of the Na-ion diffusion pathways in Na 11 Sn 2 SbSe 12 is expected to facilitate Na-ion transport compared to the sulfide-rich member of the solid solution, the opposite holds: increasing the Se fraction in Na 11 Sn 2 SbS 12−x Se x leads to a little change in the activation energy but a reduction in the ionic conductivity. We ascribe this to the lowering of the prefactor, σ 0 , in the Arrhenius relationship with increasing lattice “softening” as a function of higher Se content.

Ramos, Erika P.↗