narp-dynamics 1.4

Implementation of NARP Pulses and Their Interaction Dynamics with Two-Level Systems

from narp_functions import (
    narp_pulse,
    narp_pulse_dynamics,
    area_versus_chirp,
    notch_versus_chirp,
    notch_versus_area,
    plot_fields,
    plot_occupations,
    plot_surfaces,
)
import jax.numpy as jnp

Basic Usage

pulse = narp_pulse(
    pulse_center_frequency = 0.322, #[PHz]
    pulse_bandwidth = 0.022, #[PHz] (20 fs pulse)
    pulse_area = 6*jnp.pi,
    chirp_rate = 900, #[ps]^2
    notch_type = 0,
    notch_bandwidth = 0.022/5 #[PHz]
)
plot_fields(pulse)

rho = narp_pulse_dynamics(pulse, 0.322) # (Resonance)
plot_occupations(pulse, rho)

Parameter Space

Area versus Chirp

pulse_areas_grid, chirp_rates_grid, occupations_matrix = area_versus_chirp(
    pulse_center_frequency = 0.322, #[PHz]
    pulse_bandwidth = 0.022, #[PHz] (20 fs pulse)
    notch_type = 0,
    notch_bandwidth = 0.022/5, #[PHz]
    resolution = 25
)
plot_surfaces(chirp_rates_grid, pulse_areas_grid, occupations_matrix)

Computing occupations: 100%|██████████| 625/625 [03:04<00:00,  3.38it/s]

Notch versus Chirp

pulse_areas_grid, notch_widths_grid, occupations_matrix = notch_versus_chirp(
    pulse_center_frequency = 0.322, #[PHz]
    pulse_bandwidth = 0.022, #[PHz] (20 fs pulse)
    notch_type = 0,
    pulse_area = 8*jnp.pi,
    resolution = 25
)
plot_surfaces(notch_widths_grid, pulse_areas_grid, occupations_matrix)

Computing occupations: 100%|██████████| 625/625 [02:57<00:00,  3.52it/s]

Notch versus Area

pulse_areas_grid, notch_widths_grid, occupations_matrix = notch_versus_area(
    pulse_center_frequency = 0.322, #[PHz]
    pulse_bandwidth = 0.022, #[PHz] (20 fs pulse)
    chirp_rate = 600,
    notch_type = 0,
    resolution = 25
)
plot_surfaces(notch_widths_grid, pulse_areas_grid, occupations_matrix)

Computing occupations: 100%|██████████| 625/625 [02:53<00:00,  3.61it/s]

This project is licensed under the MIT License - see the LICENSE file for details.