Abstract
We report on design, production and characterization of an extremely broadband multilayer beamsplitter, covering wavelength range from 0.67 – 2.6 µm. The group delay dispersion has small oscillations in the above mentioned working range. We used a new algorithm with floating constants allowing us to obtain a smooth and near constant GDD. The optical element based on the beamsplitter is used for dividing a low-energy super-octave spectrum into several sub-spectral regions which are later amplified and coherently combined.
© 2016 Optical Society of America
1. Introduction
To date, advances in attosecond technology is toward generation of isolated attosecond pulses at higher flux and higher photon energies, motivated by several ground-breaking experiments, like imaging electrondynamics in solids at attosecond time-scale, and attosecond-pump attosecond-probe spectroscopy [1,2].
The cutoff energy in high harmonic generation (HHG) can be extended to photon energies above the current state of the art [3], if high-energy single-cycle laser pulses, or sub-cycle light transients in the near-infrared (NIR) or the mid-infrared (MIR) spectral ranges become available [4–6]. Generation of these high-energy waveforms calls for a technique to overcome limitations in energy scalability of the current Ti:Sa laser technology and bandwidth limitations in few-cycle optical parametric chirped pulse amplifiers (OPCPA).
OPCPA-based synthesizers keep promise [4] for generating high-energy pulses with several octave spectral bandwidths. Here a low-energy super-octave spectrum is divided into several sub-spectral regions using a broadband dielectric beamsplitter. Subsequently each spectral region is amplified in OPCPA chains individually and finally the amplified spectra are compressed to their Fourier transform limit by using broadband chirped mirrors and combined coherently with a broadband dielectric beam combiner [4,7].
In this concept design and production of dielectric beamsplitters and beam combiners used to handle such super-octave continua, are immensely complicated and require using effective numerical optimization and accurate deposition tools.
In this paper we report about the design, production and characterization of a sophisticated 96-layers beamsplitter (BS) for dividing a super-octave continuum shown in [8] into two spectral regions which will be amplified by using different harmonics of an Yb:YAG thin disk laser [9] and later coherently combined in the OPCPA-based synthesizer discussed in [4]. To make the synthesizer’s geometry simple, the BS is designed to work at the incident angle (AOI) of 45°. In this design the broadband spectrum spanning from 670 nm to 1320 nm, is transmitted, while 1580 nm−2600 nm spectral components are reflected from the BS. Flat group delay dispersion on reflection (GDD) is required by design specifications and the designed BS provides GDD value equal to −235 fs2. The setup is presented in Fig. 1.
The primary goal of this paper is to demonstrate the feasibility of the complicated multilayer coating. The secondary aim is to point out the challenges in the design-production-characterization chain of complex multilayers in the very broadband spectral range. In Section 2 we consider the theoretical designing of BS and describe produced samples as well as measurement data. In Section 3 we provide optical characterization of the produced samples. Our conclusions are presented in Section 4.
2. Design and production
In our BS we used Nb2O5 as high index and SiO2 as low index materials. We chose two types of the substrates: B260 Glass of 1 mm thickness and Suprasil of 6.35 mm thickness. Nominal wavelength dependencies of refractive indices of thin film materials are described by Cauchy formula:
where are dimensionless parameters, nm, is specified in nanometers. For Nb2O5 ; for SiO2 ; for B260 Glass ; and for Suprasil .Design of BS was performed with the help of the needle optimization technique incorporated into OptiLayer software [10,11]. A merit function estimating the closeness of the designed spectral characteristics to the target specifications is defined as:
where is a vector of layer thicknesses, and are evenly distributed wavelength points in the spectral ranges from 670 nm to 1320 nm and from 1580 nm to 2600 nm, respectively. Since design specifications required just flat GDD, a constant is included in Eq. (2), and we used the floating constants approach [12] in order to find a solution with an optimal value of this constant.As a result a 96-layer BS was obtained with fs2. The BS refractive index profile is shown in Fig. 2(a). Oblique incidence theoretical reflectance and (p-polarization) in the spectral ranges of interest are shown in Fig. 2(b). Physical thickness of the obtained design is 13586 nm. The design contains eight layers with thicknesses less than 20 nm. Deposition of such a design is a complicated task requiring not only advanced instrumental tools but also a lot of experience. We produced the designed BS using Leybold Optics magnetron sputtering Helios plant, layer thicknesses were controlled using well-calibrated time monitoring [13]. The plant is equipped with two proprietary TwinMags magnetrons and a plasma source for plasma/ion-assisted reactive middle frequency dual magnetron sputtering. The plant is also equipped with broadband monitoring (BBM) system [14] which was used in a passive mode for data acquisition only. In the deposition run we used the Suprasil substrate for the BS optical element fabrication and the cheaper B260 Glass substrate as a witness sample. As a result, two samples were obtained: BS coating on B260 Glass substrate (sample BS-B260) and BS coating on Suprasil substrate (BS-Suprasil). The transmittance scans (BBM data) were recorded after deposition of each layer of BS-B260 sample.
After the deposition, transmittance data of produced samples were measured by Perkin Elmer Lambda 950 spectrophotometer (PE) in the range from 400 nm to 2500 nm with 0.5 nm wavelength step. Reflectance of the unpolarized light () and reflectance of the p-polarized light () at AOI = 45° were taken by the Universal Reflectance Accessory (URA) in the same spectral range. Group delay (GD) and GDD of the BS-Suprasil sample were extracted from white light interferometer (WLI) data. The infra-red detector was used, AOI = 45°, light was p-polarized.
From the experimental data it was found that the reflectance is more than 89% in the spectral range from 1580 nm to 2600 nm and is less than 4% in the range from 670 nm to 1320 nm. The experimental data show that produced BS samples satisfy target specifications and can be used in the OPCPA system shown in Fig. 1.
3. Optical characterization
3.1 Processing of BBM data
In order to improve the coating production in the IR spectral range we perform the characterization of the BS samples. We start with estimating errors in layer thicknesses on the basis of BBM data. We apply the algorithm assuming quasi-random errors in layer thicknesses [11,15]. The algorithm is based on the minimization of in situ discrepancy function with respect to relative errors :
where is the vector of first thicknesses of the model coating. The second term in Eq. (3) is responsible for the reliability of the solution taking into account a priori information that relative random errors are expected to be quite low. We applied a characterization approach based on the regularization theory of solving ill-posed problems [16]. In Eq. (3) is a regularization parameter that should be selected in order to provide a balance between data fitting and stability of the solution. In our characterization process we chose .Excellent fitting of the BBM data by model data was achieved in the course of characterization process, Fig. 3(a). Estimated errors do not exceed 2.5% that is in full agreement with expected accuracy of time monitoring [17,18]. Denote the thicknesses of the model coating obtained in this way.
3.2 Indices correction
In Fig. 3(b) we compare experimental transmittance in the spectral range from 1050 nm to 2500 nm related to BS-B260 and transmittance curve related to the model coating obtained above. We observe that the correspondence is worse than in the range from 400 to 1050 nm. The deviations can be addressed to inaccuracies in the wavelength dependencies of nominal refractive indices in the IR range. At the next step of the characterization process we correct dispersion curves of Nb2O5 and SiO2. According to Ref [19]. we assume that refractive index of Nb2O5 in the broad range from 400 nm to 2500 nm can be described by the Cauchy formula (see Eq. (1)). In the case of SiO2 three-parameter Cauchy model may not provide an accurate description of the dispersion behavior in the broad spectral range from 400 nm to 2600 nm. The wavelength dependence of the SiO2 refractive index having inflection of curve (see Ref [20], Fig. 1) can be only described by more flexible six-parameter Sellmeier model:
where are dimensionless parameters, nm, is specified in nanometers.We introduce a discrepancy function in the standard way:
where is a vector of model parameters, is specified by Cauchy formula Eq. (1), , is specified by Sellmeier coefficients taken from Ref [20], is the refractive index offset, are measured normal incidence transmittance data, . We minimize discrepancy function Eq. (5) with respect to Cauchy parameters and the offset . We achieved a good fitting of experimental transmittance by the model data in the range from 1050 nm to 2500 nm. In order to obtain refractive indices in the whole spectral range from 400 nm to 2500 nm we merge nominal indices taken from the range from 400 nm to 1050 nm and estimated indices in the range from 1050 to 2500 nm. As the curve at the junction point is not continuous, we smooth the resulted wavelength dependencies by Cauchy dependence for Nb2O5 and Sellmeier dependence for SiO2. The corrected and nominal refractive indices are compared in Fig. 4(a). In this figure one can observe that the deviations between corrected and nominal refractive indices are about 0.02. While Nb2O5 refractive index mainly changed by some offset value, SiO2 refractive index changed qualitatively and the new dependence has an inflection point around 1.2-1.3 µm. It is seen also that the deviations between corrected and nominal refractive indices in the range from 400 nm to 1050 nm are almost negligible. Denote the corrected refractive indices of the model coating. In Fig. 4(b) one can observe an excellent correspondence between experimental and model T/R characteristics in the whole spectral range from 400 nm 2500 nm.3.3 Characterization of BS-Suprasil sample
In the course of the deposition process the deviations between layer thicknesses of the samples placed on the witness and on a calotte position are inevitable. The thickness non-uniformity, however, cannot exceed 1% for Helios machine. Hence, we have to correct thickness errors found for BS-B260 model coating in order to estimate layer errors in BS-Suprasil coating. We assume presence of a systematic shift in layer thicknesses and define a vector of model coating parameters. The systematic shift can be estimated on the basis of the minimization of discrepancy function (Eq. (5)) with experimental data related to the BS-Suprasil sample. The estimated systematic shift is equal to −0.8%. It is remarkable that the resulted errors in BS-Suprasil coating calculated as do not exceed 1.7%. We consider a model coating found in this way as a characterization solution for BS-Suprasil sample. In Fig. 5(a), we observe a good correspondence between measured oblique incidence reflectance as well as with the corresponding model curves. Note that we use GD and GDD data for the verification purposes only, since typically they have lower accuracy and their scales are very different from spectrophotometric data. is For additional verification of the results we compare initial discrepancy function values and discrepancies achieved for final characterization solutions (see Table 1). In this table one can see that in all cases, the final discrepancy function values are decreased in comparison the initial discrepancies. This can be considered as an indication of the reliability of the obtained characterization results.
In Fig. 5(b) we show the simulated envelopes of the input Fourier-limited pulse and pulse reflected from the BS. It is worse to note that the shape of the reflected pulse is very close to the Fourier-limit one.
4. Conclusions
An outstanding 96-layer beamsplitter for the OPCPA system was designed, produced and characterized. For the first time we have applied floating GDD constant algorithm and were able to design smooth and constant GDD in the broadband working range. Good agreement between experimental and theoretical spectral performances demonstrates the applicability of the beamsplitter in the OPCPA systems.
Funding
DFG Cluster of Excellence, “Munich Centre for Advanced Photonics,” (http://www.munich-photonics.de); European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie agreement (657596).
References and lnks
1. S. R. Leone, C. W. McCurdy, J. Burgdörfer, L. S. Cederbaum, Z. Chang, N. Dudovich, J. Feist, C. H. Greene, M. Ivanov, R. Kienberger, U. Keller, M. F. Kling, Z.-H. Loh, T. Pfeifer, A. N. Pfeiffer, R. Santra, K. Schafer, A. Stolow, U. Thumm, and M. J. J. Vrakking, “What will it take to observe processes in “real time”?” Nat. Photonics 8(3), 162–166 (2014). [CrossRef]
2. F. Krausz and M. Ivanov, “Attosecond physics,” Rev. Mod. Phys. 81(1), 163–234 (2009). [CrossRef]
3. W. Schweinberger, A. Sommer, E. Bothschafter, J. Li, F. Krausz, R. Kienberger, and M. Schultze, “Waveform-controlled near-single-cycle milli-joule laser pulses generate sub-10 nm extreme ultraviolet continua,” Opt. Lett. 37(17), 3573–3575 (2012). [CrossRef] [PubMed]
4. H. Fattahi, H. G. Barros, M. Gorjan, T. Nubbemeyer, B. Alsaif, C. Y. Teisset, M. Schultze, S. Prinz, M. Haefner, M. Ueffing, A. Alismail, L. Vámos, A. Schwarz, O. Pronin, J. Brons, X. T. Geng, G. Arisholm, M. Ciappina, V. S. Yakovlev, D.-E. Kim, A. M. Azzeer, N. Karpowicz, D. Sutter, Z. Major, T. Metzger, and F. Krausz, “Third-generation femtosecond technology,” Optica 1(1), 45–63 (2014). [CrossRef]
5. A. Gordon and F. Kärtner, “Scaling of keV HHG photon yield with drive wavelength,” Opt. Express 13(8), 2941–2947 (2005). [CrossRef] [PubMed]
6. T. Popmintchev, M.-C. Chen, D. Popmintchev, P. Arpin, S. Brown, S. Alisauskas, G. Andriukaitis, T. Balciunas, O. D. Mücke, A. Pugzlys, A. Baltuska, B. Shim, S. E. Schrauth, A. Gaeta, C. Hernández-García, L. Plaja, A. Becker, A. Jaron-Becker, M. M. Murnane, and H. C. Kapteyn, “Bright coherent ultrahigh harmonics in the keV x-ray regime from mid-infrared femtosecond lasers,” Science 336(6086), 1287–1291 (2012). [CrossRef] [PubMed]
7. O. D. Mucke, S. Fang, G. Cirmi, G. M. Rossi, S.-H. Chia, H. Ye, Y. Yang, R. Mainz, C. Manzoni, P. Farinello, G. Cerullo, and F. X. Kartner, “Toward waveform nonlinear optics using multimillijoule sub-cycle waveform synthesizers,” IEEE J. Sel. Top. Quantum Electron. 21(5), 1–12 (2015). [CrossRef]
8. H. Hanieh, Fattahi, H. Wang, A. Alismail, and F. Krausz, “Towards high-power, multi-TW light transients,” in Conference on Lasers and Electro-Optics, OSA Technical Digest (2016) (Optical Society of America, 2016), paper SM1M.6.
9. H. Fattahi, A. Alismail, H. Wang, J. Brons, O. Pronin, T. Buberl, L. Vámos, G. Arisholm, A. M. Azzeer, and F. Krausz, “High-power, 1-ps, all-Yb:YAG thin-disk regenerative amplifier,” Opt. Lett. 41(6), 1126–1129 (2016). [CrossRef] [PubMed]
10. A. V. Tikhonravov, M. K. Trubetskov, and G. W. DeBell, “Optical coating design approaches based on the needle optimization technique,” Appl. Opt. 46(5), 704–710 (2007). [CrossRef] [PubMed]
11. A. V. Tikhonravov and M. K. Trubetskov, “OptiLayer software,” http://www.optilayer.com.
12. M. K. Trubetskov, V. Pervak, and A. V. Tikhonravov, “Phase optimization of dispersive mirrors based on floating constants,” Opt. Express 18(26), 27613–27618 (2010). [CrossRef] [PubMed]
13. V. Pervak, A. V. Tikhonravov, M. K. Trubetskov, S. Naumov, F. Krausz, and A. Apolonski, “1.5-octave chirped mirror for pulse compression down to sub-3 fs,” Appl. Phys. B 87(1), 5–12 (2007). [CrossRef]
14. D. Ristau, H. Ehlers, T. Gross, and M. Lappschies, “Optical broadband monitoring of conventional and ion processes,” Appl. Opt. 45(7), 1495–1501 (2006). [CrossRef] [PubMed]
15. M. Trubetskov, T. Amotchkina, A. Tikhonravov, and V. Pervak, “Reverse engineering of multilayer coatings for ultrafast laser applications,” Appl. Opt. 53(4), A114–A120 (2014). [CrossRef] [PubMed]
16. A. N. Tikhonov and V. I. Arsenin, Solutions of Ill-Posed Problems (Winston and Sons, 1977).
17. V. Pervak, “Recent development and new ideas in the field of dispersive multilayer optics,” Appl. Opt. 50(9), C55–C61 (2011). [CrossRef] [PubMed]
18. V. Pervak, M. K. Trubetskov, and A. V. Tikhonravov, “Robust synthesis of dispersive mirrors,” Opt. Express 19(3), 2371–2380 (2011). [CrossRef] [PubMed]
19. L. Gao, F. Lemarchand, and M. Lequime, “Exploitation of multiple incidences spectrometric measurements for thin film reverse engineering,” Opt. Express 20(14), 15734–15751 (2012). [CrossRef] [PubMed]
20. I. H. Malitson, “Interspecimen comparison of the refractive index of fused silica,” J. Opt. Soc. Am. 55(10), 1205–1208 (1965). [CrossRef]