Introducing an Extensible Open-Source Toolkit Suite for Studying Second Harmonic Generation: A Case Study of Depleted Pulsed Gaussian Wave SHG

This paper introduces an extensible, open-source SHG Computational Toolkit Suite designed to overcome the limitations of existing analytical models and inaccessible experimental data by providing a coordinated collection of well-documented numerical tools for studying complex, thermally coupled second harmonic generation scenarios.

The Gist

Imagine you are trying to bake the perfect cake (creating a new color of light) using a very specific recipe (a laser crystal). For a long time, scientists have been trying to figure out exactly how this "baking" works. However, most of the old recipes were written with big, simple assumptions—like pretending the oven never gets hot or that the ingredients never run out. In reality, the oven gets hot, the ingredients change, and the process is messy and complicated.

This paper introduces a new, open-source "Digital Kitchen" (a software toolkit) that helps scientists simulate this process with much greater accuracy. Here is a breakdown of what they did, using simple analogies:

1. The Problem: The "Black Box" of Light

When you shoot a laser through a special crystal, it can double the frequency of the light, turning red light into green light (or infrared into visible light). This is called Second Harmonic Generation (SHG).

• The Old Way: Scientists used math formulas that were like "flat maps" of a mountain. They worked okay for simple hills, but they failed to capture the steep cliffs and deep valleys of real-world physics, especially when heat builds up inside the crystal.
• The Experimental Problem: To fix the math, you'd need to measure the temperature inside the crystal at every single point while the laser is firing. But you can't stick a thermometer inside a laser beam without breaking the experiment. It's like trying to measure the exact temperature of a soufflé while it's rising without opening the oven door.

2. The Solution: The "LEGO" Toolkit Suite

The authors built a Computational Toolkit Suite. Think of this not as one giant, unchangeable machine, but as a box of high-quality LEGO bricks.

• Modular: Each brick is a small, independent tool that handles a specific part of the physics (like heat, or different beam shapes).
• Extensible: If a scientist wants to study a new type of laser, they don't have to build a whole new factory. They just snap in a new LEGO brick or rearrange the existing ones.
• Open-Source: The blueprints (code) are free for anyone to see, use, and modify. This stops everyone from reinventing the wheel.

3. The Case Study: The "Depleted" Wave

To prove their new LEGO set works, they built a specific model: a pulsed Gaussian wave.

• The Analogy: Imagine a powerful water hose (the laser pulse) spraying into a sponge (the crystal).
• The "Depleted" Part: In simple models, people assume the hose keeps spraying water at the same strength all the way through. But in reality, as the sponge soaks up the water to create a new effect (the second harmonic), the hose runs dry. The water pressure drops. This is called a "depleted" pump.
• The Simulation: The authors used a method called Finite Difference Method (FDM). Imagine the crystal is a 3D grid of tiny boxes. The computer calculates what happens in each box, step-by-step, as the pulse moves through. It tracks how the "water" (fundamental light) turns into "steam" (second harmonic light) and how the pressure drops as it goes.

4. What They Found

Using their new toolkit, they simulated a specific scenario (Type II SHG in a KTP crystal) with a pulse of light that lasts 50 microseconds.

• The Result: They watched the energy transfer happen in real-time on the computer. They saw that as the pulse traveled about 5 millimeters into the crystal, almost all the original light energy was converted into the new color.
• The "Depletion" Confirmed: The original beam didn't stay strong; it got "depleted" (ran out of energy) as it gave its power to the new beam.
• The Shape: Even though the energy changed, the new beam kept the same smooth, round "Gaussian" shape as the original laser, just like a shadow that changes color but keeps its outline.

5. Why This Matters

The paper claims that this toolkit allows researchers to:

• Replicate: Run the exact same simulation to check results.
• Adapt: Tweak the settings (like changing the pulse energy or the crystal type) without rewriting the whole code.
• Extend: Add new features, like heat effects, later on.

In short, the authors didn't just solve one specific problem; they built a universal workshop where scientists can now test complex light-behavior scenarios that were previously too hard to calculate or impossible to measure directly. They proved the workshop works by successfully simulating a "running out of fuel" scenario for a laser pulse, showing exactly how the energy transforms as it travels through the crystal.

Technical Summary

Technical Summary: Introducing an Extensible Open-Source Toolkit Suite for Studying Second Harmonic Generation

Problem Statement
While Second Harmonic Generation (SHG) in nonlinear crystals is a well-established area of research, accurate modeling remains challenging in realistic regimes involving thermal effects and complex beam configurations. Existing analytical models often rely on simplifying assumptions that limit their applicability under realistic operating conditions. Furthermore, direct experimental characterization of thermal effects is impractical because it requires spatiotemporal temperature data at every point within the crystal, which is not experimentally accessible. Although numerous computational SHG models exist, there is a notable lack of a fully accessible, well-documented, and extensible resource that consolidates validated methods across diverse SHG configurations. Current tools often fail to support reproducible workflows, modular extension, or adaptation to complex nonlinear and thermo-optical scenarios.

Methodology
To address these gaps, the authors developed the SHG Computational Toolkit Suite, a coordinated collection of independent modeling toolkits hosted as a public GitHub Organization. The suite is designed with a modular architecture, allowing researchers to combine and adapt independent toolkits for specific scenarios. It provides comprehensive documentation of underlying physics and numerical implementation details, validated implementations for depleted regimes and pulsed configurations, and full source-code accessibility.

As a case study to demonstrate the suite's extensibility, the authors modeled depleted pulsed Gaussian wave SHG under a Type II configuration in a KTP crystal. The specific methodology employed includes:

• Governing Equations: The model derives from Maxwell's equations, utilizing the Helmholtz equation for a dispersive and nonlinear medium. Three coupled equations describe the evolution of two fundamental beams (ordinary and extraordinary) and the generated second-harmonic wave.
• Approximations: The study assumes ideal coupling, perfect phase matching ($\Delta k = 0$), and negligible thermal absorption for the primary demonstration, though the framework allows for the inclusion of these effects.
• Numerical Approach: The equations are solved using a Finite Difference Method (FDM) in cylindrical coordinates to exploit azimuthal symmetry. The discretization employs backward FDM for temporal derivatives, forward FDM for longitudinal spatial derivatives, and central FDM for radial derivatives.
• Boundary and Initial Conditions: The simulation uses Gaussian spatial profiles and Gaussian temporal pulse profiles (with a duration $t_p = 50 \mu s$) as boundary conditions at the crystal input ($z=0$). Initial conditions assume zero fields before the pulse arrival.
• Implementation: The solution is computed using a self-developed FORTRAN code running on Linux Ubuntu, utilizing a non-uniform radial grid to optimize efficiency.

Key Contributions

1. SHG Computational Toolkit Suite: The primary contribution is the release of a publicly accessible, modular, and documented software suite that consolidates validated SHG modeling methods. This allows researchers to replicate, adapt, and extend methods without reconstructing foundational computational infrastructure.
2. Validated Case Study: The paper presents a specific implementation of the toolkit for Type II pulsed Gaussian SHG in KTP crystals. This serves as a concrete template for adapting the suite to new research problems.
3. Reproducibility: By providing full source code, detailed implementation notes, and convergence criteria, the authors enable full reproducibility of the numerical results.

Results
The case study simulations yielded the following observations under idealized conditions (perfect phase matching, no absorption):

• Depleted Pump Dynamics: The model successfully captures the transfer of energy from the fundamental wave (FW) to the second harmonic wave (SHW). The FW efficiency drops steadily to zero along the crystal axis, while the SHW efficiency increases, indicating near-complete conversion.
• Interaction Length: Under the tested parameters (Pulse energy $E=0.45$ J, spot size $\omega_f = 80 \mu m$), nearly 90% of the FW energy is converted within approximately 5 mm of the crystal length ($L_c = 2$ cm). This short interaction length validates the necessity of a depleted-pump formalism over constant-beam approximations.
• Pulse Energy Dependence: Simulations varying pulse energy (from 0.1 J to 1.0 J) demonstrated that higher energies lead to shorter interaction lengths and faster energy conversion.
• Spatial Profiles: The generated SHW retains a transverse Gaussian profile at the output surface, consistent with the input fundamental beam profile.
• Absorption Effects: Comparisons between ideal (no absorption) and non-ideal (with absorption) scenarios showed that absorption causes a slight decrease in efficiency after the initial conversion region.

Significance
The paper positions the SHG Computational Toolkit Suite as a practical foundation for the research community. Its significance lies in:

• Accelerating Research: By providing pre-validated, modular components, the suite eliminates the need for researchers to rebuild computational infrastructure from scratch, thereby accelerating investigations into diverse SHG phenomena.
• Promoting Reproducibility: The open-source nature and comprehensive documentation address the current gap in accessible, reproducible SHG modeling resources.
• Educational Utility: The suite supports computationally supported physics education by providing executable examples that allow users to examine governing equations, run simulations, and modify parameters.
• Extensibility: The modular design facilitates future extensions to more complex scenarios, including heat generation, phase mismatch, and strongly coupled nonlinear dynamics, as well as adaptation to different beam types (e.g., Bessel-Gaussian) and crystal configurations.

The authors conclude that while this study focuses on a specific depleted pulsed Gaussian case, the toolkit's architecture is adaptable to a wide range of nonlinear optical problems, offering a reusable framework for future computational studies.