Benchmark AI / 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.

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

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 ϵ00\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 Δk2dx[ψ(x)x]2 \Delta k^2 \equiv \int_{-\infty}^{\infty} d x\left[\frac{\partial \psi(x)}{\partial x}\right]^2 , γdxψ(x)xψ(xu2+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)
$.

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.

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