Abstract
In the last decade, increasing attention has been drawn to high-precision optical experiments, which push resolution and accuracy of the measured quantities beyond their current limits. This challenge requires to place optical elements (e.g., mirrors, lenses) and to steer light beams with subnanometer precision. Existing methods for beam direction computing in resonators, e.g., iterative ray tracing or generalized ray transfer matrices, are either computationally expensive or rely on overparameterized models of optical elements. By exploiting Fermat’s principle, we develop a novel method to compute the steady-state beam configurations in resonant optical cavities formed by spherical mirrors, as a function of mirror positions and curvature radii. The proposed procedure is based on the geometric Newton method on matrix manifold, a tool with second-order convergence rate, that relies on a second-order model of the cavity optical length. As we avoid coordinates to parametrize the beam position on mirror surfaces, the computation of the second-order model does not involve the second derivatives of the parametrization. With the help of numerical tests, we show that the convergence properties of our procedure hold for non-planar polygonal cavities, and we assess the effectiveness of the geometric Newton method in determining their configurations with high degree of accuracy and negligible computational effort.
Similar content being viewed by others
References
Schreiber, K.U., Wells, J.-P.R.: Large ring lasers for rotation sensing. Rev. Sci. Instrum. 84, 041101 (2013)
Cuccato, D., et al.: Controlling the non-linear intracavity dynamics of large He–Ne laser gyroscopes. Metrologia 51, 97 (2014)
Beghi, A., et al.: Compensation of the laser parameter fluctuations in large ring-laser gyros: a Kalman filter approach. Appl. Opt. 51, 7518–7528 (2012)
Schreiber, K.U., Gebauer, A., Wells, J.-P.R.: Long-term frequency stabilization of a \(16{\rm m}^{2}\) ring laser gyroscope. Opt. Lett. 37, 1925–1927 (2012)
Bosi, F., et al.: Measuring gravitomagnetic effects by a multi-ring-laser gyroscope. Phys. Rev. D 84, 122002 (2011)
Di Virgilio, A., et al.: Performances of G-Pisa: a middle size gyrolaser. Class. Quantum Gravity 27, 084033 (2010)
Santagata, R., et al.: Optimization of the geometrical stability in square ring laser gyroscopes. Class. Quantum Gravity 32, 055013 (2015)
Belfi, J., et al.: Interferometric length metrology for the dimensional control of ultra-stable ring laser gyroscopes. Class. Quantum Gravity 31, 225003 (2014)
Di Virgilio, A., et al.: A ring lasers array for fundamental physics. C. R. Phys. 15, 866–874 (2014)
For the Romy project see: http://www.geophysik.uni-muenchen.de/ROMY/
Xingwu, Y.J.L., Meixiong, C.: Generalized ray matrices for spherical mirror reflection. Opt. Express 19, 6762–76 (2011)
Saleh, B.E.A., Carl, M.: Resonator Optics. Wiley, Hoboken (2001)
Absil, P.-A., Mahony, R., Sepulchre, R.: Optimization Algorithms on Matrix Manifolds. Princeton University Press, Princeton (2008)
Marsden, J.E., Ratiu, T.: Manifolds, Tensor Analysis, and Applications (Vol. 75), Springer Science & Business Media. Springer, Berlin (2012)
Ghatak, A.: Contemporary Optics, Springer Science & Business Media. Springer, Berlin (2012)
Ishteva, M., et al.: Differential-geometric Newton method for the best \(\text{ rank }-(R_1, R_2, R_3)\) approximation of tensors. Numer. Algorithms 51, 179–194 (2009)
Kallay, M.: A geometric Newton–Raphson strategy. Comput. Aided Geom. Des. 18, 797–803 (2001)
De Lathauwer, L., De Moor, B., Vandewalle, J.: On the best rank-1 and \(\text{ rank }-(r_1, r_2,\ldots, r_n)\) approximation of higher-order tensors. SIAM J. Matrix Anal. Appl. 21, 1324–1342 (2000)
Amat, S., Busquiera, S., Gutirrez, J.M.: Geometric constructions of iterative functions to solve nonlinear equations. J. Comput. Appl. Math. 157, 197–205 (2003)
Absil, P.-A., Ishteva, M., De Lathauwer, L.: A geometric Newton method for Oja’s vector field. J. Neural Comput. 21, 1415–1433 (2009). (MIT Press)
Absil, P.-A., Mahony, R., Trumpf, J.: An extrinsic look at the Riemannian Hessian. Geom. Sci. Inf. 8085, 361–368 (2013)
Absil, P.-A., Gallivan, K.A.: Joint diagonalization on the oblique manifold for independent component analysis. In: 2006. ICASSP 2006 Proceedings. 2006 IEEE International Conference on Acoustics, Speech and Signal Processing, Vol. 5. IEEE (2006)
Armijo, L.: Minimization of functions having Lipschitz continuous first partial derivatives. Pac. J. Math. 16, 1–3 (1966)
Nocedal, J., Wright, S.J.: Numerical Optimization. Springer, New York (2000)
Beghi, A., et al.: Shape and pose of four points in \({\mathbb{R}}^{3}\). In: 55th IEEE Conference on Decision and Control (Submitted) (2016)
Huang, W., Gallivan, K.A., Absil, P.-A.: A Broyden class of quasi-Newton methods for Riemannian optimization. SIAM J. Optim. 25, 1660–1685 (2014)
Bittencourt, T., Ferreira, O.P.: Local convergence analysis of Inexact Newton method with relative residual error tolerance under majorant condition in Riemannian Manifolds. Appl. Math. Comput. 261, 28–38 (2015)
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Sándor Zoltán Németh.
Appendix: Discussion and Convergence of Algorithm 1
Appendix: Discussion and Convergence of Algorithm 1
In Algorithm 1 at each iteration, the Newton Eq. (5) is solved for the function f, then the function \(h(x)=\left\| \text {grad}\, f(x)\right\| ^{2}\) is minimized along the computed direction. In this way, we need to compute only Hess f and grad f, avoiding the computation of Hess h, that would require to compute the third derivative of f.
Proposition 1
Algorithm 1 converges to the stationary point \(x^{*}\) of the function f with quadratic convergence rate, provided that, in a neighborhood \(\mathcal {I}(x^{*})\) of \(x^*\), \(\text {grad}\, f\ne 0\), \(\text {Hess}\, f\) is injective, and the first iterate is \(x_{0}\in {\mathcal {I}}(x^{*})\).
Proof
Let x denote a generic algorithm iterate. By hypotheses, the Newton vector \(\eta _{x}\), solution of (5), is well defined. The Riemannian gradient of h reads
By evaluating the expression \(Dh(x)[\eta _{x}]\), we get
The sequence \(\left\{ \eta _{x_{k}}\right\} \) is gradient related to \(\left\{ x_{k}\right\} \). In fact by hypothesis and (28), it holds grad\(\, h(x_{k})\ne 0\); therefore, using (29), we get \(-2\sup _{\mathcal {I}(x^{*})}h(x_{k})=\sup _{\mathcal {I}(x^{*})}Dh(x_{k})[\eta _{x_{k}}]<0\). By the smoothness of the functional Hess f and of the vector field grad f, since \({\mathcal {I}}(x^{*})\) is a compact set, we can conclude that \(\left\{ \eta _{x_{k}}\right\} \) is bounded. Hence Algorithm 1 fits in the framework of Theorem 4.3.1 and Theorem 6.3.2 [13, Chapters 4–6], stating that every accumulation point of \(\left\{ x_{k}\right\} \) is a critical point of h, so that the local quadratic convergence holds. \(\square \)
Note that the Armijo condition (7) for the function h and the direction \(\eta _{x}\) can be rewritten as
where \(y_{k}=R_{x}(\gamma _{k}\eta _{x})\), \(x=x_{k}\), \(\eta _{x}=\eta _{x_{k}}\), and k is the iteration number.
Rights and permissions
About this article
Cite this article
Cuccato, D., Saccon, A., Ortolan, A. et al. Computing Laser Beam Paths in Optical Cavities: An Approach Based on Geometric Newton Method. J Optim Theory Appl 171, 297–315 (2016). https://doi.org/10.1007/s10957-016-0981-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10957-016-0981-3