We regret that such a misleading comment [Opt. Express (2013)] has been made to our paper. First Lo states in his abstract that “However, the nonlinear Rabi model has already been rigorously proven to be undefined” to later recoil and use the contradictory statement ”(. . . ) regarding the BS model with the counter-rotating terms (. . . ) Lo and his co-authors have proven that the model is well defined only if the coupling stregth g is smaller than a critical value gc = ω/4”. While Lo focuses on the validity of the quantum optics Hamiltonians and gives a misleading assesment of our manuscript, the focus of our paper is the method to map such a set of Hamiltonians from quantum optics to photonic lattices. Our method is valid for the given class of Hamiltonians and, indeed, precaution must be exerted on the paramater ranges where those Hamiltonians are valid and where their classical simulation is feasible. These parameter ranges have to be specified in for each particular case studied. Furthermore, we gave as example the Buck-Sukumar model including counter-rotating terms which is a valid Hamiltonian for some coupling parameters.
© 2014 Optical Society of America
In a recent paper we have shown a method to simulate nonlinear Jaynes-Cummings (JC) and Rabi models in photonic lattices . Lo disputes that “the proposed classical simulation is actually not applicable to the nonlinear Rabi model and the simulation results are completely invalid” . While it is true that attention must be paid to the parameter set where the quantum models are valid and where the corresponding photonic lattices are experimentally feasible, our method to map nonlinear Rabi models into the classical propagation of light through photonic lattices is valid and we encourage the reader to consult our manuscript on this regard. Unfortunately we were not aware of the work done by Lo and his coauthors which, besides enriching our article by adding appropriate boundaries for the set of parameters where the Buck-Sukumar model including counter-rotating terms is valid, would have allowed us to cite him. Let us start by the simplified version of our Hamiltonian that is related to his comment,3]. Furthermore, setting the auxiliary function to f(n̂) = n̂ reduces the Hamiltonian to the Buck-Sukumar model including counter-rotating terms which is well defined for g < ω/4 as discussed by Lo and coauthors in .
In conclusion, the statement “Hence, the proposed classical simulation is actually not applicable to the nonlinear Rabi model and the simulation results are completely invalid” is misleading as it should not refer to the classical simulation. The mapping from the quantum models to the classical analog is correct. The quantum Rabi model with any given parameter set  and the Buck-Sukumar model with counter-rotating terms for g < ω/4  are valid as stated by Lo .
We only agree with Lo that the numerical simulation presented in Fig. 3 of our paper  should have been done for an adequate value of g in which the Buck-Sukumar Hamiltonian including counter-rotating terms is valid; for example g = 0.249ωf in Fig. 1 here. However, it is true also that for large g’s the first neighbor interaction is not valid any more for large n’s, and at least second neighbor interactions should be considered. As in every classical simulation, any particular instance of f(n̂) in our model should be studied with care in order to choose a parameter range where both the quantum model is valid and the photonic lattice is experimentally feasible.
References and links
1. B. M. Rodríguez-Lara, F. Soto-Eguibar, A. Z. Cárdenas, and H. M. Moya-Cessa, “A classical simulation of nonlinear Jaynes–Cummings and Rabi models in photonic lattices,” Opt. Express 21, 12888–128981 (2013). [CrossRef]
2. C. F. Lo, “A classical simulation of nonlinear Jaynes–Cummings and Rabi models in photonic lattices: comment,” Opt. Express (2013).
3. C. F. Lo, K. L. Liu, and K. M. Ng, “The multiquantum Jaynes–Cummings model with the counter-rotating terms,” Europhys. Lett. 42, 1–6 (1998). [CrossRef]
4. K. M. Ng, C. F. Lo, and K. L. Liu, “Exact eigenstates of the intensity-dependent Jaynes–Cummings model with the counter-rotating term,” Physica A 275, 463–474 (2000). [CrossRef]