quantumaharonovbohm2d_tests.cpp

Go to the documentation of this file.

/**
 * @file quantumaharonovbohm2d_tests.cpp
 * @brief Implements regression and invariant tests for quantumaharonovbohm2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/quantumaharonovbohm2d_tests.cpp. This file is part of Quantum Phenomena Explorer (QPE),
 * maintained by Officina Turini. See LICENSE_IT.txt / LICENSE.txt for the
 * QPE freeware source-available terms and THIRD_PARTY_NOTICES.txt for
 * third-party components.
 */
 
#include "quantumaharonovbohm2d_tests.h"
 
#include "core/math/complexfield2d.h"
#include "core/quantum/gaussianpacket2d.h"
#include "core/quantum/potentialfield2d.h"
#include "core/quantum/quantumaharonovbohm2d.h"
#include "core/quantum/splitstepfouriersolver2d.h"
#include "core/quantum/u1gaugefield2d.h"
 
#include <cmath>
#include <complex>
 
namespace qpe { namespace tests {
 
/**
 * @brief Runs quantum aharonov bohm2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runQuantumAharonovBohm2DTests(std::string &error)
{
    const double pi = 3.14159265358979323846;
    ComplexField2D field(64u, 64u, 0.20f, 0.20f);
    GaussianPacket2DParams packet;
    packet.centerX = -3.0f;
    packet.centerY = 0.0f;
    packet.sigma = 1.2f;
    packet.momentumX = 2.0f;
    packet.momentumY = 0.0f;
    packet.hbar = 1.0f;
    GaussianPacket2D::initialize(field, packet);
 
    QuantumAharonovBohm2DConfig config;
    config.enabled = true;
    config.solenoidRadius = 0.8;
    config.barrierHeight = 100.0;
    config.detectorX = 2.5;
    config.analysisEndTime = 2.0;
    config.fluxQuanta = 0.0;
 
    PotentialField2D potential;
    U1GaugeField2D gauge;
    if (!QuantumAharonovBohm2D::buildFields(field, config, potential, gauge, &error)) {
        return false;
    }
    if (!potential.matchesGrid(field) || !gauge.matchesGrid(field) ||
        potential.maximum() < 99.0f) {
        error = "Aharonov-Bohm scalar/gauge field build failed";
        return false;
    }
    for (std::size_t y = 0; y < gauge.height(); y += 7u) {
        for (std::size_t x = 0; x < gauge.width(); x += 7u) {
            if (std::abs(gauge.linkX(x, y) - U1GaugeField2D::Link(1.0f, 0.0f)) > 1.0e-6f ||
                std::abs(gauge.linkY(x, y) - U1GaugeField2D::Link(1.0f, 0.0f)) > 1.0e-6f) {
                error = "zero-flux U(1) links are not identity";
                return false;
            }
        }
    }
 
    config.fluxQuanta = 0.5;
    if (std::fabs(QuantumAharonovBohm2D::topologicalPhaseRadians(config) - pi) > 1.0e-12) {
        error = "half-flux Aharonov-Bohm phase is not pi";
        return false;
    }
    if (!QuantumAharonovBohm2D::buildFields(field, config, potential, gauge, &error)) return false;
    QuantumAharonovBohm2D::applyIncidentGaugePhase(field, config);
    if (std::fabs(field.norm2() - 1.0) > 2.0e-5) {
        error = "incident gauge phase changed wavefunction norm";
        return false;
    }
 
    SplitStepFourierSolver2D solver;
    SplitStepFourier2DParams params;
    params.mass = 1.0f;
    params.hbar = 1.0f;
    params.dt = 0.01f;
    solver.setParameters(params);
    solver.prepare(field);
    solver.preparePotential(potential);
    for (int i = 0; i < 20; ++i) solver.stepGaugeCovariant(field, potential, gauge);
    if (std::fabs(field.norm2() - 1.0) > 2.0e-3) {
        error = "gauge-covariant Cayley propagation did not preserve norm";
        return false;
    }
 
    QuantumAharonovBohmHistory2D history;
    if (!history.configure(field, config)) {
        error = "Aharonov-Bohm history configuration failed";
        return false;
    }
    history.append(solver.time(), field);
    if (history.empty() || history.samples().back().detectorProfileY.size() != field.height()) {
        error = "Aharonov-Bohm detector history append failed";
        return false;
    }
 
    error.clear();
    return true;
}
 
} }