[Paper Review] Numerical simulation of vacuum particle production: applications to cosmology, dynamical Casimir effect and time-dependent non-homogeneous dielectrics

Nuno D. Antunes|arXiv (Cornell University)|Oct 10, 2003

Quantum Electrodynamics and Casimir Effect2 references3 citations

TL;DR

This paper presents a general numerical method for simulating vacuum particle production in time- and space-dependent quantum field systems, applying it to cosmological models, the dynamical Casimir effect with oscillating mirrors, and expanding dielectric bubbles. The method achieves high accuracy by dynamically tracking moving boundaries and time-dependent optical properties, yielding precise particle spectra and production rates, with results confirming prior estimates and revealing new dependencies on expansion velocity and dielectric strength.

ABSTRACT

We develop a general numerical method aimed at studying particle production from vacuum states in a variety of settings. As a first example we look at particle production in a simple cosmological model. We apply the same approach to the dynamical Casimir effect, with special focus on the case of an oscillating mirror. We confirm previous estimates and obtain long-time production rates and particle spectra for both resonant and off-resonant frequencies. Finally, we simulate a system with space and time-dependent optical properties, analogous to a one-dimensional expanding dielectric bubble. We obtain simple expressions for the dependence of the final particle number on the expansion velocity and final dielectric constant.

Motivation & Objective

• To develop a general numerical framework for studying vacuum particle production in non-trivial spacetime and material backgrounds.
• To accurately simulate particle creation in cosmological models with time-dependent metrics.
• To model the dynamical Casimir effect with oscillating mirrors, including resonant and off-resonant frequencies.
• To investigate particle production in systems with space- and time-dependent dielectric properties, such as expanding dielectric bubbles.
• To quantify the dependence of final particle number on expansion velocity and dielectric strength in non-homogeneous media.

Proposed method

• A finite-difference time-domain (FDTD) scheme is used to solve the scalar field equation in curved or time-dependent backgrounds.
• The method dynamically tracks the position of moving boundaries (e.g., mirrors) using interpolation to estimate field and momentum values at non-lattice points.
• Boundary conditions are enforced by introducing auxiliary fictitious lattice points at the mirror's intersection with the grid, improving accuracy in both time and space.
• Spatial derivatives near the boundary are corrected using a parameter α that accounts for the mirror’s fractional position between lattice sites.
• Momentum updates at the boundary use interpolated field values and velocity information, with a time-shift parameter β to align with the mirror’s trajectory.
• The algorithm is validated against analytical solutions for uniformly expanding cavities, showing significant error reduction when dynamic tracking is applied.

Experimental results

Research questions

• RQ1 How does particle production in a cosmological model with an expanding metric compare to analytical predictions?
• RQ2 What are the long-time particle production rates and spectra in the dynamical Casimir effect with an oscillating mirror?
• RQ3 How does the particle yield depend on the expansion velocity and final dielectric constant in a time- and space-dependent dielectric medium?
• RQ4 What is the impact of numerical boundary treatment on the accuracy of particle production calculations?
• RQ5 Can a fully numerical method reproduce known analytical results for particle creation in time-dependent systems?

Key findings

• The numerical method accurately reproduces analytical results for particle production in a uniformly expanding cavity, validating its reliability.
• For the oscillating mirror case, the method confirms previous estimates of particle production rates and spectra, including both resonant and off-resonant frequencies.
• The long-time particle production rate in the dynamical Casimir effect is found to be stable and consistent with theoretical expectations.
• In the expanding dielectric bubble model, the final particle number scales non-trivially with expansion velocity and dielectric strength, with simple power-law dependencies derived from numerical data.
• The use of dynamic interpolation (with α and β) reduces numerical errors by a factor of 25 compared to standard finite-difference schemes without correction.
• High-frequency discontinuities in the field and momentum are significantly reduced by the improved boundary treatment, leading to more accurate Bogoliubov coefficient calculations.