-
Terahertz frequency combs exploiting an on-chip solution processed graphene-quantum cascade laser coupled-cavity architecture
Authors:
F. P. Mezzapesa,
K. Garrasi,
J. Schmidt,
L. Salemi,
V. Pistore,
L. Li,
A. G. Davies,
E. H. Linfield,
M. Riesch,
C. Jirauschek,
T. Carey,
F. Torrisi,
A. C. Ferrari,
M. S. Vitiello
Abstract:
The ability to engineer quantum-cascade-lasers (QCLs) with ultrabroad gain spectra and with a full compensation of the group velocity dispersion, at Terahertz (THz) frequencies, is a fundamental need for devising monolithic and miniaturized optical frequency-comb-synthesizers (FCS) in the far-infrared. In a THz QCL four-wave mixing, driven by the intrinsic third-order susceptibility of the intersu…
▽ More
The ability to engineer quantum-cascade-lasers (QCLs) with ultrabroad gain spectra and with a full compensation of the group velocity dispersion, at Terahertz (THz) frequencies, is a fundamental need for devising monolithic and miniaturized optical frequency-comb-synthesizers (FCS) in the far-infrared. In a THz QCL four-wave mixing, driven by the intrinsic third-order susceptibility of the intersubband gain medium, self-lock the optical modes in phase, allowing stable comb operation, albeit over a restricted dynamic range (~ 20% of the laser operational range). Here, we engineer miniaturized THz FCSs comprising a heterogeneous THz QCL integrated with a tightly-coupled on-chip solution-processed graphene saturable-absorber reflector that preserves phase-coherence between lasing modes even when four-wave mixing no longer provides dispersion compensation. This enables a high-power (8 mW) FCS with over 90 optical modes to be demonstrated, over more than 55% of the laser operational range. Furthermore, stable injection-locking is showed, paving the way to impact a number of key applications, including high-precision tuneable broadband-spectroscopy and quantum-metrology.
△ Less
Submitted 23 November, 2020;
originally announced November 2020.
-
Optoelectronic device simulations based on macroscopic Maxwell-Bloch equations
Authors:
Christian Jirauschek,
Michael Riesch,
Petar Tzenov
Abstract:
Due to their intuitiveness, flexibility and relative numerical efficiency, the macroscopic Maxwell-Bloch (MB) equations are a widely used semiclassical and semi-phenomenological model to describe optical propagation and coherent light-matter interaction in media consisting of discrete-level quantum systems. This review focuses on the application of this model to advanced optoelectronic devices, su…
▽ More
Due to their intuitiveness, flexibility and relative numerical efficiency, the macroscopic Maxwell-Bloch (MB) equations are a widely used semiclassical and semi-phenomenological model to describe optical propagation and coherent light-matter interaction in media consisting of discrete-level quantum systems. This review focuses on the application of this model to advanced optoelectronic devices, such as quantum cascade and quantum dot lasers. The Bloch equations are here treated as a density matrix model for driven quantum systems with two or multiple discrete energy levels, where dissipation is included by Lindblad terms. Furthermore, the one-dimensional MB equations for semiconductor waveguide structures and optical fibers are rigorously derived. Special analytical solutions and suitable numerical methods are presented. Due to the importance of the MB equations in computational electrodynamics, an emphasis is placed on the comparison of different numerical schemes, both with and without the rotating wave approximation. The implementation of additional effects which can become relevant in semiconductor structures, such as spatial hole burning, inhomogeneous broadening and local-field corrections, is discussed. Finally, links to microscopic models and suitable extensions of the Lindblad formalism are briefly addressed.
△ Less
Submitted 9 June, 2020;
originally announced June 2020.
-
Mbsolve: An open-source solver tool for the Maxwell-Bloch equations
Authors:
Michael Riesch,
Christian Jirauschek
Abstract:
The Maxwell-Bloch equations are a valuable tool to model light-matter interaction, where the application examples range from the description of pulse propagation in two-level media to the elaborate simulation of optoelectronic devices, such as the quantum cascade laser (QCL). In this work, we present mbsolve, an open-source solver tool for the Maxwell-Bloch equations. Here, we consider the one-dim…
▽ More
The Maxwell-Bloch equations are a valuable tool to model light-matter interaction, where the application examples range from the description of pulse propagation in two-level media to the elaborate simulation of optoelectronic devices, such as the quantum cascade laser (QCL). In this work, we present mbsolve, an open-source solver tool for the Maxwell-Bloch equations. Here, we consider the one-dimensional Maxwell's equations, which are coupled to the Lindblad equation. The resulting generalized Maxwell-Bloch equations are treated without invoking the rotating wave approximation (RWA). Since this full-wave treatment is computationally intensive, we provide a flexible framework to implement different numerical methods and/or parallelization techniques. On this basis, we offer two solver implementations that use OpenMP for parallelization.
△ Less
Submitted 24 July, 2021; v1 submitted 11 May, 2020;
originally announced May 2020.
-
bertha: Project Skeleton for Scientific Software
Authors:
Michael Riesch,
Tien Dat Nguyen,
Christian Jirauschek
Abstract:
Science depends heavily on reliable and easy-to-use software packages, such as mathematical libraries or data analysis tools. Developing such packages requires a lot of effort, which is too often avoided due to the lack of funding or recognition. In order to reduce the efforts required to create sustainable software packages, we present a project skeleton that ensures the best software engineering…
▽ More
Science depends heavily on reliable and easy-to-use software packages, such as mathematical libraries or data analysis tools. Developing such packages requires a lot of effort, which is too often avoided due to the lack of funding or recognition. In order to reduce the efforts required to create sustainable software packages, we present a project skeleton that ensures the best software engineering practices from the start of a project, or serves as reference for existing projects.
△ Less
Submitted 25 March, 2020; v1 submitted 3 December, 2019;
originally announced December 2019.
-
Analyzing the Positivity Preservation of Numerical Methods for the Liouville-von Neumann Equation
Authors:
Michael Riesch,
Christian Jirauschek
Abstract:
The density matrix is a widely used tool in quantum mechanics. In order to determine its evolution with respect to time, the Liouville-von Neumann equation must be solved. However, analytic solutions of this differential equation exist only for simple cases. Additionally, if the equation is coupled to Maxwell's equations to model light-matter interaction, the resulting equation set -- the Maxwell-…
▽ More
The density matrix is a widely used tool in quantum mechanics. In order to determine its evolution with respect to time, the Liouville-von Neumann equation must be solved. However, analytic solutions of this differential equation exist only for simple cases. Additionally, if the equation is coupled to Maxwell's equations to model light-matter interaction, the resulting equation set -- the Maxwell-Bloch or Maxwell-Liouville-von Neumann (MLN) equations -- becomes nonlinear. In these advanced cases, numerical methods are required. Since the density matrix has certain mathematical properties, the numerical methods applied should be designed to preserve those properties. We establish the criterion that only methods that have a completely positive trace preserving (CPTP) update map can be used in long-term simulations. Subsequently, we assess the three most widely used methods -- the matrix exponential (ME) method, the Runge-Kutta (RK) method, and the predictor-corrector (PC) approach -- whether they provide this feature, and demonstrate that only the update step of the matrix exponential method is a CPTP map.
△ Less
Submitted 10 April, 2019; v1 submitted 1 August, 2018;
originally announced August 2018.
-
Numerical Method for the Maxwell-Liouville-von Neumann Equations using Efficient Matrix Exponential Computations
Authors:
Michael Riesch,
Christian Jirauschek
Abstract:
We present a novel method to solve the Maxwell-Liouville-von Neumann (MLN) equations in an accurate and efficient way without invoking the rotating wave approximation (RWA). The method is a combination of two established concepts, namely the operator splitting method as well as the adjoint representation of the Lie algebra SU(N) (or pseudospin representation). The former concept ensures the accura…
▽ More
We present a novel method to solve the Maxwell-Liouville-von Neumann (MLN) equations in an accurate and efficient way without invoking the rotating wave approximation (RWA). The method is a combination of two established concepts, namely the operator splitting method as well as the adjoint representation of the Lie algebra SU(N) (or pseudospin representation). The former concept ensures the accuracy of the approach, but is computationally expensive. The latter concept provides an efficient representation of the problem and two optimization possibilities. We have implemented and verified both optimization approaches and demonstrate that substantial speedup can be achieved.
△ Less
Submitted 26 October, 2017;
originally announced October 2017.