Pith. sign in

REVIEW 2 major objections 4 minor 22 references

Sequential acquisition of fluorescence signals with changing fluorophore concentrations. Multivariate Curve Resolution with time measurements assistance

T0 review · 2 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Fluorescence data can change while an excitation-emission matrix is being recorded and still be resolved into clean spectra and concentration profiles.

desk verdict Time-assisted MCR on pseudo-cubes is a plausible extension of the authors' PARAFAC work, but the validation inherits too much from the initialization to count as an independent demonstration. read the letter →

arxiv 2502.19431 v1 pith:ZG6DFIXZ submitted 2025-02-14 physics.med-ph physics.bio-phphysics.ins-det

Gabriel Siano , Sofia Mora , Agustina Schenone , Leonardo Giovanini This is my paper · ORCID
classification physics.med-phphysics.bio-phphysics.ins-det
keywords fluorescencespectroscopyexcitation-emissionmatricesmultivariatecurveresolutiontimemeasurementsmissingdataimputationpsCubeschromatography-fluorescencereactionkinetics
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to show that sequentially recorded fluorescence data remain analyzable even when fluorophore concentrations change during acquisition, by abandoning the idea that each excitation-emission matrix is one bilinear block. Instead, each emission scan is cut into short time fragments, arranged in time into partially filled three-way cubes, and modeled with Multivariate Curve Resolution assisted by measured times. The paper reports that for liquid-chromatography fluorescence data the lowest spectral similarity against reference spectra was 0.9963 for excitation and 0.9995 for emission, with a validation relative error of prediction of 1.76%, while for a kinetics dataset similarity stayed above 0.99 with a 3.67% error, despite 90.9% and 96.875% missing entries respectively. If true, this extends curve resolution to settings where the standard bilinearity assumption fails, and it makes calibration and quantitation possible in flow and reaction systems with ordinary spectrofluorometers.

What carries the argument

The central object is the psCube, a three-way array built from pseudo-EEMs: each pseudo-EEM is a matrix of the dimensions of the original EEM that contains only the experimental readings from one short fragment of one emission scan at one excitation wavelength, with everything else marked missing. Working with these cubes lets the model treat a slowly changing concentration as constant inside each fragment while still sampling the whole time course at many points. The argument is carried by three devices used together: time measurements that localize each fragment, expectation-maximization imputation of missing entries from the current bilinear model with the error computed only over experimental data, and smoothness constraints on the concentration profiles evaluated with each sample's measured times.

What would settle it

Run time-assisted MCR on synthetic excitation-emission data with a known ground-truth time course and randomly delete 90% of the entries, then compare the recovered concentration profile with the true one and with what expectation-maximization would predict from the initial spectra alone; if the recovered profile matches the truth much better than it matches the imputation prior, the imputation is not steering the solution, and if it does not, the reported resolution is partly an artifact.

Watch

Extended reading notes

Core claim

The central claim is that the apparent failure of bilinearity in EEMs acquired during chromatography or kinetics is not a property of the data but of the chosen structure. Because each single fluorescence reading is recorded in milliseconds, concentrations can be treated as constant within short intervals; the paper shows that grouping readings into pseudo-EEMs and localizing those fragments with measured times produces partially filled cubes whose reshaped matrices are compatible with bilinear MCR. With non-negativity, unimodality, and time-based smoothing, the resolved excitation and emission spectra match references closely, the concentration profiles are physically meaningful, and calibration predictions are comparable to those from higher-order models applied to the same data.

Load-bearing premise

The whole strategy stands on the assumption that imputing the missing entries from the current bilinear model at each iteration does not quietly shape the final answer, because the fitting error is computed only on the roughly 3 to 9 percent of readings that were actually measured.

Editorial extensions

If this is right

  • Conventional EEMs are not the only valid unit for fluorescence modeling; the same sequential scans can be read as high-resolution time profiles.
  • Time-assisted MCR can replace or complement PARAFAC in kinetic and chromatographic systems, giving each sample its own concentration profile instead of a shared shape.
  • A very high percentage of missing data is not inherently fatal: the structural bilinearity plus time localization carried models with 90.9% and 96.875% missing entries.
  • Calibration and prediction for an analyte that photodegrades during measurement can be done with the same data used for resolution, without dedicated reaction-monitoring hardware.
  • Time measurements must be used during the iterations, not only for graphing, because irregular inter-scan gaps are absorbed by time-based smoothing.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the psCube idea generalizes, any sequentially acquired spectroscopic data with fast individual reads and a slow drift, such as infrared or Raman reaction monitoring, could be reorganized the same way and resolved by bilinear MCR.
  • The reported mean recoveries slightly above and below 100%, 101.76% for LC-EEM and 96.33% for Kin-EEM, are consistent with the imputation step gently biasing concentrations toward the model, a point the paper leaves for future work.
  • A direct stress test would be to run the same workflow on synthetic dynamic EEMs with known true profiles and randomly delete 90% of the entries, which would isolate the effect of imputation from experimental noise.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The manuscript presents a Multivariate Curve Resolution (MCR) strategy for sequentially acquired fluorescence Excitation-Emission Matrices (EEMs) recorded while fluorophore concentrations change. The authors introduce 'psCubes', partially filled three-way arrays in which each short time segment of an emission scan is treated as a row with a high fraction of missing entries, and they use Expectation Maximization imputation plus time-based smoothing of concentration profiles. The method is applied to two datasets: LC-EEM (pyridoxine with interferents, chromatographic conditions) and Kin-EEM (diclofenac photodegradation kinetics). For both datasets, the authors report high spectral similarity with reference spectra (>0.99), low validation REP (1.76% and 3.67%), and physically meaningful concentration profiles, despite 90.9% and 96.875% missing data. The central claim is that time-assisted MCR resolves such data effectively.

Significance. If the claim were validated independently, the contribution would be practically useful: it would allow a bilinear MCR framework to handle data that are not conventionally bilinear, leveraging the time information to recover profiles from highly incomplete EEMs. The paper is transparent about several limitations, including the need for good initial estimates, and it reports detailed real experimental data. However, the evaluation does not currently provide independent evidence for the method's ability to resolve the data from the observed measurements alone, and the manuscript overstates the strength of the validation. The approach is potentially valuable as a refinement or downstream analysis tool when prior multilinear solutions are available, but that narrower claim requires reframing and additional experiments.

major comments (2)
  1. [§2.3.1 and §3.3.4] The central evaluation is not independent of the PARAFAC benchmark. Section 2.3.1 states that the MCR initial estimates and the initial missing-data imputations were taken from 4-way PARAFAC models of the same datasets from previous studies [5-7]. Section 3.3.4 then admits that without these prior results, or when they are used only for imputation or only for initialization, the MCR models 'faced convergence issues' and produced unsatisfactory results. Consequently, the reported REP values (1.76% for LC-EEM, 3.67% for Kin-EEM) and spectral similarities (>0.99) may reflect information inherited from the PARAFAC initialization rather than the ability of time-assisted MCR to resolve the data from the observed measurements. The comparison to 'processing the same data as cubes' is circular because the cube solutions are the MCR starting points. To support the abstract's claim that 'time-assisted MCR resolves such data effectively,' the authors should either (a) fit the prior PARAFAC models to the calibration set only and then apply MCR to the validation samples, (b) demonstrate concretely that MCR improves on the PARAFAC initialization in validation error or profile accuracy, or (c) show that the final MCR solution is insensitive to the choice of initialization. As written, this is a load-bearing gap in the validation.
  2. [§2.2] The treatment of missing data is insufficiently specified to rule out self-consistent artifacts. With 90.909% (LC-EEM) and 96.875% (Kin-EEM) of the entries missing, the imputed values vastly outnumber the experimental ones in the matrices used by the alternating least squares updates. The paper states that 'the minimization of errors in the fitting function during the iterations only considers the experimental data and never the imputed ones,' but the described procedure (EM imputation followed by bilinear modeling at each iteration, similar to the N-way Toolbox for PARAFAC) normally minimizes the complete-data loss, which includes the imputed entries. If the algorithm instead uses a weighted least-squares objective on the observed entries, that objective and the weighting scheme need to be stated explicitly. Without such clarification, and without a sensitivity analysis with respect to the initial imputation (e.g., random or perturbed initial imputations, or alternative starting profiles), the possibility remains that the resolved profiles are largely determined by the PARAFAC-based initialization rather than by the experimental information. This is particularly concerning because Section 3.3.4 reports convergence failures for other initialization schemes.
minor comments (4)
  1. [§2.2] The acronym nVarCLFC is used before its full expansion; the expansion 'number of Variables with Constant Local Fluorophore Concentrations' should be given at the first occurrence.
  2. [§3.3.4] The statement that 'these details have not been shown here and are part of future work' refers to the negative results with alternative initializations; for a journal article, at least a brief summary of the convergence failures and their characteristics should be included in the main text or supplementary information.
  3. [§3.1 / Table S2, S3] The abbreviation 'EJCR test OK' in the calibration tables is not defined in the text; please provide a definition or reference.
  4. [§3.3.2] The argument that the success of the chosen unfolding strategy 'suggests that such interactions do not actually occur' is an inference from a single successful model; a direct test of the EX-EM independence assumption, or at least a comparison with an alternative unfolding, would strengthen this point.

Circularity Check

2 steps flagged · score 5.0 of 10

MCR success rests on PARAFAC initialization from the same data; the reported spectral similarities and predictions are refinements of prior results rather than independent resolution.

  1. fitted input called prediction [Section 2.3.1, 'Previous results']
    "To obtain initial estimates for both datasets, the EX and EM spectra resolved with 4-way PARAFAC models assisted with time measurements were used. ... The initial estimates were also used to impute the missing data before the first iteration of each MCR model."

    The PARAFAC models come from the same datasets analyzed here [5-7]. With 90.909% (LC-EEM) and 96.875% (Kin-EEM) missing entries, the initial imputation fills nearly the entire data matrix with values generated from the PARAFAC solution, and each EM iteration re-imputes missing entries from the current bilinear model. The final MCR profiles are therefore a refinement of, not an independent extraction from, the PARAFAC solution. The reported spectral similarities and validation REP are thus partly inherited from the PARAFAC initialization rather than being standalone evidence for time-assisted MCR.

  2. self citation load bearing [Section 3.3.4, 'On estimating initial profiles and imputations in MCR']
    "Furthermore, both the initial imputations and the initial profile approximations were based on prior results from PARAFAC. Without these prior results, or using them only for data imputation or only for initial estimations (not both, and in their absence replaced by mean data or SIMPLISMA profiles, respectively), the MCR models faced convergence issues. The results obtained were not satisfactory, even when trilinearity constraints were applied."

    The paper's central claim that time-assisted MCR resolves such data effectively is explicitly conditional on prior PARAFAC solutions from the authors' own previous studies [5-7]. The paper states that without those prior results MCR fails. This is not an external, machine-checked, or independently reproduced benchmark; it is a self-citation chain that supplies both the initial profiles and the initial missing-data imputation. The MCR 'prediction' is therefore load-bearing on the authors' own prior model rather than being self-contained.

full rationale

The paper is transparent that MCR is initialized with PARAFAC results from the same datasets and that removing or weakening this initialization leads to convergence failure. That makes the central demonstration partly circular: the reported spectral similarities and the claimed effective resolution are obtained by refining prior PARAFAC solutions, not by extracting the profiles from the observed data independently. The dependence is not an equation-level identity, and there is still independent content in the per-sample c(t) profiles, the time-based smoothing, the comparison against external reference spectra, and the calibration predictions for a validation set. However, those external anchors are weakened by the fact that the initial spectral shapes already encode the PARAFAC answer, and the paper explicitly says the prior results are needed for both initialization and imputation. This warrants a moderate circularity score, between the 'some self-citation with independent content' level and the 'prediction reduces by construction' level.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central claim rests mainly on domain assumptions about time-local constancy and unbiased imputation of a very large missing-data fraction, plus modeling choices (nVarCLFC, smoothing, components) that are tuned on the datasets. No new physical entities are introduced; psEEMs and psCubes are organizational constructs.

free parameters (4)
  • nVarCLFC (LC-EEM) = 32 of 32 emission variables
    Chosen by running models with groupings from full spectra down to smaller fragments and selecting where improvements stop (Section 2.2); controls time resolution and the missing-data fraction.
  • nVarCLFC (Kin-EEM) = 16 of 64 emission variables
    Same selection procedure as LC-EEM; determines that each psEEM is one quarter of an emission scan (Section 2.2).
  • Smoothing parameters for c(t) profiles = Not reported numerically; reused from refs. 5-7
    Time-based smoothing of concentration profiles uses parameters tuned in previous studies on the same datasets (Section 2.3.1); these affect how strongly profiles are regularized.
  • Number of MCR components = 4 (LC-EEM), 3 (Kin-EEM)
    Component numbers are inherited from prior PARAFAC analyses of these data (Sections 3.1-3.2); the paper does not justify them independently.
assumptions (5)
  • domain assumption Fluorophore concentrations are constant within each psEEM time window (nVarCLFC consecutive readings).
    Central to the claim that each fragment is bilinear or trilinear; if concentrations change within a window, the model is misspecified. Sections 1 and 2.2.
  • domain assumption Expectation maximization imputation of 90.9% (LC) and 96.875% (Kin) missing entries using the current bilinear model at each iteration is unbiased.
    The model relies on imputed values for almost all entries; minimizing the fit only over experimental data is assumed to prevent imputation from dominating. Sections 2.2 and 3.1-3.2.
  • domain assumption EX and EM modes are independent of the c(t) mode within a psEEM, so unfolding a psEEM into a vector is valid.
    Justifies the bilinear MCR structure after reshaping; dependencies between modes would invalidate the unfolded matrix model. Section 3.3.2.
  • domain assumption Reference spectra used for similarity calculations are correct and the sample compositions are known.
    Similarity coefficients and validation predictions are computed against these references from pure standards and HPLC cross-checks. Sections 2.1 and 3.1-3.2.
  • domain assumption No inner filter effects, quenching, or other nonlinear interactions affect the fluorescence readings.
    The paper states concentrations were chosen to avoid these effects; if present, the linear model would fail. Section 2.1.2 and Conclusions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Sequential acquisition of fluorescence signals with changing fluorophore concentrations. Multivariate Curve Resolution with time measurements assistance." pith.science (2026). https://pith.science/paper/ZG6DFIXZ

@misc{pith2026250219431,
  author       = {Pith},
  title        = {Pith review of: Sequential acquisition of fluorescence signals with changing fluorophore concentrations. Multivariate Curve Resolution with time measurements assistance},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZG6DFIXZ}},
  note         = {Machine review of arXiv:2502.19431}
}
read the original abstract

Sequential registering of fluorescence signals in conventional Excitation-Emission Matrices (EEMs), followed by modeling based on multilinear properties of the data, requires stable fluorophore concentrations throughout the acquisition of each EEM. Rapid concentration changes, as seen in chromatography or certain kinetics, can disrupt the conventional bilinearity of EEMs. This deviation depends on the relative rates of concentration changes versus spectral scanning speeds. Although entire scans can be slow, individual data points are acquired almost instantaneously, maintaining linear dependencies. Within specific time intervals, concentrations can be assumed constant. By using time-based localization, partial EEM data can be organized into partially filled cubes that allow for the correct modeling of third-order data. Additionally, Multivariate Curve Resolution requires reshaping operations. Two third-order datasets were analyzed using a time-assisted MCR implementation, following a strategy similar to one reported for chromatographic data (LC-EEM) with Parallel Factor Analysis. The second dataset comes from Diclofenac reaction kinetics (Kin-EEM). The results suggest that time-assisted MCR resolves such data effectively. The solutions found were similar to those obtained when processing the same data as cubes. Predictions from calibration models based on MCR results were comparable to those obtained when deriving higher order models from the same data. High degrees of similarity were achieved between resolved and reference spectral profiles. Both chromatographic and kinetic profiles were accurate and physically meaningful. The proposed strategy highlights the potential of incorporating time-based localization to enhance the analysis of fluorescence data in dynamic systems.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

22 extracted references · 21 canonical work pages

  1. [1]

    Jaumot, A

    J. Jaumot, A. de Juan, R. Tauler, MCR-ALS GUI 2.0: New features and applications, Chemom. Intell. Lab. Syst. 140 (2015) 1–12. https://doi.org/10.1016/j.chemolab.2014.10.003

    work page doi:10.1016/j.chemolab.2014.10.003 2015

  2. [2]

    Bro, PARAFAC

    R. Bro, PARAFAC. Tutorial and applications, Chemom. Intell. Lab. Syst. 38 (1997) 149–171. https://doi.org/10.1016/S0169-7439(97)00032-4

    doi:10.1016/s0169-7439(97)00032-4 1997

  3. [3]

    Bortolato, V.A

    S.A. Bortolato, V.A. Lozano, A.M. de la Peña, A.C. Olivieri, Novel augmented parallel factor model for four-way calibration of high-performance liquid chromatography–fluorescence excitation–emission data, Chemom. Intell. Lab. Syst. 141 (2015) 1–11. https://doi.org/10.1016/j.chemolab.2014.11.013

    work page doi:10.1016/j.chemolab.2014.11.013 2015

  4. [4]

    Tomasi, R

    G. Tomasi, R. Bro, PARAFAC and missing values, Chemom. Intell. Lab. Syst. 75 (2005) 163–180. https://doi.org/10.1016/j.chemolab.2004.07.003

    work page doi:10.1016/j.chemolab.2004.07.003 2005

  5. [5]

    Siano, L

    G.G. Siano, L. Vera Candioti, L.L. Giovanini, Chemometric handling of 26 spectral-temporal dependencies for liquid chromatography data with online registering of excitation-emission fluorescence matrices, Chemom. Intell. Lab. Syst. 199 (2020) 103961. https://doi.org/10.1016/j.chemolab.2020.103961

    work page arXiv 2020

  6. [6]

    Mora, Sofía, A

    S. Mora, Sofía, A. Schenone, G. Siano, Estudio de fotodegradación del diclofecac mediante espectroscopia de fluorescencia con datos de tercer orden., Lic. Biotec., FBCB, Universidad Nacional del Litoral, 2022. http://www.bibliotecas.unl.edu.ar/vufind/Record/308805.19.f

    work page 2022

  7. [7]

    Siano, S

    G. Siano, S. Mora, A. Schenone, Third-order data generation and time measurement-assisted modeling from Diclofenac kinetics and acquisition of excitation-emission fluorescent signals., in: Book Abstr. XIX Chemom. Anal. Chem., Universidad Nacional del Litoral, Santa Fe, Argentina, 2024: p. 82. https://drive.google.com/file/d/1SSTQkDMDg1cuOZubEUZIP0t8rCvrPwtZ/view

    work page 2024

  8. [8]

    Montemurro, G.G

    M. Montemurro, G.G. Siano, M.J. Culzoni, H.C. Goicoechea, Automatic generation of photochemically induced excitation-emission-kinetic four-way data for the highly selective determination of azinphos-methyl in fruit juices, Sens. Actuators B Chem. 239 (2017) 397–404. https://doi.org/10.1016/j.snb.2016.08.033

    work page doi:10.1016/j.snb.2016.08.033 2017

Show all 22 references
  1. [9]

    Pellegrino Vidal, A.C

    R.B. Pellegrino Vidal, A.C. Olivieri, G.A. Ibañez, G.M. Escandar, Online Third-Order Liquid Chromatographic Data with Native and Photoinduced Fluorescence Detection for the Quantitation of Organic Pollutants in Environmental Water, ACS Omega 3 (2018) 15771–15779. https://doi.o...

    2018 doi

  2. [10]

    Siano, M

    G.G. Siano, M. Montemurro, M.R. Alcaráz, H.C. Goicoechea, Open-Source Assisted Laboratory Automation through Graphical User Interfaces and 3D Printers: Application to Equipment Hyphenation for Higher-Order Data Generation, Anal. Chem. 89 (2017) 10667–10672. https://doi.org/10....

    2017 doi

  3. [11]

    Escandar, A.C

    G.M. Escandar, A.C. Olivieri, A road map for multi-way calibration models, Analyst 142 (2017) 2862–2873. https://doi.org/10.1039/C7AN00822H

    2017 doi

  4. [12]

    Gómez, M

    V. Gómez, M. Miró, M.P. Callao, V. Cerdà, Coupling of Sequential Injection Chromatography with Multivariate Curve Resolution-Alternating Least-Squares for Enhancement of Peak Capacity, Anal. Chem. 79 (2007) 7767–7774. https://doi.org/10.1021/ac071202h

    2007 doi

  5. [13]

    Tauler, Multivariate curve resolution of multiway data using the multilinearity constraint, J

    R. Tauler, Multivariate curve resolution of multiway data using the multilinearity constraint, J. Chemom. 35 (2021) e3279. https://doi.org/10.1002/cem.3279

    2021 doi

  6. [14]

    Zhang, R

    X. Zhang, R. Tauler, Flexible Implementation of the Trilinearity Constraint in Multivariate Curve Resolution Alternating Least Squares (MCR-ALS) of Chromatographic and Other Type of Data, Mol. Basel Switz. 27 (2022) 2338. https://doi.org/10.3390/molecules27072338

    2022 doi

  7. [15]

    Schneide, R

    P.-A. Schneide, R. Bro, N.B. Gallagher, Shift-invariant tri-linearity—A new model for resolving untargeted gas chromatography coupled mass spectrometry data, J. Chemom. 37 (2023) e3501. https://doi.org/10.1002/cem.3501

    2023 doi

  8. [16]

    Zhang, R

    X. Zhang, R. Tauler, Flexible Trilinearity Alignment (FTA) and Shift Invariant Transformation (SIT) Constraints in Three-Way Multivariate Curve Resolution Data Analysis, J. Chemom. 38 (2024) e3581. https://doi.org/10.1002/cem.3581

    2024 doi

  9. [17]

    Schneide, N.B

    P.-A. Schneide, N.B. Gallagher, R. Bro, Shift invariant soft trilinearity: Modelling shifts and shape changes in gas-chromatography coupled mass spectrometry, Chemom. Intell. Lab. Syst. 251 (2024) 105155. https://doi.org/10.1016/j.chemolab.2024.105155

    2024

  10. [20]

    Andersson, R

    C.A. Andersson, R. Bro, The N-way Toolbox for MATLAB, Chemom. Intell. Lab. Syst. 52 (2000) 1–4. https://doi.org/10.1016/S0169-7439(00)00071-X

    2000 doi

  11. [21]

    Models, Algorithms, and Applications, Ph.D

    Bro, Rasmus, Multi-way Analysis in the Food Industry. Models, Algorithms, and Applications, Ph.D. Thesis, University of Amsterdam, 1998

    1998

  12. [22]

    Timmerman, H.A.L

    M.E. Timmerman, H.A.L. Kiers, Three-way component analysis with smoothness constraints, Comput. Stat. Data Anal. 40 (2002) 447–470. 27 https://doi.org/10.1016/S0167-9473(02)00059-2. [23] A.C. Olivieri, H.L. Wu, R.Q. Yu, MVC3: A MATLAB graphical interface toolbox for third-orde...

    2002 doi

  13. [24]

    Zhang, G

    N. Zhang, G. Liu, H. Liu, Y. Wang, Z. He, G. Wang, Diclofenac photodegradation under simulated sunlight: Effect of different forms of nitrogen and Kinetics, J. Hazard. Mater. 192 (2011) 411–418. https://doi.org/10.1016/j.jhazmat.2011.05.038

    2011 doi

  14. [25]

    Windig, Jean

    Willem. Windig, Jean. Guilment, Interactive self-modeling mixture analysis, Anal. Chem. 63 (1991) 1425–1432. https://doi.org/10.1021/ac00014a016. 28 Supplementary Information Table of contents Table S1 Conventional MCR model. Qualitative and quantitative descriptors for the LC...

    1991 doi

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.