Struct PathIntegralMonteCarlo 

pub struct PathIntegralMonteCarlo {
    pub n_slices: usize,
    pub tau: f64,
    pub mass: f64,
    pub path: Vec<f64>,
}

A simple Path Integral Monte Carlo estimator for 1D quantum mechanics.

Uses a discrete Euclidean path (imaginary time) to estimate the ground state energy via the virial estimator.

Fields

n_slices: usize

Number of time slices.

tau: f64

Imaginary time step τ = β/N.

mass: f64

Particle mass.

path: Vec<f64>

Current path (positions at each time slice).

Implementations

impl PathIntegralMonteCarlo

pub fn new(n_slices: usize, beta: f64, mass: f64) -> Self

Create a new PIMC estimator with all positions set to zero.

pub fn kinetic_action(&self) -> f64

Euclidean (kinetic) action for the current path. S_E = Σ_k m/(2τ) * (x_{k+1} − x_k)²

pub fn potential_action_harmonic(&self, omega: f64) -> f64

Harmonic potential action contribution. S_V = τ * Σ_k ½ ω² x_k²

pub fn total_action_harmonic(&self, omega: f64) -> f64

Total Euclidean action (kinetic + harmonic potential).

pub fn virial_energy_estimate(&self, omega: f64) -> f64

Estimate the ground state energy via the virial theorem (harmonic case). E_0 ≈ ½ω² ⟨x²⟩ + ½/(2m τ²) corrections (simplified).

pub fn set_path(&mut self, path: Vec<f64>)

Set the path to a given configuration.

pub fn winding_number(&self) -> f64

Compute the winding number (only meaningful for a ring).

Trait Implementations

impl Clone for PathIntegralMonteCarlo

fn clone(&self) -> PathIntegralMonteCarlo

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for PathIntegralMonteCarlo

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.