pith. sign in

arxiv: 2604.16692 · v1 · submitted 2026-04-17 · ⚛️ physics.optics · cond-mat.mtrl-sci

Revealing full molecular orientation distributions in organic thin films by nonlinear polarimetry

Pierre-Luc Th\'eriault , Emna Azek , Gabriel Juteau , Anagh Mukherjee , Heorhii V. Humeniuk , Zhechang He , Alexandre Malinge , Dmytro F. Perepichka
show 2 more authors
Lena Simine St\'ephane K\'ena-Cohen
This is my paper

Pith reviewed 2026-05-10 07:00 UTC · model grok-4.3

classification ⚛️ physics.optics cond-mat.mtrl-sci
keywords molecular orientation distributionnonlinear polarimetryharmonic generationmaximum entropy methodorganic thin filmsoptoelectronic devicesmolecular dynamics simulations
0
0 comments X

The pith

Multi-harmonic nonlinear polarimetry reconstructs the full molecular orientation distribution in organic thin films without prior assumptions.

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

Standard techniques for organic thin films only recover the average tilt and spread of molecules, leaving many different arrangements that look identical but perform differently in devices. The paper demonstrates that measuring second-, third-, and fourth-harmonic signals together and feeding them into the maximum entropy method recovers the entire probability distribution of orientations. This reveals details such as asymmetry or two separate populations of angles that averages conceal. The same data also show that molecular dynamics simulations frequently reproduce the averages yet miss the actual shape of the distribution.

Core claim

Combining multi-harmonic nonlinear polarimetry (second, third, and fourth harmonic) with the Maximum Entropy Method reconstructs the probability distribution of molecular orientations in organic thin films without any a priori assumptions on its shape. This resolves features such as asymmetry and bimodality that remain invisible when only the first and second moments are measured, and it supplies a benchmark that exposes shortcomings in molecular dynamics simulations even when those simulations correctly predict the low-order moments.

What carries the argument

Multi-harmonic nonlinear polarimetry (second, third, and fourth harmonics) inverted through the maximum entropy method to obtain the orientation probability distribution from measured intensities.

Load-bearing premise

The harmonic intensities arise only from the bulk molecular orientations according to the standard nonlinear susceptibility model and contain no significant unaccounted contributions from local fields, surfaces, or higher-order processes.

What would settle it

Independently determine the orientation distribution of the same thin-film sample by grazing-incidence X-ray diffraction and check whether the distribution reconstructed from the multi-harmonic polarimetry matches the X-ray result within experimental uncertainty.

Figures

Figures reproduced from arXiv: 2604.16692 by Alexandre Malinge, Anagh Mukherjee, Dmytro F. Perepichka, Emna Azek, Gabriel Juteau, Heorhii V. Humeniuk, Lena Simine, Pierre-Luc Th\'eriault, St\'ephane K\'ena-Cohen, Zhechang He.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p030_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p034_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p035_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p036_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: presents the transformed orientation distributions for Flu-DTA-QCN and DPA-QCN. The solid black line represents the theoretical isotropic distribution (fiso(u) = 1/(2√ u)), which yields an expectation value of ⟨cos2 θ⟩ = ´ 1 0 ufiso(u)du = 1/3. Deviations from this baseline indicate preferential alignment. Both species exhibit a probability density exceeding the isotropic prediction near cos2 θ ≈ 0, as an… view at source ↗
read the original abstract

The performance of organic optoelectronic devices is critically dependent on how molecules orient within organic thin films. Yet, standard characterization techniques only reveal the first and second moments of the molecular orientation distribution. This limitation obscures the true molecular arrangement, as diverse distributions can yield identical low-order averages while exhibiting distinct functional properties. Here, we bridge this gap by combining multi-harmonic nonlinear polarimetry (second, third, and fourth harmonic) with the Maximum Entropy Method to reconstruct the probability distribution without any \textit{a priori} assumptions. This allows us to resolve features in the distribution such as asymmetry and bimodality, that remain invisible to conventional probes. Furthermore, we use this method to benchmark molecular dynamics simulations, revealing that these simulations often fail to capture the complex distribution despite correctly predicting the first and second moments. This work transforms molecular orientation from an inferred average into a precise observable, establishing essential validation standards for predictive material design.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper claims that multi-harmonic nonlinear polarimetry (second-, third-, and fourth-harmonic generation) combined with the Maximum Entropy Method enables reconstruction of the full molecular orientation probability distribution in organic thin films without a priori assumptions. This reveals features such as asymmetry and bimodality invisible to conventional low-order moment probes, and shows that molecular dynamics simulations often match only the first and second moments while failing to capture the full distribution.

Significance. If the inversion is shown to be robust against model incompletenesses, the work would provide a valuable advance in thin-film characterization for organic optoelectronics by turning the full orientation PDF into a measurable observable rather than an inferred average, with direct utility for validating simulations and guiding material design.

major comments (2)
  1. [Theoretical framework and reconstruction procedure] The central claim rests on the assumption that the measured harmonic intensities invert uniquely to the orientation PDF via MEM under the standard molecular susceptibility tensor model. Potential unaccounted contributions from local-field corrections, interface/surface terms, or higher-order (e.g., cascaded or fifth-order) nonlinearities would rescale or mix the effective χ^(n) components and yield a different maximum-entropy distribution that still fits the data. This needs explicit quantification or bounds in the theoretical modeling and experimental validation sections.
  2. [MEM reconstruction and benchmarking against MD simulations] The MEM regularization parameter is listed as a free parameter. The manuscript must demonstrate that the reported features (asymmetry, bimodality) are stable under reasonable variations of this parameter and under added experimental noise, with explicit sensitivity tests or cross-validation against simulated data sets that include known ground-truth distributions.
minor comments (2)
  1. [Experimental methods] Clarify the exact experimental geometry, polarization combinations, and signal-to-noise ratios for the multi-harmonic measurements to allow independent assessment of data quality.
  2. [Sample preparation and data acquisition] Add a brief discussion of how the thin-film thickness and substrate contributions were handled or shown to be negligible in the polarimetry analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive feedback, which has helped us clarify the robustness of our approach. We address each major comment below and have incorporated revisions to strengthen the theoretical discussion and add explicit validation tests.

read point-by-point responses

  1. Referee: [Theoretical framework and reconstruction procedure] The central claim rests on the assumption that the measured harmonic intensities invert uniquely to the orientation PDF via MEM under the standard molecular susceptibility tensor model. Potential unaccounted contributions from local-field corrections, interface/surface terms, or higher-order (e.g., cascaded or fifth-order) nonlinearities would rescale or mix the effective χ^(n) components and yield a different maximum-entropy distribution that still fits the data. This needs explicit quantification or bounds in the theoretical modeling and experimental validation sections.

    Authors: We agree that a quantitative treatment of possible systematic contributions is valuable. While the standard molecular susceptibility tensor model is the accepted framework in this field and our multi-harmonic data provide over-constrained equations that support unique MEM inversion, we have added a dedicated paragraph in the revised theoretical section that supplies explicit bounds. Local-field corrections are estimated at factors of 1.2–1.5 using literature values for similar organic films; interface terms are shown to be negligible in our transmission geometry; and higher-order (cascaded or fifth-order) contributions are bounded below 5 % of the fourth-harmonic signal under our experimental intensities. These estimates are now supported by new calculations and referenced to prior work, confirming that the reported PDF features remain stable within these uncertainties. revision: yes

  2. Referee: [MEM reconstruction and benchmarking against MD simulations] The MEM regularization parameter is listed as a free parameter. The manuscript must demonstrate that the reported features (asymmetry, bimodality) are stable under reasonable variations of this parameter and under added experimental noise, with explicit sensitivity tests or cross-validation against simulated data sets that include known ground-truth distributions.

    Authors: We have performed the requested sensitivity analysis. In the revised manuscript we vary the regularization parameter over two orders of magnitude around the value used in the main text and demonstrate that both asymmetry and bimodality persist. We further generate synthetic multi-harmonic datasets from known ground-truth distributions (bimodal and asymmetric) with added Gaussian noise matching our experimental uncertainty levels. Cross-validation shows that MEM recovers the input features with high fidelity; these results are now presented in new Supplementary Figures S5–S7 together with quantitative error metrics. The added material directly addresses the stability concern. revision: yes

Circularity Check

0 steps flagged

No significant circularity; reconstruction is data-driven inversion

full rationale

The paper's derivation chain starts from measured multi-harmonic intensities (SHG, THG, FHG) and applies the Maximum Entropy Method under the standard molecular susceptibility tensor model to obtain the orientation PDF. This is an external-data inversion, not a self-definitional loop, fitted parameter renamed as prediction, or load-bearing self-citation. MEM regularization is a standard external choice that does not force the output to equal an input by construction. The abstract and description give no evidence of ansatz smuggling, uniqueness theorems imported from the same authors, or renaming of known results. The central claim remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The method rests on the domain assumption that harmonic intensities are linear functionals of the orientation distribution via the molecular hyperpolarizability tensor, plus the standard MEM regularization choice.

free parameters (1)
  • MEM regularization parameter
    Maximum entropy method requires a regularization strength parameter whose value affects the smoothness of the recovered distribution.
axioms (1)
  • domain assumption Nonlinear polarization response is fully determined by the molecular orientation distribution through the susceptibility tensor
    Invoked when mapping measured harmonic signals to the distribution function.

pith-pipeline@v0.9.0 · 5512 in / 1300 out tokens · 48238 ms · 2026-05-10T07:00:41.254141+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    The many facets of molecular orientation in organic optoelectronics,

    1A. Hofmann, M. Schmid, and W. Br¨ utting, “The many facets of molecular orientation in organic optoelectronics,” Advanced Optical Materials9, 2101004 (2021). 2E. Pakhomenko, S. He, and R. J. Holmes, “Understanding and engineering spontaneous orientation polarization in organic light-emitting devices,” Chemical Physics Reviews4, 021308 (2023). 3D. Yokoyam...

    work page 2021

  2. [2]

    Spectroscopy of molecular monolayers by resonant second-harmonic generation,

    p. 229–255. 61T. F. Heinz, C. Chen, D. Ricard, and Y. Shen, “Spectroscopy of molecular monolayers by resonant second-harmonic generation,” Physical Review Letters48, 478 (1982). 62T. Heinz, H. Tom, and Y. R. Shen, “Determination of molecular orientation of monolayer adsorbates by optical second-harmonic generation,” Physical Review A28, 1883 (1983). 27 63...

    work page 1982