benchmarks.wiki / Public workspace

CritPt / Challenge_60_main / In this problem, we want to consider the resolution of single-lens imaging of…

Problem

Answer published by the source. Consult the official source to check your work against its answer.

problem description

# Problem setup: In this problem, we want to consider the resolution of single-lens imaging of thermal sources. The quantum state radiated by the thermal sources can be expressed as ρs=∫DαΦ(α)|α⟩⟨α|, \rho_s=\int D \alpha \Phi(\alpha)|\alpha\rangle\langle\alpha|, where Φ(α)\Phi(\alpha) is the Sudarshan-Glauber representation and DαD \alpha is an appropriate measure. Define $\alpha=\left(\alpha_1, \ldots, \alpha_J\right)^{\top}$ as a column vector of complex field amplitudes for JJ optical spatial modes on the source plane and |α⟩|\alpha\rangle as a multimode coherent state with amplitude α\alpha. For thermal sources, it is standard to assume Φ\Phi to be a zero-mean complex Gaussian given by Φ(α)=1det(πΓs)exp(−α†Γ−1α) \Phi(\alpha)=\frac{1}{\operatorname{det}(\pi \Gamma_s)} \exp \left(-\alpha^{\dagger} \Gamma^{-1} \alpha\right) where $\alpha^{\dagger}=\left(\alpha_1^*, \ldots, \alpha_J^*\right)$ denotes the complex transpose of α\alpha. We consider a simple two point source case \begin{equation} [\Gamma_s]_{uv}=\epsilon_0\delta_{uv}[\delta_{uu_1}+\delta_{uu_2}] \end{equation} where ϵ0→0\epsilon_0\to 0, and u1u_1, u2u_2 are the positions of the two point sources. If the quantum state ρs\rho_s passes through a single lens which has point spread function ψ(x)\psi(x). We want to quantify the performance of a naive direct imaging method and the fundamental limit for resolution. # Main problem: Define θ=13u1+23u2\theta=\frac{1}{3}u_1+\frac{2}{3}u_2, calculate the quantum Fisher information of estimating θ\theta per each measured photon, express the answer using Δk2≡∫−∞∞dx[∂ψ(x)∂x]2 \Delta k^2 \equiv \int_{-\infty}^{\infty} d x\left[\frac{\partial \psi(x)}{\partial x}\right]^2 , γ≡∫−∞∞dx∂ψ(x)∂xψ(x−u2+u1) \gamma \equiv \int_{-\infty}^{\infty} d x \frac{\partial \psi(x)}{\partial x} \psi\left(x-u_2+u_1 \right) .
Plain-text mathematical notation (without MathML)

# Problem setup:
In this problem, we want to consider the resolution of single-lens imaging of thermal sources. The quantum state radiated by the thermal sources can be expressed as

ρ_(s)=∫DαΦ(α)|α⟩⟨α|,




where Φ(α) is the Sudarshan-Glauber representation and Dα is an appropriate measure. Define  $\alpha=\left(\alpha_1, \ldots, \alpha_J\right)^{\top}$ as a column vector of complex field amplitudes for J optical spatial modes on the source plane and |α⟩ as a multimode coherent state with amplitude α.
 For thermal sources, it is standard to assume Φ to be a zero-mean complex Gaussian given by

Φ(α)=(1)/(det(πΓ_(s)))exp(−α^(†)Γ^(−1)α)

where $\alpha^{\dagger}=\left(\alpha_1^*, \ldots, \alpha_J^*\right)$ denotes the complex transpose of α.  We consider a simple two point source case

\begin{equation}
[\Gamma_s]_{uv}=\epsilon_0\delta_{uv}[\delta_{uu_1}+\delta_{uu_2}]
\end{equation}
where ϵ₀→0, and u₁, u₂ are the positions of the two point sources. If the quantum state ρ_(s) passes through a single lens which has point spread function ψ(x). We want to quantify the performance of a naive direct imaging method and the fundamental limit for resolution.

# Main problem:

Define θ=(1)/(3)u₁+(2)/(3)u₂, calculate the quantum Fisher information of estimating θ per each measured photon, express the answer using
Δk²≡∫_(−∞)^(∞)dx[(∂ψ(x))/(∂x)]²,
γ≡∫_(−∞)^(∞)dx(∂ψ(x))/(∂x)ψ(x−u₂+u₁).
Original LaTeX notation

# Problem setup:
In this problem, we want to consider the resolution of single-lens imaging of thermal sources. The quantum state radiated by the thermal sources can be expressed as

$$
\rho_s=\int D \alpha \Phi(\alpha)|\alpha\rangle\langle\alpha|,
$$




where $\Phi(\alpha)$ is the Sudarshan-Glauber representation and $D \alpha$ is an appropriate measure. Define  $\alpha=\left(\alpha_1, \ldots, \alpha_J\right)^{\top}$ as a column vector of complex field amplitudes for $J$ optical spatial modes on the source plane and $|\alpha\rangle$ as a multimode coherent state with amplitude $\alpha$.
 For thermal sources, it is standard to assume $\Phi$ to be a zero-mean complex Gaussian given by

$$
\Phi(\alpha)=\frac{1}{\operatorname{det}(\pi \Gamma_s)} \exp \left(-\alpha^{\dagger} \Gamma^{-1} \alpha\right)
$$

where $\alpha^{\dagger}=\left(\alpha_1^*, \ldots, \alpha_J^*\right)$ denotes the complex transpose of $\alpha$.  We consider a simple two point source case

\begin{equation}
[\Gamma_s]_{uv}=\epsilon_0\delta_{uv}[\delta_{uu_1}+\delta_{uu_2}]
\end{equation}
where $\epsilon_0\to 0$, and $u_1$, $u_2$ are the positions of the two point sources. If the quantum state $\rho_s$ passes through a single lens which has point spread function $\psi(x)$. We want to quantify the performance of a naive direct imaging method and the fundamental limit for resolution.

# Main problem:

Define $\theta=\frac{1}{3}u_1+\frac{2}{3}u_2$, calculate the quantum Fisher information of estimating $\theta$ per each measured photon, express the answer using
$
\Delta k^2 \equiv \int_{-\infty}^{\infty} d x\left[\frac{\partial \psi(x)}{\partial x}\right]^2
$,
$
\gamma \equiv \int_{-\infty}^{\infty} d x \frac{\partial \psi(x)}{\partial x} \psi\left(x-u_2+u_1 \right)
$.

code template

Code

import sympy as sp

Delta_k_sq, gamma = sp.symbols('Delta_k_sq gamma')

def answer(Delta_k_sq, gamma):
    r"""
    Return the expression of the quantum Fisher information in Sympy format.

    Inputs
    ----------
    Delta_k_sq: sympy.Symbol, $\Delta k^2 \equiv \int_{-\infty}^{\infty} d x\left[\frac{\partial \psi(x)}{\partial x}\right]^2$
    gamma: sympy.Symbol, $\gamma \equiv \int_{-\infty}^{\infty} d x \frac{\partial \psi(x)}{\partial x} \psi\left(x-u_2+u_1 \right)$

    Outputs
    ----------
    QFI: sympy.Expr, quantum Fisher information of estimating $\theta=\frac{1}{3}u_1+\frac{2}{3}u_2$ per each measured photon,
         $u_1$ and $u_2$ are the positions of the two point sources
    """

    # ------------------ FILL IN YOUR RESULTS BELOW ------------------
    QFI = ...  # a SymPy expression of inputs
    # ---------------------------------------------------------------

    return QFI

Discussion

Discussion

No discussion posts on this page yet. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.

See answer Answer published by the source

Artifacts

Code, notes and reproducible work shared by participants. Files are served from a separate origin.

No artifacts on this page yet. Share reproducible code or notes in a contribution. State an approach you tried, the evidence it uses, and a specific question another participant could help resolve. Use the posting template.

Source and history

Official source

initial import