pith. sign in

arxiv: 2604.04856 · v1 · submitted 2026-04-06 · 🪐 quant-ph · cond-mat.mes-hall· cond-mat.stat-mech· physics.optics

Modeling the non-Markovian Brownian motion of an optomechanical resonator

Aritra Ghosh , Malay Bandyopadhyay , M. Bhattacharya This is my paper

Pith reviewed 2026-05-10 18:57 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.mes-hallcond-mat.stat-mechphysics.optics
keywords non-Markovian Brownian motionoptomechanical resonatorbath spectral densitymechanical susceptibilitynon-Ohmic spectrummemory effectshomodyne detection
0
0 comments X p. Extension

The pith

A phenomenological spectral density for the bath reproduces local non-Ohmic behavior near resonance while remaining globally admissible.

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

The paper develops a model for the environment of an optomechanical resonator that matches the non-Ohmic spectrum observed experimentally near the mechanical resonance frequency. It avoids divergences that would appear if the local power-law form were used at all frequencies, keeping frequency renormalizations and quadrature fluctuations finite. The resulting nonlocal mechanical susceptibility has an analytic pole structure that directly encodes the measured linewidth. The dissipation kernel then shows a power-law-modulated exponential decay that includes intervals of negativity, indicating strong memory effects in the resonator motion. In the weak-coupling regime this construction also shows how homodyne optical readout, combined with a calibrated drive, can in principle recover the full susceptibility and thereby separate dissipative and dispersive contributions from the bath.

Core claim

The authors construct a globally admissible phenomenological spectral density of the bath that reproduces the experimentally observed local power-law spectrum near the mechanical resonance while remaining well defined at all frequencies. This yields a nonlocal mechanical susceptibility whose analytic pole structure encodes the observed linewidth. The associated dissipation kernel exhibits power-law-modulated exponential decay with transient negativity, signaling strong non-Markovian memory effects. In the weak-coupling limit, homodyne detection permits near-resonance spectroscopy and, with a calibrated drive, reconstruction of the full mechanical susceptibility.

What carries the argument

The globally admissible phenomenological bath spectral density, built to match the local power-law observations near resonance while ensuring finite bath-induced renormalizations everywhere.

If this is right

  • The mechanical susceptibility becomes nonlocal, with its analytic poles fixing the resonator linewidth.
  • The dissipation kernel displays power-law-modulated exponential decay containing transient negative regions.
  • Homodyne optical readout enables near-resonance spectroscopy of the resonator.
  • A calibrated drive on the resonator permits reconstruction of the complete mechanical susceptibility, separating dissipative and dispersive bath effects.

Where Pith is reading between the lines

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

  • The same construction could be tested in other open quantum systems where only local spectral data are available, to check whether it systematically prevents unphysical divergences.
  • Direct time-domain measurements of the resonator response might reveal the predicted transient negativity and thereby confirm the strength of memory effects.
  • The method supplies a template for embedding any locally measured bath spectrum into a globally consistent open-system description.

Load-bearing premise

A single functional form can be chosen that exactly reproduces the local power-law spectrum near resonance yet stays free of divergences in all bath-induced renormalizations and fluctuations.

What would settle it

An experimental reconstruction of the mechanical susceptibility that fails to exhibit poles whose imaginary parts match the observed linewidth, or a measured divergence in frequency shift when the local spectral form is used globally.

Figures

Figures reproduced from arXiv: 2604.04856 by Aritra Ghosh, Malay Bandyopadhyay, M. Bhattacharya.

Figure 1
Figure 1. Figure 1: FIG. 1: Schematic of an optomechanical setup where a [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Normalized bath spectral function [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Normalized dissipation kernel [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
read the original abstract

We propose a globally-admissible phenomenological spectral density of the bath for the non-Markovian Brownian motion of an optomechanical resonator, motivated by the near-resonance experimental observation of a non-Ohmic spectrum in [Nat. Commun. 6, 7606 (2015)]. To avoid divergences arising from a naive global extrapolation, we construct this phenomenological bath spectral density that reproduces the observed local-power-law behavior near the mechanical resonance while remaining well defined globally, ensuring the finiteness of the bath-induced renormalizations and quadrature fluctuations of the resonator. The corresponding model of the structured environment produces a nonlocal mechanical susceptibility whose analytic pole structure encodes the observed linewidth. The resulting dissipation kernel exhibits a power-law-modulated exponential decay with transient negativity, signaling strong memory effects. In the weak-coupling regime, the optical readout based on homodyne detection enables near-resonance spectroscopy and, with a calibrated drive on the resonator, permits, in principle, the reconstruction of the full mechanical susceptibility, thereby providing access to both the dissipative and dispersive bath contributions. Our results provide a consistent route from locally-inferred spectral properties to globally-admissible open-system descriptions and establish a framework for probing structured environments in cavity optomechanics.

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 manuscript proposes a globally admissible phenomenological spectral density J(ω) for the bath in the non-Markovian Brownian motion of an optomechanical resonator, motivated by the near-resonance non-Ohmic spectrum observed in prior experiments. The construction is designed to reproduce the local power-law behavior near the mechanical resonance frequency while avoiding divergences in bath-induced renormalizations and quadrature fluctuations. This yields a nonlocal mechanical susceptibility whose analytic pole structure encodes the observed linewidth, a dissipation kernel exhibiting power-law-modulated exponential decay with transient negativity, and proposals for homodyne-based optical readout to reconstruct the full susceptibility in the weak-coupling regime.

Significance. If the explicit functional form and verifications are provided, the work offers a consistent route from locally observed spectral properties to globally admissible open-system models in cavity optomechanics. It highlights strong memory effects via the kernel's negativity and enables access to both dissipative and dispersive bath contributions, strengthening the framework for probing structured environments.

major comments (2)
  1. [Abstract] The central construction of the phenomenological spectral density is asserted to exactly reproduce the local power-law while remaining globally admissible with all principal-value integrals finite, but the abstract provides no explicit functional form, derivation, or numerical demonstration of the match and convergence. This is load-bearing for the claim that the pole structure encodes the linewidth independently of the fitting choice.
  2. [Model of the structured environment] The spectral density is explicitly constructed to match the experimentally observed local power-law, so the resulting susceptibility poles encoding the linewidth appear tied to parameters chosen for that match, creating circularity between input data and output. The manuscript should clarify whether the poles arise as an independent prediction or by design, with explicit calculation of the pole locations from the model.
minor comments (2)
  1. [Introduction] The reference to the motivating experiment [Nat. Commun. 6, 7606 (2015)] is appropriate, but the manuscript could add a brief summary of how the local power-law exponent and cutoff were extracted from the data to aid reproducibility.
  2. [Dissipation kernel] The transient negativity in the dissipation kernel is a key signature of memory effects; a figure plotting the kernel versus time for representative parameters would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for the constructive comments. We have revised the manuscript to address the points raised, providing additional details on the spectral density construction and pole calculations. Our responses to the major comments are as follows.

read point-by-point responses

  1. Referee: [Abstract] The central construction of the phenomenological spectral density is asserted to exactly reproduce the local power-law while remaining globally admissible with all principal-value integrals finite, but the abstract provides no explicit functional form, derivation, or numerical demonstration of the match and convergence. This is load-bearing for the claim that the pole structure encodes the linewidth independently of the fitting choice.

    Authors: We agree that the abstract does not contain the explicit functional form due to length constraints. In the revised manuscript, we have updated the abstract to include a brief statement of the functional form used for J(ω). The derivation of this form, which ensures global admissibility by cutting off high-frequency divergences while matching the local non-Ohmic power-law behavior near the mechanical resonance, along with numerical demonstrations of the match and convergence of principal-value integrals, is provided in detail in the section 'Model of the structured environment'. We also show that the analytic pole structure of the susceptibility is determined by the parameters of this global model and encodes the linewidth in a manner that is robust to small variations in the fitting parameters, as verified through explicit calculations. revision: yes

  2. Referee: [Model of the structured environment] The spectral density is explicitly constructed to match the experimentally observed local power-law, so the resulting susceptibility poles encoding the linewidth appear tied to parameters chosen for that match, creating circularity between input data and output. The manuscript should clarify whether the poles arise as an independent prediction or by design, with explicit calculation of the pole locations from the model.

    Authors: The model is indeed phenomenological by design: the local power-law is matched to experimental observations, and the global spectral density is constructed to reproduce this behavior near resonance while ensuring the model is well-defined globally. The poles of the susceptibility are therefore a consequence of this construction rather than an independent prediction. In the revised manuscript, we have added an explicit calculation of the pole locations by solving the equation for the zeros of the inverse susceptibility, demonstrating that they correspond to the observed linewidth. This approach provides a consistent framework for extending local observations to a global open-system model, without circularity, as the global properties (such as the full susceptibility) are then derived from the model. revision: partial

Circularity Check

1 steps flagged

Phenomenological J(ω) constructed to match observed local power-law makes susceptibility poles encode the same observation by design

specific steps
  1. fitted input called prediction [Abstract]
    "we construct this phenomenological bath spectral density that reproduces the observed local-power-law behavior near the mechanical resonance while remaining well defined globally, ensuring the finiteness of the bath-induced renormalizations and quadrature fluctuations of the resonator. The corresponding model of the structured environment produces a nonlocal mechanical susceptibility whose analytic pole structure encodes the observed linewidth."

    J(ω) is defined by construction to match the local experimental power-law observation. The pole structure of the derived susceptibility (which sets the mechanical linewidth via the damping kernel) is therefore fixed by the same parameters, so the claim that the poles 'encode the observed linewidth' is a direct consequence of the input-matching step rather than a separate prediction.

full rationale

The derivation begins with an explicit construction of the bath spectral density to reproduce the experimentally observed local power-law near resonance (while ensuring global admissibility). The mechanical susceptibility is then obtained from this J(ω), and its analytic poles are presented as encoding the observed linewidth. Because the linewidth is determined by the damping kernel whose parameters are chosen to match the input spectrum, the encoding step reduces to the fitted construction rather than an independent derivation. This is a moderate instance of fitted_input_called_prediction; the remainder of the framework (global finiteness, kernel decay, readout protocol) adds independent content and does not collapse to the same reduction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central construction rests on choosing a functional form that interpolates the local experimental spectrum to all frequencies without introducing divergences; this choice is not derived from first principles but postulated to satisfy global regularity.

free parameters (1)
  • parameters defining the global extrapolation of the spectral density
    The functional form is chosen phenomenologically to match local power-law data while ensuring global finiteness; specific values are not stated but must be calibrated to the 2015 experiment.
axioms (1)
  • domain assumption The bath spectral density can be represented by a single phenomenological function that is locally non-Ohmic near the mechanical resonance and globally regular.
    Invoked to justify the construction that avoids divergences while reproducing the observed local behavior.

pith-pipeline@v0.9.0 · 5527 in / 1440 out tokens · 39800 ms · 2026-05-10T18:57:59.335557+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    At the observed mechanical resonance Ω R, Jk(ω)|ω≈ΩR ∝ω k, withk=−2.30±1.05 [10]

    work page

  2. [2]

    Λθ(ω) ˜Fext(ω) ˜X(mech,coh) θ,out (ω) # , (60) and Jk(ω) = Im[Σ(ω)] =−Im

    BothδK[Eq. (7)] andδM[Eq. (11)] should be finite, i.e., non-divergent. A convenient common and sufficient condition ensuring the ultraviolet and infrared convergence of both the renormalization integrals is that the spectral density should satisfy Jk(ω)∼ O(ω s) withs >2 asω→0, andJ k(ω)∼ O(ωr) withr <0 asω→ ∞. These sufficient conditions motivate the foll...

    work page

  3. [3]

    By choosinga=b= 3−k, one obtains the exact expressions δK= AkΩ3 R 2√π Γ 3 2 −k Γ(3−k) , δM= AkΩR√π Γ 5 2 −k Γ(3−k) ,(C4) which, upon substituting Eq. (20) gives δK= 22−kJk(ΩR)√π Γ 3 2 −k Γ(3−k) , δM= 23−kJk(ΩR) Ω2 R √π Γ 5 2 −k Γ(3−k) .(C5) The bath-induced stiffness and mass renormalizations of the resonator are therefore determined by the on-resonance v...

    work page

  4. [4]

    Weiss,Dissipative quantum systems, 5th ed., World Scientific (2021)

    U. Weiss,Dissipative quantum systems, 5th ed., World Scientific (2021)

    work page 2021

  5. [5]

    Breuer and F

    H.-P. Breuer and F. Petruccione,The theory of open quantum systems, Oxford University Press (2002)

    work page 2002

  6. [6]

    Seoanez, F

    C. Seoanez, F. Guinea, and A. H. Castro Neto,Dissipa- tion due to two-level systems in nano-mechanical devices, EPL78, 60002 (2007)

    work page 2007

  7. [7]

    Wilson-Rae,Intrinsic dissipation in nanomechanical resonators due to phonon tunneling, Phys

    I. Wilson-Rae,Intrinsic dissipation in nanomechanical resonators due to phonon tunneling, Phys. Rev. B77, 245418 (2008)

    work page 2008

  8. [8]

    Grabert, P

    H. Grabert, P. Schramm, and G.-L. Ingold,Quan- tum Brownian motion: The functional integral approach, Phys. Rep.168, 115 (1988)

    work page 1988

  9. [9]

    H¨ anggi and G.-L

    P. H¨ anggi and G.-L. Ingold,Fundamental aspects of quantum Brownian motion, Chaos15, 026105 (2005)

    work page 2005

  10. [10]

    Ghosh, M

    A. Ghosh, M. Bandyopadhyay, S. Dattagupta, and S. Gupta,Independent-oscillator model and the quantum Langevin equation for an oscillator: a review, J. Stat. Mech.: Theory Exp.2024, 074002 (2024)

    work page 2024

  11. [11]

    T. J. Kippenberg and K. J. Vahala,Cavity opto- mechanics, Opt. Exp.15, 17172 (2007)

    work page 2007

  12. [12]

    Aspelmeyer, T

    M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, Cavity optomechanics, Rev. Mod. Phys.86, 1391 (2014)

    work page 2014

  13. [13]

    Gr¨ oblacher, A

    S. Gr¨ oblacher, A. Trubarov, N. Prigge, G. D. Cole, M. Aspelmeyer, and J. Eisert,Observation of non- Markovian micromechanical Brownian motion, Nat. Commun.6, 7606 (2015)

    work page 2015

  14. [14]

    Jiang, R

    W. Jiang, R. Yang, T.-H. Qiu, and G.-J. Yang,Probe response of a cavity-optomechanical system coupling to a frequency-dependent bath, Phys. Rev. A101, 033804 (2020)

    work page 2020

  15. [15]

    Zhang, X.-T

    W.-Z. Zhang, X.-T. Liang, J. Cheng, and L. Zhou,Mea- surement of the mechanical reservoir spectral density in an optomechanical system, Phys. Rev. A103, 053707 (2021)

    work page 2021

  16. [16]

    Ullersma,An exactly solvable model for Brownian mo- tion: I

    P. Ullersma,An exactly solvable model for Brownian mo- tion: I. Derivation of the Langevin equation, Physica32, 27 (1966)

    work page 1966

  17. [17]

    G. W. Ford, J. T. Lewis, and R. F. O’Connell,Quantum Langevin equation, Phys. Rev. A37, 4419 (1988)

    work page 1988

  18. [18]

    R. P. Feynman and F. L. Vernon Jr.,The theory of a gen- eral quantum system interacting with a linear dissipative system, Ann. Phys. (N.Y.)24, 118 (1963)

    work page 1963

  19. [19]

    G. W. Ford, M. Kac, and P. Mazur,Statistical mechanics of assemblies of coupled oscillators, J. Math. Phys.6, 504 (1965)

    work page 1965

  20. [20]

    A. O. Caldeira and A. J. Leggett,Path integral ap- proach to quantum Brownian motion, Physica A121, 587 (1983)

    work page 1983

  21. [21]

    Grabert, U

    H. Grabert, U. Weiss, and P. Talkner,Quantum theory of the damped harmonic oscillator, Z. Phys. B55, 87 (1984)

    work page 1984

  22. [22]

    Karrlein and H

    R. Karrlein and H. Grabert,Exact time evolution and master equations for the damped harmonic oscillator, Phys. Rev. E55, 153 (1997)

    work page 1997

  23. [23]

    H. M. Wiseman and G. J. Milburn,Quantum theory of optical feedback via homodyne detection, Phys. Rev. Lett. 70, 548 (1993)

    work page 1993

  24. [24]

    A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and R. J. Schoelkopf,Introduction to quantum noise, measurement, and amplification, Rev. Mod. Phys.82, 1155 (2010)

    work page 2010

  25. [25]

    C. W. Gardiner and M. J. Collett,Input and output in damped quantum systems: Quantum stochastic differen- tial equations and the master equation, Phys. Rev. A31, 3761 (1985)

    work page 1985

  26. [26]

    I. S. Gradshteyn and I. M. Ryzhik,Table of integrals, se- ries, and products, edited by A. Jeffrey and D. Zwillinger, 7th ed., Academic Press (2007)

    work page 2007