Search NASA⌕ Search

SEARCH · Search NASA

Results for “MERLIN”

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 19 records

Enabling machine learning-ready HPC ensembles with Merlin

With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computing (HPC) environment. Here, we present Merlin, a workflow framework to enable large ML-friendly ensembles of scientific HPC simulations. By augmenting traditional HPC with distributed compute technologies, Merlin aims to lower the barrier for scientific subject matter experts to incorporate ML into their analysis. As a producer–consumer workflow model, Merlin enables multi-machine, cross-batch job, dynamically allocated yet persistent workflows capable of utilizing surge-compute resources. Key features of Merlin are a flexible HPC-centric interface, low per-task overhead, multi-tiered fault recovery, and a hierarchical sampling algorithm that allows for $\mathscr{O}$(N) task execution and $\mathscr{O}$(N ln N) task queuing to ensembles of millions of tasks. In addition to Merlin’s design, we test the algorithm’s performance in an HPC center and demonstrate the ability to enqueue 40 million simulations in 100 s, with a 30 millisecond per-task overhead that is independent of ensemble size. Finally, we describe some example applications that Merlin has enabled on leadership-class HPC resources, such as the ML-augmented optimization of nuclear fusion experiments and the calibration of infectious disease models to study the progression of and possible mitigation strategies for COVID-19.

97 MATHEMATICS AND COMPUTING↗

merlin , an improved framework for the reconstruction of high-quality genome-scale metabolic models

Abstract Genome-scale metabolic models have been recognised as useful tools for better understanding living organisms’ metabolism. merlin (https://www.merlin-sysbio.org/) is an open-source and user-friendly resource that hastens the models’ reconstruction process, conjugating manual and automatic procedures, while leveraging the user's expertise with a curation-oriented graphical interface. An updated and redesigned version of merlin is herein presented. Since 2015, several features have been implemented in merlin, along with deep changes in the software architecture, operational flow, and graphical interface. The current version (4.0) includes the implementation of novel algorithms and third-party tools for genome functional annotation, draft assembly, model refinement, and curation. Such updates increased the user base, resulting in multiple published works, including genome metabolic (re-)annotations and model reconstructions of multiple (lower and higher) eukaryotes and prokaryotes. merlin version 4.0 is the only tool able to perform template based and de novo draft reconstructions, while achieving competitive performance compared to state-of-the art tools both for well and less-studied organisms.

Capela, João (ORCID:0000000212352922)↗

Medium Energy Electron Flux in Earth's Outer Radiation Belt (MERLIN): A Machine Learning Model

The radiation belts of the Earth, filled with energetic electrons, comprise complex and dynamic systems that pose a significant threat to satellite operation. While various models of electron flux both for low and relativistic energies have been developed, the behavior of medium energy (120–600 keV) electrons, especially in the MEO region, remains poorly quantified. At these energies, electrons are driven by both convective and diffusive transport, and their prediction usually requires sophisticated 4D modeling codes. In this paper, we present an alternative approach using the Light Gradient Boosting (LightGBM) machine learning algorithm. The Medium Energy electRon fLux In Earth's outer radiatioN belt (MERLIN) model takes as input the satellite position, a combination of geomagnetic indices and solar wind parameters including the time history of velocity, and does not use persistence. MERLIN is trained on >15 years of the GPS electron flux data and tested on more than 1.5 years of measurements. Tenfold cross validation yields that the model predicts the MEO radiation environment well, both in terms of dynamics and amplitudes o f flux. Evaluation on the test set shows high correlation between the predicted and observed electron flux (0.8) and low values of absolute error. The MERLIN model can have wide space weather applications, providing information for the scientific community in the form of radiation belts reconstructions, as well as industry for satellite mission design, nowcast of the MEO environment, and surface charging analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-scale VLBI observations of the candidate host galaxy of GRB 200716C

We present the discovery and the subsequent follow up of radio emission from SDSS J130402.36+293840.6 (J1304+2938), the candidate host galaxy of the gamma-ray burst (GRB) GRB 200716C. The galaxy is detected in the RACS (0.89 GHz), the NVSS, the Apertif imaging survey, and the FIRST (1.4 GHz), the VLASS (3 GHz), and in public LOFAR (130–170 MHz), WISE (3.4–22 μm), and SDSS (z, i, r, g, u filters) data. The luminosity inferred at 1.4 GHz is (5.1 ± 0.2) × 10 30 erg s -1 Hz -1 . To characterise the emission and distinguish between different components within the galaxy, we performed dedicated, high-sensitivity and high-resolution observations with the European VLBI Network (EVN) + e-MERLIN at 1.6 and 5 GHz. We did not detect any emission from a compact core, suggesting that the presence of a radio-loud active galactic nucleus (AGN) is unlikely, and therefore we ascribe the emission observed in the public surveys to star-forming regions within the galaxy. We confirm and refine the redshift estimate, z = 0.341 ± 0.004, with a dedicated Telescopio Nazionale Galileo (TNG) spectroscopic observation. Finally, we compiled a list of all the known hosts of GRB afterglows detected in radio and computed the corresponding radio luminosity: if GRB 200716C belongs to J1304+2938, this is the third most radio-luminous host of a GRB, implying one of the highest star-formation rates (SFRs) currently known, namely SFR ~ 324±61 M ⊙ yr -1 . On the other hand, through the analysis of the prompt emission light curve, recent works suggest that GRB 200716C might be a short-duration GRB located beyond J1304+2938 and gravitationally lensed by an intermediate-mass black hole (IMBH) hosted by the galaxy. Neither the public data nor our Very Long Baseline Interferometry (VLBI) observations can confirm or rule out the presence of an IMBH acting as a (milli-)lens hosted by the galaxy, a scenario still compatible with the set of radio observations presented in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

TranSyT , an innovative framework for identifying transport systems

The importance and rate of development of genome-scale metabolic models have been growing for the last few years, increasing the demand for software solutions that automate several steps of this process. However, since TRIAGE’s release, software development for the automatic integration of transport reactions into models has stalled. Here, in this paper, we present the Transport Systems Tracker (TranSyT). Unlike other transport systems annotation software, TranSyT does not rely on manual curation to expand its internal database, which is derived from highly curated records retrieved from the Transporters Classification Database and complemented with information from other data sources. TranSyT compiles information regarding transporter families and proteins, and derives reactions into its internal database, making it available for rapid annotation of complete genomes. All transport reactions have GPR associations and can be exported with identifiers from four different metabolite databases. TranSyT is currently available as a plugin for merlin v4.0 and an app for KBase.

59 BASIC BIOLOGICAL SCIENCES↗

Young but fading radio sources: searching for remnants among compact steep-spectrum radio sources

The incidence of young but fading radio sources provides important information on the life cycle of radio emission in radio-loud active galactic nuclei. Despite its importance for constraining the models of radio source evolution, there are no systematic studies of remnants in complete samples of young radio sources. We report results of the study of 18 compact steep-spectrum (CSS) radio sources, selected from the statistically complete B3-VLA (Very Large Array) CSS sample, characterized by a steep optically thin spectrum (α ≥ 1.0) and no core detection in earlier studies. Our deep multifrequency VLA, pc-scale Very Long Baseline Array, and enhanced Multi Element Remotely Linked Interferometer Network (e-MERLIN) observations allowed us to locate the core component in 10 objects. In three CSS sources, there is no clear evidence of present-time active regions, suggesting that they are likely in a remnant phase. Among sources with core detection, we find three objects that have no clear active regions (hotspots) at the edges of the radio structure, suggesting that the radio emission may have just restarted. Our results support a power-law distribution of the source ages, although the poor statistics prevents us from setting solid constraints on the percentage of remnants and restarted sources in subpopulations of radio sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Mutagenesis of Human Cytomegalovirus Glycoprotein L Disproportionately Disrupts gH/gL/gO over gH/gL/pUL128-131

Herein, cell-free and cell-to-cell spread of herpesviruses involves a core fusion apparatus comprised of the fusion protein glycoprotein B (gB) and the regulatory factor gH/gL. The human cytomegalovirus (HCMV) gH/gL/gO and gH/gL/pUL128-131 facilitate spread in different cell types. The gO and pUL128-131 components bind distinct receptors, but how the gH/gL portions of the complexes functionally compare is not understood. We previously characterized a panel of gL mutants by transient expression and showed that many were impaired for gH/gL-gB-dependent cell-cell fusion but were still able to form gH/gL/pUL128-131 and induce receptor interference. Here, the gL mutants were engineered into the HCMV BAC clones TB40/e-BAC4 (TB), TR, and Merlin (ME), which differ in their utilization of the two complexes for entry and spread. Several of the gL mutations disproportionately impacted gH/gL/gO-dependent entry and spread over gH/gL/pUL128-131 processes. The effects of some mutants could be explained by impaired gH/gL/gO assembly, but other mutants impacted gH/gL/gO function. Soluble gH/gL/gO containing the L201 mutant failed to block HCMV infection despite unimpaired binding to PDGFRα, indicating the existence of other important gH/gL/gO receptors. Another mutant (L139) enhanced the gH/gL/gO-dependent cell-free spread of TR, suggesting a “hyperactive” gH/gL/gO. Recently published crystallography and cryo-electron microscopy studies suggest structural conservation of the gH/gL underlying gH/gL/gO and gH/gL/pUL128-131. However, our data suggest important differences in the gH/gL of the two complexes and support a model in which gH/gL/gO can provide an activation signal for gB.

59 BASIC BIOLOGICAL SCIENCES↗

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Barge Site - Avian Radar System / Derived Data

This is a combined data set of 67,410 bird/bat tracks from an avian radar system deployed on a research barge (MERLIN True3D, DeTect, Panama City, Florida, USA) and concurrent wind measurements from two scanning lidars (WindCube v2.1, Vaisala, Vantaa, Finland, and Halo XR+, Halo Photonics, Lannion, France). The research barge (16.5 m x 61 m) was deployed as part of the Wind Forecast Improvement Project (WFIP-3) off the northeast coast of the United States south of Massachusetts (40.9 deg N, 70.79 deg W). This data set comprises 5 weeks of data between August 27th 2024 and September 27th 2024. Radar data were provided by DeTect and Lidar data were accessed through the Wind Data Hub (wfip3/barg.WINDPROF.z01.a0) The data have been filtered and sorted into two size groups ("big" and "small") based on a clustering approach. See Snortland, A., Clerc, J., Hein, C., & Cotter, E. (2025). Wind as Driver of Bird and Bat Abundance, Flight Direction, Altitude, and Speed on the North Atlantic Shelf. arXiv preprint arXiv:2511.14983 for complete details. Data are provided in 2 files: "Birds" and "Birds_hourly" Birds: This file contains information about each of the 67,410 flying animal tracks detected by the radar during the data collection period, including parameters measured by the radar and wind information interpolated from the lidar wind measurements. We note that the raw radar dataset contained 301,618 tracks; tracks in this processed dataset were filtered based on the requirements described in Snortland et al. (2025). Birds_hourly: This file contains timeseries of the number of tracks detected per hour over the course of the data collection period, including wind conditions and sun position for each hour. These data were used for generalized additive modeling in Snortland et al. (2025).

17 WIND ENERGY↗

Pyrolysis Molecular Beam Mass Spectrometry_Analysis_of_Natural_Variants_of_Poplulus_Trichocarpa_Leaves

Select leaves from natural variants of Poplar (Populus Trichocarpa) grown in a greenhouse at Oak Ridge National Laboratory were analyzed by Pyrolysis-Molecular Beam Mass Spectrometry (Py-MBMS). Leaves were harvested, cryomilled and kept frozen until analysis. Py-MBMS analysis was conducted using approximately 4 mg of biomass and each sample was analyzed in duplicate. A Frontier PY2020 unit pyrolyzed samples at 500°C for 30 s in 80 µL deactivated stainless steel cups. An Extrel Super-Sonic MBMS Model Max 1000 was used to collect mass spectral data fromm/z30 to 450 at 17 eV and processed using Merlin Automation software (V3). Spectral ion intensities were normalized to the total ion chromatogram signal for each sample for analysis of spectral variance. Lignin content (wt %) was estimated based on relative responses from standards of known Klason lignin content using mean-normalized ion intensities ofm/z120, 124 (G), 137 (G), 138 (G), 150 (G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S) where G indicates guaiacyl-derived ions, S indicates syringyl-derived ions, and other ions either derive from other lignin monomers or multiple sources. Ratios of S and G lignin monomer units (S/G) were obtained by dividing the sum of S-based ions by the sum of G-based ions using mean-normalized ion intensities.

CBI↗

Pyrolysis_Molecular_Beam_Mass_Spectrometry_Analysis_of_Specific_Switchgrass_Genotypes

Select natural variant switchgrass genotypes grown in Tifton, GA were analyzed by Pyrolysis-Molecular Beam Mass Spectrometry (Py-MBMS). Biomass was harvested, milled, several genotypes were analyzed with and without being destarched and extracted with ethanol prior to analysis (indicated with -DE if destarched and extracted). Py-MBMS analysis was conducted using approximately 4 mg of biomass and each sample was analyzed in duplicate. A Frontier PY2020 unit pyrolyzed samples at 500°C for 30 s in 80 µL deactivated stainless steel cups. An Extrel Super-Sonic MBMS Model Max 1000 was used to collect mass spectral data fromm/z30 to 450 at 17 eV and processed using Merlin Automation software (V3). Spectral ion intensities were normalized to the total ion chromatogram signal for each sample for analysis of spectral variance. Lignin content (wt %) was estimated based on relative responses from standards of known Klason lignin content using mean-normalized ion intensities ofm/z120, 124 (G), 137 (G), 138 (G), 150 (G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S) where G indicates guaiacyl-derived ions, S indicates syringyl-derived ions, and other ions either derive from other lignin monomers or multiple sources. Ratios of S and G lignin monomer units (S/G) were obtained by dividing the sum of S-based ions by the sum of G-based ions using mean-normalized ion intensities.

CBI↗

Pyrolysis_Molecular_Beam_Mass_Spectrometry_Analysis_of_hybrid_cross_of_Populus_tremula_x_P_alba_717-1B4_and_overexpression_of_a_lectin_receptor-like_kinase_(PtLecRLK1)

Stem tissues from the hybrid poplarPopulus tremula × P. albaclone 717-1B4 and from lectin receptor-like kinase overexpression lines PP7 and PP19 were individually colonized with the ectomycorrhizal fungiLaccaria bicolorstrain S238N,Hyaloscypha finlandicastrain PMI746, orUmbelopsis vinaceastrain PMI3018, as well as with a mixed fungal inoculum; non-inoculated plants served as controls. Plants were grown in a greenhouse at Oak Ridge National Laboratory and harvested in January 2025. Stem samples were analyzed using Pyrolysis–Molecular Beam Mass Spectrometry (Py-MBMS). Stems were harvested, debarked, dried, milled, destarched and ethanol extracted prior to analysis. Py-MBMS analysis was conducted using approximately 4 mg of wood from biomass and each sample was analyzed in duplicate. A Frontier PY2020 unit pyrolyzed samples at 500°C for 30 s in 80 µL deactivated stainless steel cups. An Extrel Super-Sonic MBMS Model Max 1000 was used to collect mass spectral data fromm/z30 to 450 at 17 eV and processed using Merlin Automation software (V3). Spectral ion intensities were normalized to the total ion chromatogram signal for each sample for analysis of spectral variance. Lignin content (wt %) was estimated based on relative responses from standards of known Klason lignin content using mean-normalized ion intensities ofm/z120, 124 (G), 137 (G), 138 (G), 150 (G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S) where G indicates guaiacyl-derived ions, S indicates syringyl-derived ions, and other ions either derive from other lignin monomers or multiple sources. Ratios of S and G lignin monomer units (S/G) were obtained by dividing the sum of S-based ions by the sum of G-based ions using mean-normalized ion intensities.

CBI↗

The Energetics of the Central Engine in the Powerful Quasar 3C 298

The compact steep-spectrum radio source 3C 298 (redshift of 1.44) has the largest 178 MHz luminosity in the Third Cambridge Revised Catalogue (3CR); its radio lobes are among the most luminous in the universe. The plasma state of the radio lobes is modeled with the aid of interferometric radio observations (in particular, the new Low Frequency Array observation and archival MERLIN data) and archival single-station data. It is estimated that the long-term time-averaged jet power required to fill these lobes with leptonic plasma is $\overline{Q}\,\approx 1.28\pm 0.51\times {10}^{47}\,\mathrm{erg}\,{{\rm{s}}}^{-1}$, rivaling the largest time-averaged jet powers from any quasar. Supporting this notion of extraordinary jet power is a 0.5–10 keV luminosity of ≈ 5.2 × 10 46 erg s –1 , comparable to luminous blazars, yet there is no other indication of strong relativistic beaming. We combine two new high signal-to-noise ratio optical spectroscopic observations from the Hobby-Eberly Telescope with archival Hubble Space Telescope, Two Micron All Sky Survey, and Galaxy Evolutionary Explorer data to compute a bolometric luminosity from the accretion flow of L bol ≈ 1.55 ± 0.15 × 10 47 erg s –1 . The ratio, $\overline{Q}/{L}_{\mathrm{bol}}\approx 1$, is the approximate upper limit for quasars. Characteristic of a large $\overline{Q}/{L}_{\mathrm{bol}}$, we find an extreme-ultraviolet (EUV) spectrum that is very steep (the "EUV deficit" of powerful radio quasars relative to radio-quiet quasars), and this weak ionizing continuum is likely a contributing factor to the relatively small equivalent widths of the broad emission lines in this quasar.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of the OECD/NEA PWR MOX/UO{sub 2} core transient benchmark using a coupled whole core pin-by-pin route in WIMS

This paper investigates the new CAMELOT coupling route developed in WIMS, which uses prepared tabulated cross-sections combined with the flux solver MERLIN module and the integrated sub-channel thermal hydraulics solver ARTHUR module to solve for the coupled core state. This coupling route is applied to the OECD/NEA pressurized water reactor (PWR) core transient benchmark exercise, involving 3-D neutronics and thermal-hydraulics analysis of a full-sized PWR core featuring MOX fuel at various core states, as well as simulating a control rod ejection transient. Analysis of the results for each part of the benchmark demonstrate that the high-resolution pin-by-pin analysis used in WIMS closely matches the other participants in the benchmark, with differences of between 20-200 pcm for k-effective, 1-2.5% RMS difference for assembly-averaged powers, 2-8% difference in rod worth calculations, < 200 ppm for calculations of critical boron concentration, zero difference in delayed neutron fraction calculation and a very close match for power response to the rod ejection transient. This is the first application of the WIMS-CAMELOT approach for the transient analysis of a full-sized reactor and as such further improvements such as reducing computational cost through parallelization, optimisation and memory reduction, are currently in active development. This paper demonstrates the accuracy and flexibility of the CAMELOT coupling route, providing a straightforward and consistent user image that can be easily used for a variety of modelling problems involving coupled neutronics and thermal-hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Probing axionlike particles with multimessenger observations of neutron star mergers

Axion-like particles (ALPs) can be copiously produced in binary neutron star (BNS) mergers through nucleon-nucleon bremsstrahlung if the ALP-nucleon couplings 𝑔 𝑎⁢𝑁 are sizable. The ALP-photon coupling 𝑔 𝑎⁢𝛾 may trigger conversions of ultralight ALPs into photons in the magnetic fields of the merger remnant and of the Milky Way. This effect would lead to a potentially observable short gamma-ray signal, in coincidence with the gravitational-wave signal produced during the merging process. This event could be detected through multimessenger observation of BNS mergers employing the synergy between gravitational-wave detectors and gamma-ray telescopes. Here, we study the sensitivity of current and proposed MeV gamma-ray experiments to detect such a signal. As an explicit example, we consider ALP couplings related as in the Kim-Shifman-Vainshtein-Zakharov axion model, and show that in this case the proposed instruments can reach a sensitivity down to 𝑔 𝑎⁢𝛾 ≳ few ×10 −13 GeV −1 for 𝑚 𝑎 ≲ 10 −9 eV, comparable with the SN 1987A limit.

Axion-like particles↗

Current–Induced Ion Concentration Polarization at a Perfect Ion–Exchange Patch in an Infinite Insulating Wall

This research examines theoretically the ion concentration polarization, ion fluxes, and electrostatic fields near an ion exchange patch in the wall of an electrified fluidic channel. These phenomena are important in related microfluidic ion-preconcentration systems. Under an electric field, counterions enter the ion-exchange patch at one side and leave at the other, with salt depletion occurring near the entrance and accumulation near the exit. Furthermore, the high patch conductivity and the concentration profiles lead to local electric field perturbations that may facilitate preconcentration. This study includes analytical expressions of ion concentrations and electrochemical potentials at small to moderate electric fields, as well as numerical simulations. Additionally, a simple matrix of polynomial coefficients (obtained via fitting of numerical data) enables analytical calculation of the 2-dimensional concentration profiles at all electric fields within the range investigated in the numerical simulations. This is possible because a single dimensionless parameter controls this problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗