pith. sign in

arxiv: 2604.25425 · v1 · submitted 2026-04-28 · 🪐 quant-ph · cond-mat.supr-con· physics.optics

Elucidating mechanism of optical cavities in superconducting strip single photon detectors using transmission line and impedance models

Hiroki Kutsuma , Taro Yamashita This is my paper

Pith reviewed 2026-05-07 16:59 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.supr-conphysics.optics
keywords SSPDoptical cavitiesimpedance matchingabsorptancetransmission line modelsuperconducting detectors
0
0 comments X

The pith

SSPDs with optical cavities maximize absorptance when input impedance matches the impedance of the input medium.

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

The paper uses transmission line models to derive exact analytical formulas for light absorptance in superconducting strip single photon detectors equipped with single-side, double-side, or multilayer optical cavities. These formulas match results from numerical simulations to high accuracy. An impedance model then shows that maximum absorptance occurs exactly when the combined input impedance of the SSPD and cavity equals the impedance of the incident medium. This matching supplies a direct design rule for cavity geometry and material choices. The same principle is stated to extend to other superconducting detectors such as microwave kinetic inductance detectors and transition-edge sensors.

Core claim

SSPDs with optical cavities achieve the maximum absorptance when their input impedance matches the impedance of the input medium.

What carries the argument

The impedance model, which identifies maximum absorptance as the condition of impedance matching between the SSPD-cavity structure and the input medium.

If this is right

  • Analytical absorptance formulas enable rapid design of single-side, double-side, and dielectric multilayer cavities without repeated numerical simulations.
  • Impedance matching supplies an explicit criterion for choosing cavity thickness, dielectric constants, and strip geometry.
  • The same matching principle applies to the design of microwave kinetic inductance detectors and transition-edge sensors.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Impedance matching may function as a general optimization rule for cavity-enhanced superconducting sensors across different wavelengths.
  • Measuring input impedance on fabricated devices with intentionally varied cavity parameters would test whether the ideal matching condition survives real fabrication tolerances.
  • Adding small resistive or dispersive terms to the models could quantify how far real devices depart from the predicted maximum absorptance.

Load-bearing premise

The transmission line and impedance models capture all relevant optical and electrical behavior of the real SSPD-cavity system without significant unmodeled losses, dispersion, or fabrication imperfections.

What would settle it

Fabricate or simulate an SSPD-cavity device whose measured input impedance deviates from the input-medium impedance and show that absorptance is not at its predicted maximum.

Figures

Figures reproduced from arXiv: 2604.25425 by Hiroki Kutsuma, Taro Yamashita.

Figure 1
Figure 1. Figure 1: Schematics of SSPDs with optical cavities. (a) SSPD with single-side optical cavity. (b) SSPD with double-side optical cavity. (c) SSPD with dielectric multi-layer optical cavity. The two dielectric layers with different refractive indices are stacked in N periods. where γw, ηw, and dw are the propagation constant, characteristic impedance, and thickness of the wire layer, respectvely. Because we have set … view at source ↗
Figure 2
Figure 2. Figure 2: Comparisons of analytical and simulation results. In each figure, the results represented by the dotted lines are calculated using the analytical formulae given by the transmission line model, and those represented by dashed (solid) lines are obtained from the simulations using RCWA (FEM). (a) Dependence of the absorptance on the thickness of the wire layer of SSPDs with single-side optical cavities. (b) D… view at source ↗
Figure 3
Figure 3. Figure 3: Comparisons between analytical and simulation results. In each figure, the results represented by the dotted lines are calculated using the analytical formula given by the transmission line model, and those represented by the dashed (solid) lines are obtained from simulations using RCWA (FEM). (a) Dependence of the absorptance on the thickness of the dielectric layer of SSPDs with single-side optical cavit… view at source ↗
Figure 4
Figure 4. Figure 4: Dependence of the absorptance and ratio of the input impedance to the impedance of the input media on the thickness of the wire layer. The upper panel of each figure shows the absorptance as a function of the thickness of the wire layer, determined by Eq. (9) in (a), Eq. (17) in (b), and Eq. (34) in (c). The lower panel shows the ratio of the absolute value of the input impedance of SSPDs with the cavities… view at source ↗
read the original abstract

We clarified the physical mechanism of superconducting strip single photon detectors (SSPDs) with optical cavities by using transmission line and impedance models. By introducing the transmission line model, we derived the analytical formulae for the absorptance of SSPDs with optical cavities. We compared the absorptance obtained from the analytical formulae for SSPDs with single-side, double-side, and dielectric multi-layer optical cavities against the results of numerical simulations. The comparison showed that the results were nearly identical. By introducing the impedance model, it was clearly shown that the SSPDs with optical cavities achieved the maximum absorptance when their input impedance of the SSPDs with optical cavities matched the impedance of the input medium. The design concepts proposed in this study are applicable to other superconducting detectors, such as microwave kinetic inductance detectors and transition-edge sensors.

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

1 major / 2 minor

Summary. The manuscript derives closed-form analytical expressions for the absorptance of superconducting strip single-photon detectors (SSPDs) that incorporate optical cavities, using a transmission-line (transfer-matrix) representation of the optical stack. It reports near-exact numerical agreement between these formulae and independent simulations for single-side, double-side, and dielectric multilayer cavity geometries. The work further shows, via an impedance model, that peak absorptance occurs precisely when the computed input impedance of the SSPD-cavity system equals the impedance of the incident medium; the design insight is stated to generalize to other superconducting detectors such as MKIDs and TESs.

Significance. If the derivations and comparisons hold, the paper supplies a transparent, model-based explanation—rooted in standard electromagnetic impedance matching—for why optical cavities enhance absorption in SSPDs. The explicit provision of analytical formulae together with direct side-by-side validation against numerical simulations is a clear strength, offering a route to design optimization that does not rely solely on black-box simulation. This contributes constructively to the SSPD literature by making the underlying mechanism explicit rather than phenomenological.

major comments (1)
  1. The central claim that the analytical formulae agree 'nearly identically' with simulations for three distinct cavity types is load-bearing for the paper's validation strategy, yet the manuscript provides neither an explicit error metric (e.g., maximum relative deviation or RMS error) nor the precise layer thicknesses, refractive indices, and strip parameters used in the comparisons. Without these, the robustness of the agreement cannot be independently assessed.
minor comments (2)
  1. A compact table listing peak absorptance values (analytical vs. simulated) together with relative errors for each cavity configuration would make the 'nearly identical' statement quantitative and easier to evaluate.
  2. The impedance-matching condition is derived directly from the reflection coefficient in the same transmission-line model; a brief remark clarifying that this is a mathematical identity within the model (rather than an independent physical discovery) would prevent any misreading.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of our work and the recommendation for minor revision. We address the major comment below.

read point-by-point responses

  1. Referee: The central claim that the analytical formulae agree 'nearly identically' with simulations for three distinct cavity types is load-bearing for the paper's validation strategy, yet the manuscript provides neither an explicit error metric (e.g., maximum relative deviation or RMS error) nor the precise layer thicknesses, refractive indices, and strip parameters used in the comparisons. Without these, the robustness of the agreement cannot be independently assessed.

    Authors: We agree that explicit quantitative error metrics and the precise parameters are necessary for independent verification. In the revised manuscript we will supply the layer thicknesses, refractive indices, and strip parameters for each of the three cavity geometries. We will also report the maximum relative deviation and RMS error between the analytical formulae and the numerical simulations for every case, either in the main text or in a supplementary note. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivations rest on standard EM models

full rationale

The paper derives closed-form absorptance expressions via the standard transmission-line (transfer-matrix) representation of the optical stack and verifies them against independent numerical simulations for single-side, double-side, and multilayer cavities, showing near-exact agreement. The statement that maximum absorptance occurs when input impedance matches the incident-medium impedance follows directly from the reflection coefficient formula inherent to any transmission-line model, which is a general property of electromagnetism rather than a self-defined or fitted result. No load-bearing self-citations, uniqueness theorems from prior author work, or ansatzes smuggled via citation are present; the central claims remain self-contained against external benchmarks and do not reduce to their inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard electromagnetic transmission-line theory and circuit impedance concepts applied to the detector geometry; no new physical entities are postulated.

axioms (2)
  • domain assumption The SSPD-cavity structure can be accurately represented as a transmission line for optical frequencies.
    Invoked to derive analytical absorptance formulas.
  • standard math Maximum power transfer occurs when input impedance matches the source impedance.
    Used to explain peak absorptance condition.

pith-pipeline@v0.9.0 · 5446 in / 1269 out tokens · 42636 ms · 2026-05-07T16:59:27.350084+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

47 extracted references · 45 canonical work pages

  1. [1]

    N.et al.Picosecond superconducting single-photon optical detector.Appl

    Gol’tsman, G. N.et al.Picosecond superconducting single-photon optical detector.Appl. Phys. Lett.79, 705–707, DOI: 10.1063/1.1388868 (2001)

    work page doi:10.1063/1.1388868 2001

  2. [2]

    Tversky and D

    Zhong, H.-S.et al.Quantum computational advantage using photons.Science370, 1460–1463, DOI: 10.1126/science. abe8770 (2020)

    work page doi:10.1126/science 2020

  3. [3]

    Field test of quantum key distribution in the tokyo qkd network,

    Sasaki, M.et al.Field test of quantum key distribution in the Tokyo QKD Network.Opt. Express19, 10387–10409, DOI: 10.1364/OE.19.010387 (2011)

    work page doi:10.1364/oe.19.010387 2011

  4. [4]

    Photonics10, 441–444, DOI: https: //doi.org/10.1038/nphoton.2016.74 (2016)

    Kobayashi, T.et al.Frequency-domain Hong–Ou–Mandel interference.Nat. Photonics10, 441–444, DOI: https: //doi.org/10.1038/nphoton.2016.74 (2016). 12/15

    work page doi:10.1038/nphoton.2016.74 2016

  5. [5]

    Express29, 11728–11738, DOI: 10.1364/OE.423142 (2021)

    Endo, M.et al.Quantum detector tomography of a superconducting nanostrip photon-number-resolving detector.Opt. Express29, 11728–11738, DOI: 10.1364/OE.423142 (2021)

    work page doi:10.1364/oe.423142 2021

  6. [6]

    Photonics7, 210–214, DOI: 10.1038/nphoton.2013.13 (2013)

    Marsili, F.et al.Detecting single infrared photons with 93% system efficiency.Nat. Photonics7, 210–214, DOI: 10.1038/nphoton.2013.13 (2013)

    work page doi:10.1038/nphoton.2013.13 2013

  7. [7]

    China Physics, Mech

    Zhang, W.et al.NbN superconducting nanowire single photon detector with efficiency over 90% at 1550 nm wavelength operational at compact cryocooler temperature.Sci. China Physics, Mech. & Astron.60, 1–10, DOI: 10.1007/s11433-017-9113-4 (2017)

    work page doi:10.1007/s11433-017-9113-4 2017

  8. [8]

    V ., Nerem, R

    Reddy, D. V ., Nerem, R. R., Nam, S. W., Mirin, R. P. & Verma, V . B. Superconducting nanowire single-photon detectors with 98% system detection efficiency at 1550 nm.Optica7, 1649–1653, DOI: /10.1364/OPTICA.400751 (2020)

    work page doi:10.1364/optica.400751 2020

  9. [9]

    Express28, 36884–36891, DOI: 10.1364/OE.410025 (2020)

    Hu, P.et al.Detecting single infrared photons toward optimal system detection efficiency.Opt. Express28, 36884–36891, DOI: 10.1364/OE.410025 (2020)

    work page doi:10.1364/oe.410025 2020

  10. [10]

    APL Photonics6, 036114, DOI: 10.1063/5.0039772 (2021)

    Chang, J.et al.Detecting telecom single photons with (99.5+0.5 −2.07)% system detection efficiency and high time resolution. APL Photonics6, 036114, DOI: 10.1063/5.0039772 (2021)

    work page doi:10.1063/5.0039772 2021

  11. [11]

    Photonics Res.9, 958–967, DOI: 10.1364/PRJ.419514 (2021)

    Xu, G.-Z.et al.Superconducting microstrip single-photon detector with system detection efficiency over 90% at 1550 nm. Photonics Res.9, 958–967, DOI: 10.1364/PRJ.419514 (2021)

    work page doi:10.1364/prj.419514 2021

  12. [12]

    V ., Otrooshi, N., Nam, S

    Reddy, D. V ., Otrooshi, N., Nam, S. W., Mirin, R. P. & Verma, V . B. Broadband polarization insensitivity and high detection efficiency in high-fill-factor superconducting microwire single-photon detectors.APL Photonics7, 051302, DOI: 10.1063/5.0088007 (2022)

    work page doi:10.1063/5.0088007 2022

  13. [13]

    Miki, S.et al.Superconducting coincidence photon detector with short timing jitter.Appl. Phys. Lett.112, 262601, DOI: 10.1063/1.5037254 (2018)

    work page doi:10.1063/1.5037254 2018

  14. [14]

    Korzh, B.et al.Demonstration of sub-3 ps temporal resolution with a superconducting nanowire single-photon detector. Nat. Photonics14, 250–255, DOI: 10.1038/s41566-020-0589-x (2020)

    work page doi:10.1038/s41566-020-0589-x 2020

  15. [15]

    Colangelo, M.et al.Impedance-matched differential superconducting nanowire detectors.Phys. Rev. Appl.19, 044093, DOI: 10.1103/PhysRevApplied.19.044093 (2023)

    work page doi:10.1103/physrevapplied.19.044093 2023

  16. [16]

    Express22, 16267– 16272, DOI: https://doi.org/10.1364/OE.22.016267 (2014)

    Yang, X.et al.Superconducting nanowire single photon detector with on-chip bandpass filter.Opt. Express22, 16267– 16272, DOI: https://doi.org/10.1364/OE.22.016267 (2014)

    work page doi:10.1364/oe.22.016267 2014

  17. [17]

    & Tokura, Y

    Shibata, H., Shimizu, K., Takesue, H. & Tokura, Y . Ultimate low system dark-count rate for superconducting nanowire single-photon detector.Opt. Lett.40, 3428–3431, DOI: 10.1364/OL.40.003428 (2015)

    work page doi:10.1364/ol.40.003428 2015

  18. [18]

    Zhang, W.et al.Fiber-coupled superconducting nanowire single-photon detectors integrated with a bandpass filter on the fiber end-face.Supercond. Sci. Technol.31, 035012, DOI: 10.1088/1361-6668/aaa6b4 (2018)

    work page doi:10.1088/1361-6668/aaa6b4 2018

  19. [19]

    Chiles, J.et al.New constraints on dark photon dark matter with superconducting nanowire detectors in an optical haloscope.Phys. Rev. Lett.128, 231802, DOI: 10.1103/PhysRevLett.128.231802 (2022)

    work page doi:10.1103/physrevlett.128.231802 2022

  20. [20]

    Zhang, W.et al.A 16-pixel interleaved superconducting nanowire single-photon detector array with a maximum count rate exceeding 1.5 GHz.IEEE Trans. Appl. Supercond.29, 1–4, DOI: 10.1109/TASC.2019.2895621 (2019)

    work page doi:10.1109/tasc.2019.2895621 2019

  21. [21]

    V .et al.Gigahertz detection rates and dynamic photon-number resolution with superconducting nanowire arrays

    Resta, G. V .et al.Gigahertz detection rates and dynamic photon-number resolution with superconducting nanowire arrays. Nano Lett.23, 6018–6026, DOI: 10.1021/acs.nanolett.3c01228. (2023)

    work page doi:10.1021/acs.nanolett.3c01228 2023

  22. [22]

    Craiciu, I.et al.High-speed detection of 1550 nm single photons with superconducting nanowire detectors.Optica10, 183–190, DOI: 10.1364/OPTICA.478960 (2023)

    work page doi:10.1364/optica.478960 2023

  23. [23]

    Hao, H.et al.A compact multi-pixel superconducting nanowire single-photon detector array supporting gigabit space-to- ground communications.Light. Sci. Appl.13, 25, DOI: 10.1038/s41377-023-01374-1 (2024)

    work page doi:10.1038/s41377-023-01374-1 2024

  24. [24]

    Verevkin, A.et al.Detection efficiency of large-active-area NbN single-photon superconducting detectors in the ultraviolet to near-infrared range.Appl. Phys. Lett.80, 4687–4689, DOI: 10.1063/1.1487924 (2002)

    work page doi:10.1063/1.1487924 2002

  25. [25]

    M.et al.Nanowire single-photon detector with an integrated optical cavity and anti-reflection coating.Opt

    Rosfjord, K. M.et al.Nanowire single-photon detector with an integrated optical cavity and anti-reflection coating.Opt. Express14, 527–534, DOI: 10.1364/OPEX.14.000527 (2006)

    work page doi:10.1364/opex.14.000527 2006

  26. [26]

    J.et al.Constriction-limited detection efficiency of superconducting nanowire single-photon detectors.Appl

    Kerman, A. J.et al.Constriction-limited detection efficiency of superconducting nanowire single-photon detectors.Appl. Phys. Lett.90, 101110, DOI: 10.1063/1.2696926 (2007)

    work page doi:10.1063/1.2696926 2007

  27. [27]

    & Wang, Z

    Miki, S., Takeda, M., Fujiwara, M., Sasaki, M. & Wang, Z. Compactly packaged superconducting nanowire single-photon detector with an optical cavity for multichannel system.Opt. Express17, 23557–23564, DOI: 10.1364/OE.17.023557 (2009). 13/15

    work page doi:10.1364/oe.17.023557 2009

  28. [28]

    & Wang, Z

    Miki, S., Yamashita, T., Fujiwara, M., Sasakiand, M. & Wang, Z. Multichannel SNSPD system with high detection efficiency at telecommunication wavelength.Opt. Lett.35, 2133–2135, DOI: https://doi.org/10.1364/OL.35.002133 (2010)

    work page doi:10.1364/ol.35.002133 2010

  29. [29]

    Gaggero, A.et al.Nanowire superconducting single-photon detectors on GaAs for integrated quantum photonic applications. Appl. Phys. Lett.97, 151108, DOI: 10.1063/1.3496457 (2010)

    work page doi:10.1063/1.3496457 2010

  30. [30]

    J.et al.Compact cryogenic self-aligning fiber-to-detector coupling with losses below one percent.Opt

    Miller, A. J.et al.Compact cryogenic self-aligning fiber-to-detector coupling with losses below one percent.Opt. Express 19, 9102–9110, DOI: 10.1364/OE.19.009102 (2011)

    work page doi:10.1364/oe.19.009102 2011

  31. [31]

    Clem, J. R. & Berggren, K. K. Geometry-dependent critical currents in superconducting nanocircuits.Phys. Rev. B84, 174510, DOI: 10.1103/PhysRevB.84.174510 (2011)

    work page doi:10.1103/physrevb.84.174510 2011

  32. [32]

    & Wang, Z

    Miki, S., Yamashita, T., Terai, H. & Wang, Z. High performance fiber-coupled NbTiN superconducting nanowire single photon detectors with Gifford-McMahon cryocooler.Opt. Express21, 10208–10214, DOI: 10.1364/OE.21.010208 (2013)

    work page doi:10.1364/oe.21.010208 2013

  33. [33]

    Fang, Y .-Q.et al.InGaAs/InP single-photon detectors with 60% detection efficiency at 1550 nm.Rev. Sci. Instrum.91, DOI: 10.1063/5.0014123 (2020)

    work page doi:10.1063/5.0014123 2020

  34. [34]

    Lett.33, 639–641, DOI: 10.1364/OL.33.000639 (2008)

    Kamada, H.et al.Efficient and low-noise single-photon detection in 1550 nm communication band by frequency upconversion in periodically poled LiNbO3 waveguides.Opt. Lett.33, 639–641, DOI: 10.1364/OL.33.000639 (2008)

    work page doi:10.1364/ol.33.000639 2008

  35. [35]

    Express16, 10750–10761, DOI: 10.1364/OE.16.010750 (2008)

    Anant, V .et al.Optical properties of superconducting nanowire single-photon detectors.Opt. Express16, 10750–10761, DOI: 10.1364/OE.16.010750 (2008)

    work page doi:10.1364/oe.16.010750 2008

  36. [36]

    Yamashita , author S

    Yamashita, T., Miki, S., Terai, H. & Wang, Z. Low-filling-factor superconducting single photon detector with high system detection efficiency.Opt. Express21, 27177–27184, DOI: 10.1364/OE.21.027177 (2013)

    work page doi:10.1364/oe.21.027177 2013

  37. [37]

    Moharam, M. G. & Gaylord, T. K. Rigorous coupled-wave analysis of planar-grating diffraction.JOSA71, 811–818, DOI: 10.1364/JOSA.71.000811 (1981)

    work page doi:10.1364/josa.71.000811 1981

  38. [38]

    J., Clough, R

    Turner, M. J., Clough, R. W., Martin, H. C. & Topp, L. J. Stiffness and deflection analysis of complex structures.J. Aeronaut. Sci.23, 805–823 (1956)

    1956

  39. [39]

    Driessen, E. F. C.et al.Impedance model for the polarization-dependent optical absorption of superconducting single- photon detectors.Eur . Phys. J. - Appl. Phys.47, 10701, DOI: 10.1051/epjap/2009087 (2009)

    work page doi:10.1051/epjap/2009087 2009

  40. [40]

    Low Temp

    Kouwenhoven, K.et al.Model and measurements of an optical stack for broadband visible to near-infrared absorption in TiN MKIDs.J. Low Temp. Phys.209, 1249–1257, DOI: 10.1007/s10909-022-02774-0 (2022). 41.Pozar, D. M.Microwave engineering: theory and techniques, 4th Edition(John Wiley & Sons, 2011)

    work page doi:10.1007/s10909-022-02774-0 2022

  41. [41]

    & Miki, S

    Yabuno, M., China, F., Terai, H. & Miki, S. Superconducting wide strip photon detector with high critical current bank structure.Opt. Quantum1, 26–34, DOI: 10.1364/OPTICAQ.497675 (2023)

    work page doi:10.1364/opticaq.497675 2023

  42. [42]

    Conversions between S, Z, Y , H, ABCD, and T parameters which are valid for complex source and load impedances.IEEE Trans

    Frickey, D. Conversions between S, Z, Y , H, ABCD, and T parameters which are valid for complex source and load impedances.IEEE Trans. Microw. Theory Tech.42, 205–211, DOI: 10.1109/22.275248 (1994)

    work page doi:10.1109/22.275248 1994

  43. [43]

    Aspnes, D. E. Local-field effects and effective-medium theory: A microscopic perspective.Am. J. Phys50, 704–709 (1982)

    1982

  44. [44]

    Computer Physics Communications , author=

    Liu, V . & Fan, S. S4: A free electromagnetic solver for layered periodic structures.Comput. Phys. Commun.183, 2233–2244, DOI: 10.1016/j.cpc.2012.04.026 (2012)

    work page doi:10.1016/j.cpc.2012.04.026 2012

  45. [45]

    & Yamashita, T

    Kutsuma, H. & Yamashita, T. Source codes for ’elucidating mechanism of optical cavities in superconducting strip single photon detectors using transmission line and impedance models’, DOI: https://doi.org/10.5281/zenodo.19722815 (2025)

    work page doi:10.5281/zenodo.19722815 2025

  46. [46]

    Nature425, 817–821 (2003) https://doi.org/10.1038/nature02037

    Day, P. K., LeDuc, H. G., Mazin, B. A., Vayonakis, A. & Zmuidzinas, J. A broadband superconducting detector suitable for use in large arrays.Nature425, 817–821, DOI: 10.1038/nature02037 (2003)

    work page doi:10.1038/nature02037 2003

  47. [47]

    Irwin, K. D. An application of electrothermal feedback for high resolution cryogenic particle detection.Appl. Phys. Lett. 66, 1998–2000, DOI: 10.1063/1.113674 (1995). Acknowledgements Hiroki Kutsuma and Taro Yamashita acknowledge the Center for Heterogeneous Quantum/Material Fusion Technologies and Center for Key Interdisciplinary Research at Tohoku Unive...

    work page doi:10.1063/1.113674 1998